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Abstract 


— 1  A  theoretical  model  of  a  general  production  process  is  constricted.  A  production  process 
is  egarded  as  a  network  of  jointly  operating,  interrelated  activities  which  use  system  exogenous 
inputs  of  goods  and  services  to  produce  outputs.  The  production  model  displays  explicitly  the 
intermediate  product  transfers  between  activities  and  incorporates  the  time-varying  aspects  of 
production  directly.  The  primitive  elements  which  are  taken  to  be  common  to  all  production 
processes  are  the  activity  production  functions  and  the  flows  of  products,  goods  and  services. 
To  enhance  clarity  and  rigor,  the  model  is  developed  axioms ticaily,  i.e.,  properties  on  the  primi¬ 
tive  elements  which  are  conjectured  to  be  true  in  order  to  facilitate  the  theory  are  identified. 

■  cau " 4  •(■■■  .  , 

The  general  model  extends  previous  axiomatic^  models  of  production) used  in  economic 
theory.  Specifically,  laws  of  production  and  Shephard's  Duality  Theorem  are  proved  using 
axioms  of  the  general  model.  Moreover,  the  general  model  provides  guidelines  as  to  what 
en tails  a  satisfactory  model  of  production  so  that  it  may  be^tfiizetfto  study  models  and  solu¬ 
tions  of  specific  production  planning  problems.  To  illustrate,  the  general  model  is  used  to  sys- 

.  •  c 

tematically  analyze  a  heuristic  solution  proposed  by  Leachman  and  Boysen'tTHSJr f or  the  prob¬ 
lem  of  multi-project  resource-use  planning  and  to  show  bow  their  approach  can  be  extended  and 
improved.  j'fT|  ^  £  •. 
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1.  INTRODUCTION 


In  production  planning,  a  decision-maker  identifies  objectives,  available  choices,  and  the 
result  of  each  choice.  Once  the  particular  problem  has  been  described,  a  model  is  then  formu¬ 
lated  as  a  means  of  accomplishing  the  objectives.  This  includes  modeling  the  result  of  each 
choice  available  to  the  decision-maker.  In  production  planning,  a  model  of  this  choice-result 
relationship  is,  either  indirectly  or  directly,  a  model  of  the  actual  production  process.  Thus, 
modeling  the  actual  production  process  is  a  necessary  and  essential  aspect  of  production  plan- 
ning. 

Because  the  production  planner  has  to  develop  a  model  which  is  tractable  for  analysis,  he 
often  makes  certain  assumptions  which  are  not  specifically  identified  or  adequately  justified. 
Hence,  the  development  of  a  model  for  a  particular  production  process  is  primarily  an  art.  To 
elevate  the  model-building  to  a  more  scientific  level  would  require  the  modeler  to  identify  all 
of  the  assumptions  taken,  either  implicitly  or  explicitly,  and  justify  each  one.  Therein  lies  a 
fundamental  problem  confronting  research  in  the  area  of  production  planning:  how  does  one 
know  what  were  the  implicit  and  explicit  assumptions  taken?  There  are  no  general  guidelines  to 
help  the  production  planner  identify  the  assumptions,  and  worse,  there  is  no  mechanism  to  sys¬ 
tematically  analyze  the  limitations  such  assumptions  necessarily  impose. 

Another  fundamental  problem  confronting  research  in  the  area  of  production  planning  is 
the  analysis,  or  lack  thereof,  of  heuristic  solutions.  A  rigorous  analysis  of  a  heuristic  solution 
would  be  one  which  provides  a  rational  basis  or  logical  foundation  for  the  proposed  methodol¬ 
ogy.  This  is  impossible  unless  ail  of  the  assumptions  of  the  original  model  are  identified. 

As  a  first  step  towards  elevating  the  art  of  production  modeling  and  the  analysis  of  heuris¬ 
tic  solutions  to  a  science,  we  construct  a  theoretical  model  of  a  general  production  process 
(Chapter  2).  A  production  process  is  regarded  as  a  network  of  jointly  operating,  interrelated 
production  activities  which  use  system  exogenous  inputs  of  goods  and  services  in  production  to 
produce  final  outputs.  The  production  model  displays  explicitly  the  intermediate  product 
transfers  between  activities.  We  identify  the  primitive  elements  which  we  take  to  be  common  to 


all  production  processes.  They  are  the  activity  production  functions  and  the  flows  of  goo—  and 
services  through  time.  Die  model  incorporates  ihe  time-varying  aspects  of  a  production  process 
directly.  An  explicit  description  of  how  :  ;pv'  .s  transformed  into  output  is  presented.  Exten¬ 
sions  and  limitations  of  the  present  model  are  also  discussed. 

To  enhance  clarity  and  rigor,  we  develop  the  model  axiomattcaUy.  That  is,  properties  on 
the  primitive  elements  which  we  take  to  be  true  in  order  to  facilitate  the  theory  are  identified. 
The  chief  reasons  for  axiomatizing  the  theoretical  model  are  (1)  it  discloses  many  of  the  hidden 
assumptions,  (2)  it  displays  the  structure  of  the  theory,  (3)  the  key  concepts  and  hypotheses 
are  identified,  (4)  the  consequences  of  changes  in  the  foundations  are  better  realized,  and  (S) 
the  shortcomings  of  the  theory  can  be  spotted  and  corrected.1 

Since  a  specific  model  of  a  production  process  necessarily  imposes  assumptions  on  the 
primitive  elements,  the  general  model  facilitates  the  identification  of  the  assumptions.  To  illus¬ 
trate,  we  describe  in  our  framework  the  production  models  implicit  in  the  production  planning 
techniques  of  Material  Requirements  Planning  (MRP)  and  the  ordinary  Critical  Path  Method 
(CPM).  These  descriptions  will  clearly  reveal  the  implicit  assumptions  about  the  production 
processes  made  by  production  planners  who  use  such  techniques. 

Once  the  implicit  assumptions  on  the  primitive  elements  have  been  identified,  it  is  possi¬ 
ble  to  analyze  any  proposed  heuristic  solution  offered  to  solve  a  particular  problem.  To  illus¬ 
trate,  we  use  the  general  model  as  a  tool  to  provide  a  systematic  analysis  of  a  heuristic  solution 
proposed  by  Leachman  and  Roy  sen  [1983]  for  the  problem  of  multi-project  resource-use  planning 
for  a  multi-project  production  system  (Chapter  4).  The  problem  is  to  determine  explicit 
resource  allocations  through  time  to  projects  to  insure  that  schedules  are  met.  Our  systematic 
analysis  not  only  provides  a  logical  foundation  for  their  approach  but  more  importantly  shows 
how  their  approach  can  be  extended  and  improved.  The  analyses  carried  out  in  this  chapter 
illustrate  the  value  of  using  a  general  conceptual  framework  of  a  production  system  to  evaluate 

proposed  heuristic  solutions  to  production  planning  problems. 

1  Adapted  from  Buns*  (1973) 


Theoretical  models  of  "general"  production  processes  have  been  developed  before.  Recog¬ 
nizing  the  benefits  of  axiomatizatioa,  Shephard  [1970a]  developed  an  axiomatic  description  of 
a  steady-state  production  technology.2  The  primitive  element  in  his  model  was  the  correspon¬ 
dence  which  modeled  the  input-to-output  relationship.  The  axioms  were  therefore  imposed  on 
the  correspondences.  Realizing  that  a  steady-state  framework  did  not  ^corporate  the  dynamic 
aspects  of  production  directly,  Shephard  and  Fare  [1980]  extended  the  framework  to  model 
dynamic  production  systems.  Again,  the  primitive  element  was  the  correspondence  which 
modeled  the  input-to-  output  relationship.  Functions  of  time  were  taken  to  model  the  flows  of 
goods  and  services.  The  axiomatic  description  for  this  model  was  virtually  the  same  as  the 
steady-state  model  except  that  certain  mathematical  axioms  were  employed  to  facilitate  the 
theory.  Shephard  et.  al.  [1977]  developed  the  first  network  model  of  production  and  later 
presented  an  axiomatic  description  of  this  model  in  1981. 3 

The  axiomatic  models  oi  production  constructed  in  the  past  are  not  useful  for  the 
development  or  evaluation  of  dynamic  production  planning  models.  However,  these  models 
have  proved  to  be  useful  for  the  development  of  steady-state  cost  and  production  functions.4 
More  importi  tly,  perhaps,  the  previous  axiomatic  mode’s  enable  one  to  prove  laws  of  produc¬ 
tion.  By  proving  that  laws  of  production  hold  from  the  axiomatic  description,  the  question  of 
the  validity  of  the  laws  is  reduced  to  the  question  of  whether  the  axioms  in  the  axiomatic 
framework  are  appropriate.  We  continue  this  worthwhile  task  by  proving,  in  our  general  setting, 
tv  o  variants  of  the  Law  of  Diminishing  Returns  as  formulated  by  Shephard.5  In  addition,  we  dis¬ 
cuss  technical  efficiency  and  provide  two  different  proofs  of  Shephard's  Duality  Theorem  which 
be.cer  explains  this  famous  theorem  (Chapter  3). 

Since  our  model  explicitly  defines  the  correspondence  which  models  the  input-to-ou.put 
relationship,  our  model  is  more  descriptive  of  the  production  process  than  the  past  axiomatic 

2  See  also  Shephard  11953),  11970M. 

3  See  also  Heckman  and  Shephard  (19831. 

4  See,  for  example.  Htnoch  and  Rothschild  (19721. 

5  See  Shephard  and  Fare  (1980). 


frameworks  mentioned  above.  In  addition,  since  the  choice  of  primitive  elements  dictates  the 
nature  of  the  axiomatic  description,  our  axiomatic  description  is  completely  different  from  past 
frameworks.  Some  of  the  axioms  which  were  taken  in  the  earlier  frameworks  were  completely 
analytical,  i.e,  it  would  not  be  possible  to  verify  their  validity  by  experimentation.  We  believe 
our  axiomatic  description  can  be  verified  through  experimentation  and  is  easier  to  justify. 

In  summary,  the  funeral  framework  of  production  introduced  here  provides  guidelines  as 
to  what  entails  a  satisiactory  model  of  production.  The  framework  is  "general"  enough  to  extend 
previous  axiomatic  models  useful  for  developing  steady-state  cost  and  production  functions  and 
for  understanding  laws  of  production  in  economic  theory.  Moreover,  it  is  ’general"  enough  to 
study  models  and  solutions  of  specific  production  planning  problems. 


2.  THE  GENERAL  MODEL 


In  this  chapter,  we  develop  a  general  model  of  a  production  process.  Section  2.1  presents 
the  frame*vork  of  the  model  Section  2.2  devcloos  the  model  axiomatically.  Section  2.3  illus¬ 
trates  the  generality  of  the  model.  Section  2.4  discusses  the  limitations  of  the  present  model 
and  describes  how  the  model  could  be  adapted  to  St  more  specific  cares. 

2.1.  The  Framework  of  the  Model 

In  Section  2.111,  we  present  a  conceptual  framework  common  to  all  network  models  of 
production.1  In  Section  2.1.2,  we  describe  hew  we  cnoose  to  model  the  flows  of  goods  and  ser¬ 
vices.  Sections  2.1.3  and  2.1.4  develop  the  modet  of  the  transformation  of  input  to  output  at 
the  activity  and  network  levels. 

2.1.1.  Production  Networks:  A  Conceptual  Frair-work 

A  production  system  is  modeled  as  a  directed  network,  the  nodes  of  which  represent  primi¬ 
tive  production  activities.  (An  example  of  a  Production  Network  is  shown  in  Figure  (2-1).)  Prim¬ 
itive  production  activities  tre  those  within  which  the  intermediate  product  transfers  need  not  be 
considered  for  the  purposes  at  hand.  The  nodes  are  connected  by  directed  arcs  to  indicate  pos¬ 
sible  transfers  of  intermediate  and  final  products.  (Cycles  are  permitted.)  System  exogenous 
inputs  such  as  labor  services,  machine  and  facility  services,  energy  and  fuels,  etc.,  are  treated  as 
transfers  from  cn  initial  node  A<y  For  a  system  with  N  producing  activities,  final  outputs  are 
taken  as  delivered  to  node  AN+\.  Thus,  a  production  system  is  regarded  as  a  jointly  operating, 
finite  number  of  interrelated  primitive  production  activities  A\,  A  2, . . .  ,AN  which  use  system 
exogenous  inputs  of  goodt  and  services  in  production  to  produce  final  outputs. 

Note  that  the  production  model  displays  explicitly  the  intermediate  product  transfers.  This 
display  is  essential  for  dynamic  models  of  production  since  final  output  evolves  as  the  evolu¬ 
tionary  flow  of  intermediate  products  to  finJ  products. 
s  Saa,  for  axamplt,  Shephard  at.  d.  11977),  Shephard  U9S1J,  tad  Hackman  and  Shephard  (1913). 


FIGURE  (2-1) 

EXAMPLE  OF  A  PRODUCTION  NETWORK 
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2.1.2.  Modeling  the  Flews  of  Goods  tod  Services 

The  flows  of  goods  and  services  (input  and  output)  should  have  a  truly  dynamic  character. 
A  flow,  therefore,  will  be  taken  to  be  an  element  of  an  appropriate  subset  of  the  set  of  nonne¬ 
gative  functions  defined  on  the  nonnegative  part  of  the  real  line.  Each  How  will  be  referred  to 
as  a  time-rate  history .2 

There  are  two  fundamental  types  of  flows.  The  first  and  more  common  type,  called  a  cm- 
tinuous  flow,  is  one  for  which  x(t)  represents  the  rate-quantity  per  unit  time-at  time  t.  The 
second  type,  called  an  event-based  flow,  is  one  for  which  x(t)  is  a  numerical  representation  of 
an  event  at  time  t.  One  example  of  an  event-based  flow,  suitable  for  project-oriented  produc¬ 
tion  systems  (see  Section  2.1)  is  when 

1 1  if  r  is  the  project  completion  time 
x(t)~{o  otherwise. 

Another  example,  suitable  for  batch  transfen  of  intermediate  products,  is  when  x(r)  indicates 
the  quantity  transferred  at  time  r. 

2.1.3.  The  Activity’s  Dynamic  Prodnctlon  Conesp— deace 

How  we  choose  to  model  the  relationship  of  input  into  output  at  the  activity  level  is 
described  in  this  section.  The  description  of  the  input-output  relationship  at  the  network  level 
is  presented  in  the  next  section. 

To  produce  output,  each  activity  utilizes  system  exogenous  inputs  and  intermediate  pro¬ 
ducts.  The  outputs  may  be  intermediate  products  used  as  inputs  by  other  activities,  or  final 
products,  or  mixtures  of  both  as  in  the  case  of  spare  parts. 

As  notation,  let 


1  Axioms  for  Um  Bowi  ot  foods  and  mtvkm  tn  praaaowd  m  Sactioo  (2.2.1). 
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x0,  —  Oto,,  .  . .  ,*0,)  denote  a  vector  of  n  time-rate  histories  of  system  exogenous 

inputs  allocated  to  the  /**  activity,  i— 1,2, . A\ 

V,j  «•  (Vj,  . . . ,  K")  denote  a  vector  of  ,n  time-rate  histories  of  transfers  of  outputs 
from  the  i*  activity  to  they'*  activity,  i— 1,2 . N,j— 1,2, . . .  ,A+1, 

V,  -  (K\  . . . .  K")  denote  a  vector  of  m  time-rate  histories  of  outputs  of  the  /'* 

activity.  /— 1,2 . A. 

Note  that  in  the  foregoing  representation'  constants  n  and  m  were  taken  for  convenience.  That 
is,  a  component  of  X& ,  V0  or  V,  may  be  zero.3 

Abstract'y,  the  activity’s  dynamic  production  correspondence,  denoted  by  I,,  is  taken  &s  a 
correspondence  (set- valued  mapping)  V,  — »  Lt(V,)  where  informally  one  interprets  the  state¬ 
ment  that 

L.  £  €  L'(K) 

\  J-t  I 

to  mean  that  the  /*  activity  A,  may  produce  V,  if  allocated  X&  as  system  exogenous  input  (over 

N 

time)  and  £  VJt  as  intermediate  product  transfers  from  other  activities  (over  time). 

j~i 

,  To  define  the  correspondences  L,  more  formally,  one  needs  to  differentiate  between  the 
concepts  of  allocation  and  application.  Since  an  activity  may  be  allowed  to  dispose  or  store  its 

N  | 

inputs,  what  is  allocated  to  the  activity  u  input  x&.  X  m*y  not  be  what  is  actually 

I 

applied,  or  consumed,  as  input  into  the  production  process.  As  noution,  let 

X  “  (y,1 . >,")  denote  a  vector  oi  .?  time-rate  histories  of  system  exogenous 

inputs  applied  into  the  production  process  of  A,,  /-l, 2, ....  N, 

W,  -  (H',1 . W'7")  denote  a  vector  of  m  time-rate  histories  of  intermediate  pro¬ 

duct  transfer  inputs  applied  into  the  production  process  of  A,,  '•1,2, . . . ,  N, 

,  1  Th«  tcnn ‘zero’ here  refers  to  Um  function  nidi  Uut 
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T,  -  (r,1. ....  77,  r;+1 . T?*m)  denote  a  vector  of  n-tm  time-rate  histories  of 

disposal  for  both  input  and  output  for  A,,  /— 0,1, 2 . Af,.N+l. 

I  N 

Given  a  particular  allocation  of  input  nr0,,  £  K>(  those  choices  for  applications  of  input 

l  J-i 

y, ,  H'  and  disposal  T,  which  are  feasible  are  those  for  which  the  following  inventory  balance 
constraints  are  satisfied:4 

if  the  j'h  system  exogenous  input  is  activity  storable, 

I 

0  <  cj  +  J  (4,  -  Tj-yflJn  <  Cj,  Vr  €  (2.!) 

if  the  k'h  product  is  activity  storable, 

t 

0  <  b'  +  J  ( X  K~  V+k~  WUn  <  Vr  €  *+  (2.2) 

*o  J 

if  the  j **  system  exogenous  input  is  not  activity  storable, 

0  -xL-Tf-rf  (2.3) 

if  the  k*  product  is  not  activity  storable, 

0”  Z  77+*~  W  (2.4) 

j 

where  cj,  oft  R+  represent  initial  stocks,  if  any,  and  C/,fi,*€  R+  U  {<»)  represent 
the  constant  capacity  levels.  It  is  understood  that  Tj  ( 7 *)  is  zero  if  the  /*  system 
exogenous  input  (k'A  product)  is  not  disposable. 

The  model  at  the  activity  level  assumes  the  existence  of  a  production  Junction,  denoted  by 
F„  which  takes  a  vector  of  inputs  applied  to  the  production  process  (y(,W')  into  realized  vector 

4  The  measure  M  will  be  defined  precisely  in  Axiom  t  next  section.  The  set  on  which  we  are  intefriting 

is  I0.fl. 
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(of  dimension  m)  of  outputs  of  the  activity  obtained  through  production,  F,(y,,W,).s  Since 
activities  may  be  allowed  to  dispose  or  store  products  produced  F, (y,,W,)  need  not  be  equal  to 
K-  By  the  expression  "F,  (y^W,)  is  enough  to  support  output  level  V'  we  mean  that  the  fol¬ 
lowing  inventory  balance  constraints  are  satisfied: 
if  the  klh  product  is  activity  storable. 


/ 

0  <  r*  +  F*<y, %W,)-T?+*-  V*\dp.  <  B*.  Vr  €  + , 


(2.5) 


if  the  Ar1*  product  is  not  activity  storable. 


0-  f*0'„K')-77+*~  Vf. 


(2.6) 


Thus,  to  say  that 

J- 1 


€  L,(V,)  we  mean  that  there  exist  an  application  vector 


(k.H',)  and  disposal  vector  T,  such  that  (2.1M2.6)  are  satisfied.  We  now  turn  to  describing 
the  input-output  transformation  at  the  network  level. 


2.1.4.  The  Network  Dynamic  Production  Correspondence 

Let  u  —  (u*. . . .  ,um)  denote  a  vector  of  final  output  time-rate  histories.  The  network 
dynamic  production  correspondence,  denoted  by  LN,  is  a  correspondence  u  —  LN(u)  which, 
loosely  described,  is  the  set  of  all  vectors  of  system  exogenous  input  rate  histories  x  that  when 
appropriately  allocated  to  the  activities  may  produce  the  vector  u  of  final  output  rate  histories. 

To  be  specific,  x  €  LN(w)  means  that  we  can  find  allocations  of  system  exogenous  inputs 
to  the  activities,  x0,'5,  and  allocations  of  intermediate  products  to  activities  obtained  from  other 
activities,  the  VtJ  ’s ,  such  that  the  following  inventory  balance  constraints  are  satisfied: 

On  the  input  side, 

if  the  J *  system  exogenous  input  is  system  storable, 

*  Axioms  for  the  Activity  Production  Function!  ut  pram iMd  ia  Section  (2.2.2). 
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l 

0  <  V ■+  I  {*>-  n- 2,xi)dn  <  #  ,  Vr  €  A+ 

Wft  I 


(2.7) 


if  the  j,h  system  exogenous  input  is  not  system  storable. 


(2.8) 


where  b>  €  is  the  initial  stock,  if  any,  and  BJ  €  *+  U  M  represents  the  constant 
capacity.  It  is  understood  that  Ti  -0  if  the  j*  system  exogenous  input  is  not  sys¬ 
tem  disposable. 

On  the  output  side, 

if  the  k*  product  is  system  storable, 

I 

0  <  +  JlX  Vb+ 1  -  nt\  ~«kW <  Bfi+t ,  Vr  €  J?4  ,  (2.9) 

•'o  * 

if  the  product  is  not  system  disposable, 

0-XK£*+i  --Tfft-u*  (2-10) 

’  i 

where  b£+i  €  is  the  initial  stock,  if  any,  and  fljv+)  €  /?+U  {»)  represents  the 
constant  capacity.  It  is  understood  that  7jy+{  -0  if  the  k1*  product  is  not  system 
disposable. 

F.nally,  one  needs  to  insure  that  the  individual  activities  can  produce  what  is  required  of  them, 
so  one  adds 

*Wi>.  (2.11) 

Two  important  comments  are  in  ordet  with  respect  to  the  definition  of  the  correspondence 
LN.  First,  we  have  tacitly  assumed,  for  each  /, 


2 


k  -  I  y,j .  (2.i2) 

j- 1 

This  convenient  assumption  is  not  restrictive  because: 

(1)  Storage  and  disposal  of  intermediate  products  are  already  allowed  at  the  activity 
level, 

(2)  If  the  produr  ion  system  is  such  that  it  is  more  appropriate  to  dispose  or  store 
interinediate  products  at  the  system  levei-for  example,  products  are 
warehoused-then  the  acts  of  disposal  and  storage  may  be  simply  modeled  as  a 
separate  activity.  If  the  acts  of  storing  and  disposing  were  costly,'  then  model¬ 
ing  such  acts  as  activities  would  be  appropriate. 

N 

The  second  comment  concerns  the  expression  £  VJt  (the  vector  sum  of  intermediate  product 

J-  > 

transfer  inputs  from  ail  activities  into  activity  /).  It  is  dear  from  this  expression  that  our  model 
does  not  incorporate  transfer  or  shipment  lags.  To  incorporate  this  time  lag,  one  could  define 
Vj,  ( t )  as  what  Aj  sends  A,  at  time  t  and  introduce  Kj’(r)  as  what  A,  receives  from  Aj  at  time  t. 
For  example,  if  there  were  a  constant  time  lag  lJt  for  shipment  then 

K»-l  o  jro««,. 

Finally,  we  make  a  useful  definition.  If  x  €  LN(u),  there  are  many  possible  collections 
of  flows  xo/s,  V,j's,  y,'s  W,'s,  and  T/s  of  goods  and  services  which  satisfy  (2.1)r(2.11).  By 
the  expression  ’a  feasible  flow  for  input  x  to  support  output  level  u\  we  mean  one  such  collec¬ 
tion  of  flows.  On  occasion  we  will  simply  say  *a  feasible  flow  to  support  output  level  u*  if 
reference  to  a  specific  x  is  not  required  and  a  ’feasible  flow”  if  reference  to  a  spetific  u  is  not 
required. 
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2.2.  An  Axiomatic  Presentation  of  the  General  Model 

!n  Section  2.1,  we  provided  the  framework  of  the  general  model.  We  first  introduced  pro¬ 
duction  networks  as  a  conceptual  model.  Then  we  modeled  the  flows  of  goods  and  services  so 
that  they  would  be  truly  dynamic.  Finally,  we  described  how  input  is  transformed  into  output. 

From  the  discussion  in  Section  2.1,  it  is  clear  that  two  types  of  primitive  elements  comprise 
the  model: 

(1)  the  flows  of  goods  and  services,  and 

(2)  the  activity  production  functions. 

The  purpose  of  this  section  is  to  define  precisely  these  elements.  As  mentioned  in  the  intro¬ 
duction,  the  best  way  of  being  precise  is  through  the  axiomatic  method.  In  Section  2.2.1,  we 
develop  and  justify  the  axioms  taken  for  the  flows  of  goods  and  services.  In  Section  2.2.2,  we 
develop  and  justify  the  axioms  taken  for  the  activity  production  functions. 

We  remind  the  reader  that  it  is  our  attitude  that  an  axiom  is  not  an  a  priori  truth  but 
rather  a  scientific  hypothesis  conjectured  in  order  to  facilitate  a  theory.  The  justification  of  each 
axiom  rests  mainly  on  the  assertion  that  the  property  imposed  is  a  property  one  would  expect  to 
observe  in  the  future  or  one  has  observed  in  the  past. 

2.2.1.  Axioms  for  the  Flows  of  Goods  and  Services 

2.2.1. 1.  Axiom  1:  The  Underlying  Space  of  Flows  of  Goods  and  Services 

It  was  argued  in  Section  2.1  that  to  represent  the  truly  dynamic  character  of  production,  a 
flow  of  a  good  or  service  should  be  modeled  as  a  nonnegative  function  of  time.  That  is,  if  z 
represents  a  flow  then 

z:  /?+ 

where  /?+— {r:  t  >  0}  models  the  time  axis. 
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The  mathematical  operation  of  integration  was  employed  to  define  the  inventory  balance 
constraints  (2.1)-(2.11)  presented  in  Section  2.1.  To  formally  integrate  a  flow,  it  requires  one 
to  first  define  a  measurable  space  with  respect  to  which  the  flow  is  measurable.  If  B(R+) 
denotes  the  Lebesgue  cr-field  of  R  restricted  to  R+,  then  we  take  C R+,2f  (/?+))  to  be,  the 
measurable  space.  Second,  one  must  select  an  appropriate  measure  p  on  B  (R+). 

The  measure  p  must  account  for  both  continuous  and  event-based  flows.  Since  a  con¬ 
tinuous  flow  measures  the  quantity  per  unit  time  continuously  through  time,  Lebesgue  meas¬ 
ure,  denoted  by  X,  is  a  suitable  measure  which  may  be  used  to  integrate  a  continuous  flow. 
However,  Lebesgue  measure  is  not  suitable  for  event-based  flows.  Events,  as  we  choose  to 
think  of  them,  do  not  occur  continuously  through  time.  That  is,  the  set  on  which  an  event- 
based  flow  is  positive  has  Lebesgue  measure  equal  to  zero.  A  suitable  measure  which  may  be 
used  to  integrate  an  event-based  flow  is  a  counting  measure,  denoted  by  of  the  following  type: 

A  countably  infinite  index  set  T  —  {f*)“  C  BtR+)  u  assumed  to  exist  for  which  if 
B€B(R+), 


j  \b  n  r\  if  Is  n  rl  <  • 

v(B)  —  j  oo  otherwise 
where  \A  I  denotes  the  cardinality  of  the  set  A . 

Each  point  in  T  indicates  a  possible  time  of  ah  event  So,  to  account  for  both  continuous  and 
event-based  flows  the  measure  p  is  taken  to  be  of  the  form  X  +  r.  We' are  now  ready  to  state 
Axiom  1. 

Axiom  1 

If  z  represents  a  flow  of  a  good  or  service,  then  z  €  L?m+,B(Jl+) ,X+v)  where  L? 
denotes  the  nonnegative  orthant  of  Lm.  If  T  denotes  the  index  set  associated  with  v,  then  T  is 
assumed  to  satisfy  the  following  properties: 
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(1)  0-/0<'l<'2<  •••  <4,  <*„+!<  ••• 

(2)  inf[tk-tk-i}>  0 

k  >  1 

(3)  sup  {tk  -/*-/}<  0*  . 

*  >  1 

We  make  3  comments  about  Axiom  1: 

(A)  The  restrictions  on  T  were  imposed  to  allow  us  to  view  lk  —  (r*_j,  tk)  as  a  period  in  pro¬ 
duction  planning,  the  points  in  T  as  time-grid  points,  and  T  as  the  time  grid.  The  assump¬ 
tion  of  a  time  grid  T  satisfying  properties  (l)-(3)  above  is  always  assumed  in  discrete¬ 
time  planning  and  control. 

(B)  Since  a  flow  is  defined  up  to  sets  of  measure  zero,  it  is  understood  that  a  continuous  flow 
is  constrained  to  be  one  for  which  2(f) -0  if  t  €  T  and  an  event-based  flow  is  one  for 
which  r(f)  —  0  if  /  i  T. 

(C)  Axiom  1  states  that  a  flow  is  initially  constrained  to  be  a  function  which  is  (i)  nonnega¬ 
tive,  (ii)  measurable,  and  (iii)  essentially  bounded.  The  set  of  functions  satisfying  these 
three  properties  is  large.  It  includes  functions  which  are  not  flows  of  goods  and  services 
on.  expects  to  occur  in  production.  Therefore,  additional  constraints  in  the  following 
axioms  r.eed  to  be  imposed  to  define  what  constituted  an  acceptable  flow.  The  constraints 
to  be  imposed  will  apply  tc  each  flow  type.  By  the  expression  "flow  type”,  we  refer  to  a 
particular  class  of  flows.  For  example,  one  of  the  flow  types  is  the  class  of  flows 
corresponding  to  the  allocation  of  the  j *  input  to  the  i*  activity  which  we  have  denoted 
in  Section  2.1  by  the  symbol  x^,. 

2.2.1. 2.  Axiom  2:  Limiting  the  Shape  of  the  Flows  of  Goods  and  Services 

Each  flow  is  (essentially)  bounded  (in  norm)  over  the  infinite  horizon  and  hence  in  each 
period.  We  further  insist  that  a  flow’s  bound  in  a  period  is  bounded  by  a  function  of  the  cumu¬ 
lative  amount  of  the  flow  in  that  period. 


The  thrust  of  Axiom  2  is  to  limit  flows  which  are  sharply  'peaked.*  That  is.  Axiom  2  in 
effect  limits  the  deviation  between  the  bound  of  a  flow  in  a  period,  llz-l^tL,  and  the  mean 


l 


zdfi 


value  of  the  flow  in  the  period,  — rrr-  An  example  of  a  flow  which  is  sharply  peaked  is 
graphically  depicted  in  Figure  (2-2). 

The  restriction  to  be  imposed  in  Axiom  2  is  made  parametrically.  The  parameters  are  not 
constants  but  continuous  non-decreasing  nonnegative  scalar  valued  functions  on  R+. 


Axiom  2 

For  the  /**  flow  type,  there  exists  families  of  parameters  [gflZ- lt  (A*}*Li,  so  that  if  z  is  a 
flow  of  this  type  then  for  all  k 

(1)  if /a//*  <  A  then  llr-l/JL  <  gj(A) 

‘k 

(2)  if  f  zdn  <  A  then  llz-l(.  }ll«.  <  h'k(A)x 

In  discrete-time  planning  and  control,  the  assumption  almost  always  taken  is  that  each 
flow  is  a  continuous  flow  which  is  a  step-junction  associated  with  the  time  grid.  That  is,  a  flow  is 
assumed  to  be  constant  on  each  period  Ik.  Since  our  time  grid  T  does  not  change,  we  will  refer 
to  these  functions  as  step-functions.  The  justification  for  the  structure  of  the  property  imposed 
in  Axiom  2  is  that  if  we  properly  select  the  parameters  {«(}“,  {Aj}"  then  Axiom  2  may  be  used 
to  limit  a  flow  to  the  set  of  step-functions.  Moreover,  we  may  seiect  parameters  in  such  a  way 
so  as  to  restrict  the  range  of  a  flow. 

Proposition  (24) 

1  Of  course,  this  reduces  to  the  statement  that  x(tt)  <  A  —  x(ik)  <  k>L4). 
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Let  S  be  a  non-empty  closed  subset  of  /?+.  Tb'tn  there  exists  parameters  {«*}”,  {A*)" 
such  that  if  z  is  a  ilow  satisfying  Axiom  2  (with  these  parameters)  then  2  must  be  a  step  func¬ 
tion  whose  range  ;s  contained  in  5. 

Proof  of  Proposition  (2.1) 

Since  Sc  is  open  in  R ,  write  it  as  a  countable,  disjoint  union  of  intervals  /„  -  (a„,b„), 

i.e.,  Sc«*  (J  I„.2 
«- 1 

For  each  k,  define  gk  and  hk  as  follows: 

tuj  (S 

(rib'4') 

a„  +  - — - —  if  X(/*)a)l  <  A  <k(Jk)b„  ,  am  >  0,  b„  <  *» 

I  b„  —  a„ 

gtCA )-]  (  a  r  ai3) 

I  aTaT  I 

- — r — if  X(/t)a„  <A<  X(/t)6„ ,  a,<0,  b„  <00 
b, , 

V4  ifX(/*)a,  <A  <X(/4)f,,*,-oo. 

! 

hkU)-'AA.  (114) 

Figure  (2-3)  exhibits  a  pictorial  representation  of  the  graph  a'  a  typical  gk.  Hire,  the  set 
5— (2,4,6).  Inspection  of  (2.13)  and  (2.14>  shows  that  each  gk  and  A*  are  accepts  ale  parame¬ 
ters. 

Since  the  nonnegative  a*'j  and  b„’s  necessarily  lie  in  5,  it  is  a  simple  matter  to  check  that 

\a  €i?+|  fo(X(A)i4)<4-/f+  (2.15) 


1  Royden  (1968].  p.  39. 
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(2.16) 


Let  2  be  a  flow  satisfying  (1)  and  (2)  of  Axiom  2.  By  (2)  and  our  choice  of  hk's,  we 
have  that,  for  each  k, 

z(tt)  <  A* (r (r*))  -  VizG*) 

which  implies  that  z (tk) - 0  for  each  k.  Hence,  z  must  be  a  continuous  flow. 

By  (1)  and  (2.1S),  we  also  have  that,  for  each  k,  z  must  satisfy  the  following  inequalities: 


<  llz  l,tILA(/*) 


<  tk ( f  zdn)kOk)  ty  (!) 

<  f  zdfi  by  12.15). 

Hence,  the  above  inequalities  must  be  equalities  which  implies  that 


f  *dli 

t*  (JT  zdfi)  m  ”jy“y  for  each*,  and 


zU)  — 


f  zdn 
x(/») 


if  f  €  /*  . 


(2.17) 


(2.18) 


In  view  of  (2.16)  and  (2.17),  it  is  immediate  that  the  constant 


r 

w 


iS. 


By  (2.18),  the  result  follows.* 
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2.2.I.3.  Axiom  3:  Existence  of  a  Lower  Bound 

Axiom  3  essentially  states  that  If  a  flow  is  positive,  then  it  must  be  some  minimal  value. 

Axiom  3 

For  the  /*  flow  type,  there  exists  an  «,  >  0  such  that  if  2  is  a  flow  of  this  type  then 

Ml0<2<€,}-0. 


In  production  planning  problems,  usually  only  the  assumption  of  nonnegativity  on  the 
flows  is  imposed.  This  assumption  is  made  in  order  to  facilitate  the  implementation  of  algo¬ 
rithms  used  to  solve  the  formulation  of  the  problem.  We  believe,  however,  that  flows  one 
would  observe  or  have  observed  in  production  satisfy  Axiom  3. 

2. 2. 1.4.  Axiom  4:  Limitation  on  Set-op  rimes 

Given  a  How  z,  a  set-up  time  for  z,  loosely  worded,  is  a  time  r  at  which  z  is  either  "start¬ 
ing  up  again"  or  "stopping."  If  z  €  I*  IR*, £(/?+), m).  then  we  define  the  set  of  set-up  times 
for  z  to  be 

S(z)  -  |t  €  1  if  I  is  an  open  interval,  r  €/  C  then 

M(/n  (z>0»  >  0  ,  and mI/H  (z-0)}  >  oj.  (2.19) 

Axiom  4  insists  that  the  time  between  any  two  set-up  times  for  a  flow  is  at  least  some  minimal 
value. 

Axiom  4 

For  the  i*  flow  type,  there  exists  a  6,  >  0  so  that  if  z  is  a  flow  of  this  type  then 
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Before  we  proceed  to  present  the  axioms  taken  for  the  activity  production  functions,  it 
will  be  useful  to  introduce  some  notation.  We  let  L?  (jt)  denote  L?  Ut+tBLR+),ii).  We  let 
the  generic  symbol  L,  denote  that  subset  of  L“  <ii)  which  satisfies  Axioms  2-4  for  flow  type  /. 
When  referring  to  a  particular  flow  type,  we  will  use  the  symbol  wnich  refers  to  flows  of  that 
type  (for  example,  L^). 

i 

2.2.2.  Axioms  far  the  Activity  Protection  F auctions 

22.21.  Axiom  S:  Cl  os  art  of  the  Domain 

Let  L,  —  L. i  x  •  •  •  x  L„.,  Lw  “  L~i  x  •  •  •  x  L_,„.  If  D,  denotes  the  domain  of 
the  function  F,,  then  D,  C  L,  x  L*.  If  there  are  additional  constraints  linking  the  domains  of 
the  inputs  applied  in  production  then  D,  may  be  a  proper  subset  of  L,  x  L*  .3  For  example,  in 
the  Leontief-type  input-output  models  of  production  the  inputs  applied  into  production,  if  non 
zero,  are  assumed  to  be  proportional.  Another  example  is  when  an  activity  utilizes  one 
machine  to  produce  several  similar  types  of  products.  In  this  example,  the  exogenous  inputs 
applied  into  production  include  the  rate  of  machine  hours  applied  to  each  type  of  product.  The 
functions  which  define  such  rates  are  linked  in  that  no  two  of  them  are  positive  at  the  same 
time. 

Axiom  5 

*  Production  ijrsums  sain  •bora  coaxrabits  tra  knpoood  on  tbs  rapicitiom  of  (roots  botonsiro  to 
acUvNM*.  w«  cfaoow  oot  to  tocorporato  ms  footura  la  our  ptoarei  modal  Such  am  will  hsvo  to 
bo  raominod  oo  an  individual  bom. 
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Endow  L*  with  the  relative  weak-star  topology  for  each  k.4  It  is  assumed  that  D,  is  closed 
in  the  product  topology  on  L,,  x  Lw  . 

It  is  necessary  for  mathematical  reasons  that  the  domain  for  each  production  function  is 
closed.  It  is  not  possible  to  examine  each  instance  of  additional  constraints  linking  the  domains 
of  the  inputs  applied  in  production  and  show  closure  in  each  case.  Hence,  Axiom  S  is  imposed. 
However,  it  will  be  shown  that  if  D,  -  Ly  x  L„  then  Dt  is  closed  in  the  weak-star  product 
topology  on  (Proposition  (3.4.2),  Appendix,  Ch.  2).  Furthermore,  the  domains 

associated  with  the  two  examples  given  above  are  also  shown  to  satisfy  Axiom  S  (Proposition 
(3.4.6),  Appendix,  Ch.  3).  The  purpose  of  examining  these  examples  and  showing  that  Axiom 
5  holds  in  each  case  is  to  give  a  plausible  basis  for  accepting  Axiom  S. 

1.2.1.2.  Axiom  6:  Nall  Activities  are  Excluded 

An  obvious  property  to  impose  on  each  production  function  is  to  insure  that  each  activity 
is  capable  of  producing  at  least  one  product. 

j 

Axiom  6 

For  all  Ait  F,  on  D,  is  non  trivial.  That  is,  3k  and  a  €  D,  s.t.  llf?(yj,H'j)IL  >  0. 

2.2.2.3.  Axiom  7:  Existence  a  T  System  Essential  Inpats 

A  subset  of  system  exogenous  inputs  is  called  essential  if  when  none  of  these  inputs  are 
applied  in  production  then  no  oitput  of  any  kind  is  possible.  Typically,  in  most  production  sys¬ 
tems  labor  services  is  an  example  of  a  set  of  esMntial  inputs. 

4  A  basic  open  natjhbofhood  K  of  taro  in  this  topotofy  has  ths  fotowms  form:  3k1 . />*€£'  tad 

•a  s  >  0  such  that  N-{t  t  A"' I  |J  <«,  V/).  N  Is  usually  denoted  by  A(y' . #";«>. 


24 


Axiom  7 

There  exists  a  set  £,  0* £  C  {1,2, . . .  ,n},  |£|<n,  such  that  for  each  <4,  if 
O’, , W', )  €  D,  with y,  -0  fory  €£  then  F,(ylt H'J-O. 

Any  V.  Mtisfying  the  above  property  will  be  called  an  essential  set  of  system  exogenous 
inputs. 


2.2.2.4.  Adorn  S:  Boundedness  of  the  Production  Function 

If  the  rates  of  application  associated  with  an  essential  set  of  system  exogenous  inputs 
applied  into  an  activity's  production  process  were  bounded,  then  we  maintain  that  ti,r  output 
rates  realized  through  production  would  be  bounded  regardless  of  the  magnitude  of  the  rates  of 
application  of  the  intermediate  product  input.  This  axiom  suggests  itself  from  experience. 


Axiom  8 

Let  £  be  an  essential  set  of  system  exogenous  inputs.  For  C  €  A+,  let 

Fc-0'€L,(|  nuutliyiL^C}. 

Then  for  each  A, 


2.2.2.S.  Axiom  t:  Input-Output  Continuity 

Let  O’i.B'i).  O':,**' j)  be  two  input  vectors  belonging  to  the  domain  D,  of  F,.  Let  A  >  0. 
Fix  the  horizon  to  the  finite  interval  [0,A].  If  the  difference  in  the  cumulative  amount  of  each 
input  in  each  period  between  the  two  input  vectors  is  sufficiently  small,  then  we  maintain  that 
the  difference  in  the  cumulative  amount  of  each  output  over  the  horizon  between  the  realized 
output  vectors  of  production  must  necessarily  be  small. 


Axiom  9 


For  each  A ,  and  V/i  €  /?+,  if  O'i.W'j)  €  D, ,  then  V«  >  0,  38  >  0  so  that  if  (y2,  W2)  €  D, 
satisfies 

(i)  max  |  J  (H  -yi)dn  |  <  8,  1 

'  7(nio,*! 

(ii)  max  |  J  -  W^dp.  |  <  8 ,  1  <  *  <  m  , 

1  J/nlo.Al 

then 

I  f  I  <  «, 

[S.A1 


There  are  examples  of  production  functions  which  do  not  satisfy  Axiom  9.  One  example 
is  the  production  function  associated  with  a  chemical  process  which  produces  a  successful  reac¬ 
tion  only  when  a  certain  level  or  critical  threshold  of  input  is  reached.  To  incorporate  these 
activity  production  functions  into  an  framework  one  may  either  choose  to  modify  the  domain 
D,  so  that  levels  below  the  critical  threshold  no  longer  belong  to  the  domain  or  modify  the  pro- 
d  ction  function  itself.  We  believe  that  either  modification  can  be  made  without  seriously 
alTecting  the  model  of  the  activity’s  production  process. 

2.2.2. 6.  Axiom  10:  Efficiency  of  the  Proportion  Fanction 

If  the  vector  of  intermediate  product  inputs  applied  into  an  activity’s  production  process  is 
fixed,  then  we  maintain  that  there  exists  a  bound  on  the  cumulative  amount  of  exogenous 
inputs  applied  in  production  before  the  activity  operates  inefficiently .  By  the  expression  "activity 
operates  inefficiently*  we  mean  that  an  activity  could  produce  at  least  as  much  cumulative  out¬ 
put  while  using  "less*  exogenous  input. 

The  bound  on  the  cumulative  amount  of  exogenous  inputs  applied  into  production  is 
assumed  to  be  a  function  of  the  vector  of  intermediate  products  applied  into  production.  It  is 
also  assumed  that  this  function  itself  is  bounded  in  a  manner  similar  to  Axiom  8. 
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Axiom  10 

For  each  A,  there  exists  a  g, :  Lw  x  R+  —  /?+  such  that  VA  €  R+  if  (y,,IF)  *  L)  with 


max  f  g,(W%h) 

(fol 

then  there  exists  a  y2  <  such  that  (y2,  W)  €  D,  and  such  that 
/  F?(y2,W)dv.  >  f  F*(yltW)dv,  1 

ft*)  (S.A) 

Let  H'c  -•  { Ff'  €  Lh-  |  max  II  W'*|l„  <  C).  Then  each  g;  is  assumed  to  satisfy  the  following 
property:  VA  € 


2.2. 2.7.  Axiom  11:  Past  Production  Not  Affected  by  Fatal*  Inputs 

Let  (yj.FF),  (yltW)  be  two  input  vectors  in  the  domain  D,  such  that  for  some  A  >  0. 


That  is,  is  the  restriction  of  y\  to  the  horizon  (0,Al.  Since  no  future  input  can  affect  the  out¬ 
put  already  generated  up  to  time  A,  then  we  maintain  that 

l<k*m 

even  though  mathematically  (yi,W)*i  (yu IP). 


Axiom  11 

For  each  A,  and  VA  €  R+  if  (y,W )  €  D,  then 


2.2. 2.8.  Axiom  12:  Cumulative  Production  Limited  in  Finite  Horizon 

We  maintain  that  in  a  finite  amount  of  time,  an  activity  cafinot  produce  an  infinite 
amount  of  output  regardless  of  the  inputs  applied. 

Axiom  12 

For  each  A,  and  VA  € 

sup  f  Fj‘{y,W)dft  <  oo,  1  <  k  <  m  . 
fbi 

We  make  an  important  observation  about  the  implication  of  Axiom  12.  If  in  a  finite 
amount  of  time  only  a  finite  amount  of  cumulative  output  may  be  realized  through  production, 
then  the  cumulative  amounts  of  the  allocations  and  applications  of  intermediate  product  input 
associated  with  any  feasible  flow  must  be  finite  in  a  finite  horizon.  By  Axiom  2,  this  implies 
that  in  a  finite  horizon  the  flows  of  the  allocations  and  applications  of  intermediate  product 
input  associated  with  any  feasible  flow  are  uniformly  bounded  in  norm.  Thus,  a  function 
B:  R+  —  R+  exists  such  that  B(h)  is  such  a  bound.  The  function  B(h)  is  a  function  of  the 
horizon  10, A 1  and  the  parameters  used  to  define  Axioms  2  and  12. 
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2.3.  Special  Cases  of  the  General  Model 

In  this  section,  we  illustrate  the  model's  generality.  In  Section  2.3.1,  we  use  the  descrip¬ 
tion  of  the  general  model  to  define  precisely  the  Dynamic  Linear  Activity  Analysis  Model  of 
Production.  In  Section  2.3.2,  we  show  how  the  Traveling  Salesman  Problem  may  be  embedded 
into  the  framework  of  the  general  model.  In  Section  2.3.3,  we  discuss  Material  Requirements 
Planning.  Finally,  in  Section  2.3.4,  we  discuss  single-project  production  systems. 

2.3.1.  The  Dynamic  Linear  Activity  Analysis  Model  (DLAAM) 

The  Dynamic  Linear  Activity  Analysis  Model  of  Production  (DLAAM)  is  ah  extension  of  the 
dynamic  Leontief  input-output  models  of  production.1  All  of  these  models  assume  a  particular 
form  for  the  domain  of  the  activity  production  functions  and  the  production  functions  them¬ 
selves. 

The  Activity  Production  Fnnctions 

Each  domain  D,  is  assumed  to  have  the  following  form:  there  exist  constants  a,,, 
>  —  1,2 . /i,a*,*  — 1,2 . msuchthat 

A  -  I  3z,  such  that  y> -  a^z, ,  Wf  -  «*z,} . 

In  other  words,  if  the  applications  of  inputs  are  positive,  then  they  must  be  proportional.,  thus, 
they  may  be  indexed  in  terms  of  one  profile,  z,  ,  called  the  intensity  curve.  The  intensity  curve  z, 
is  assumed  to  be  a  bounded  step-function  associated  with  an  equal-length  period  time  grid. 

The  production  function  F,  is  a  function  of  In  view  of  the  form  of  the  domain 

A  we  may  write  the  production  function  as  a  function  of  the  intensity  curve  denoted  by  Ft  (z,). 

Each  F,  (z,)  is  assumed  to  have  the  following  form;  there  exists  constants  c*,  k  •» 1,2 . m 

such  that 

1  First  developed  in  Shepherd,  Al-Ayet,  Leachman  119771  and  later  adapted  by  Leerhman  {19121.  See 
also  Leontief  119511,  Koopmant  119511,  and  Morsenitera  119541. 


In  other  words,  the  outputs  generated  are  assumed  proportional  (if  positive)  and  indexed  by  the 
same  profile  as  the  one  used  to  index  the  inputs. 

Inventory  Calculations 

In  the  general  model  presented  in  Chapter  2,  the  inventory  of  a  good  or  service  was  con¬ 
strained  to  be,  nonnegative  and  less  then  capacity  at  all  points  in  time.  DLAAM  relaxes  this 
restriction  in  that  the  calculations  of  inventory  only  occur  at  the  time  grid  points.  For  a  finite 
horizon,  this  relaxation  implies  that  the  level  sets  LN(u)  may  be  described  in  terms  of  a  finite 
set  of  linear  inequalities.  (The  cumulative  amounts  of  each  flow  in  each  period  become  the 
variables.)  Thus,  formulations  of  production  planning  problems  which  assume  DLAAM  as  a 
model  of  production  permit  the  use  of  linear  programming. 

The  general  inventory  balance  constraints  imposed  on  the  intermediate  product  transfers 
are: 

f  x  y^r-tx)dn  >  f  X  »  V.,  >  f, ,  1  <  *  <  m  .  (2.20) 

Ui,tr]  i  if |,  1,1  < 

A  lag  of  one  period  is  incorporated  in  (2.20)  to  insure  feasibility  since  the  calculations  of  the 
inventory  only  occur  at  the  time  grid  points.  Note  that  DLAAM  does  not  explicitly  mode!  each 

V,f  but  only  the  sums  K*X  and  ^/“X  K/<-  Th*s  is  because  DLAAM  does  not  impose 

<  j 

any  side  constraints  on  the  V,j's.  Without  side  constraints  on  the  F)/s,  these  variables  become 
redundant  In  fact,  most  production  models  do  hot  explicitly  mention  the  F/s  precisely 
because  they  do  not  impose  side  constraints  on  the  Fj/s.  For  a  further  discussion  of  DLAAM 
and  its  relationship  to  Leontief-type  input-output  models  of  production  ste  Leachman  [1984]. 


2.3.2.  The  Traveling  Salesman  Problem 

Let  H  be  a  directed  network  with  N  nodes,  N  >  3.  The  traveling  salesman  problem  is  to 
find  a  minimum-length  Hamiltonian  cycle,  i.e.,  a  cycle  passing  through  each  node  exactly  once. 
As  notation,  let  d^jy  —  length  of  arc  (ij).  We  assume  d(tJ)  is  a  positive  integer.  We  let  A 
denote  the  set  of  arcs  in  H. 

In  this  section,  we  will  construct  a  production  system  such  that,  for  an  appropriate,  choice 
of  the  final  output  vector  u,  each  feasible  flow  for  u  corresponds  to  a  Hamiltonian  cycle  and 
vice-versa.  Moreover,  what  is  being  produced  by  this  production  system  is  the  length  of  the  Hamil¬ 
tonian  cycle.  The  development  also  illustrates  the  very  general  ways  in  which  the  primitive  ele¬ 
ments  can  be  defined. 

We  first  proceed  to  construct  the  production  network  G  associated  with  the  production 
system.  The  nodes  in  G  correspond  to  the  arcs  in  H\  that  is,  if  arc  ( ij )  exists  in  H  then  node 
A  aj)  exists  in  G.  We  also  add  a  "sink"  node  A  w+uv+1).  There  is  an  arc  from  A (,j)  to  A  (kJ)  if 
(ij)  (*,/)  and  j-k.  /.i.'o,  for  each  (ij)  s*  (N+lyN+l)  an  arc  from  A^jy  to  ^(jv+ijv+i)  *s 
added. 

To  define  the  production  system  associated  with  G,  we  need  to  define  the  limitations  on 
the  flow  types.  We  list  these  below. 

(1)  The  time  grid  T  is  the  set  of  whole  numbers  {0,1,2,  •••}.' 

(2)  All  flows  are  event-based. 

(3)  the  range  of  f-ach  flow  belongs  to  {0,1}. 

(4)  Except  at  the  sink  node,  no  inventories  of  any  kind  are  allowed. 

(5)  No  disposal  of  any  kind  is  allpwed. 

(6)  There  is  only  one  exogenous  input  and  only  one  product. 

(7)  Attention  is  restricted  to  the  finite  horizon  [0,A]  where  h  «•  £  d\,j\. 

(ijHA 

The  additional  constraints  imposed  on  the  flows -other  dian  those  constraints  already 
imposed  by  the  general  model-are  listed  below: 

ff  X  Ju.nW  - 1  (2.21) 

•foj  l<*,/>*</V+t,#'*l>  I 


where  T(W'<*i/))-*min{n  |  H^(4  /)(n)-  1],  A  routine  check  shows  that  all  flow  types  satisfy  the 
axiomatic  system  presented  in  Section  3.2. 

Before  we  state  and  prove  Proposition  (2.2),  it  will  be  useful  for  the  reader  to  interpret 
the  statement 


“  1 *  (ij)SA  ,  (J,k)  €  A 

to  mean  that  *at  time  t  we  have  visited  node  j  from  node  /  and  will  now  visit  node  k  from 
node./." 


Proposition  (2.2) 

Let  u  —  l)*).  Each  feasible  flow  to  support  output  u  for  the  production  system 
corresponds  to  a  Hamiltonian  cycle  in  H.  Conversely,  to  each  Hamiltonian  cycle  in  H  we  may 
associate  a  feasible  flow  to  support  output  level  u . 


Proof  of  Proposition  (2.2) 

We  first  show  that  a  feasible  flow  corresponds  to  a  Hamiltonian  cycle.  Assume  a  feasible 
flow  to  support  u  exists  and  select  one.  By  (2.21)  there  exists  exactly  one  (*,/)  €  A  such  that 

y{kJ)m-  1.  ,  (2.25) 
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Since  node  /  must  satisfy  (2.22),  it  follows  that  there  exists  an  m,  1  <  m  <  N,  and  a  time  /“ 
such  that 


As  a»i  easy  consequence  of  constraint  (2.22),  the  path  as  generated  by  our  fcrward-through- 
trnr  a  gument  must  be  finite.  (Loosely  worded,  Constraint  2.22  insures  that  a  node  in  H  can 
nv:  be  visited  more  than  once  for  any  feasible  flow.)  Hence,  ihe  path  must  be  of  the  form 

(k,l)  — *  (.\m)  —  ( m,n )  —  •  •  •  —  (w,x)  —  (xvv)  —  Ov)  • 

The  forward-through-time  argument  clearly  indicates  that  Since  no  inventory  is 

allowed  except  at  the  sink  node  and  (yj)  is  the  last  node  in  the  path  it  is  easy  to  see  that 
(yj)  —  W+1JV+1).  By  (2.23),  it  follows  that  Or  j')—  (*,/).  Hence,  it  is  easy  to  verify  that 

+  ••  ‘  “  1  (2.33) 

+  +  •  •  •  +<*(. .,*))  *  1.  (2.34) 

We  see  that  the  forward-through-time  argument  constructs  from  the  feasible  flow  a  cycle 
in  H  which  we  symbolize  by 

k  —  l  •—  m  —  n  —  — w  —  Ac. 

We  must  now  argue  that  this  cycle  is  in  fact  a  Hamiltonian  cycle.  Intuitively,  Constraint  2.22 
insures  that  a  node  in  H  has  to  be  visited  at  least  once  for  any  feasible  flow. 

To  argue  more  formally,  let  /  denote  a  node  in  H  not  in  the  cycle  constructed  above. 
Since  /  must  satisfy  Constraint  (2.22),  it  follows  that  for  some  h,j  and  t~ 

1.  .  (2.35) 

Instead  of  using  a  forwards- through-time  argument  starting  at  r“,  we  now  use  a  backwards- 
through-time  argi  -vents.  Constraint  (2.35)  implies  that 

—  d{,j))  -*.l. 

By  (2.21),  y(*,o-0,  hence  by  (2.36)  —  1.  This  means  that  for  some  g 


(2.36) 


V(s  .*).(*./)  w  —d(,j))  —  1 . 


(2.37) 


It  is  easy  to  see  how  to  propagate  the  backwards-through-time  argument. 

The  path  in  H  generated  from  the  backwards-through-time  argument,  by  (2.22),  must 
terminate.  Furthermore,  no  node  in  this  path  can  be  a  node  in  the  cyde  we  generated  from  the 
forward-through-time  argument  (again,  by  (2.22)).  If  ( a,b )  denotes  the  first  arc  in  the  path 
generated  by  the  backwards-through-time  argument,  it  is  easy  to  check  (and  intuitively  dear) 
that 

><..»)(£»- 1.  (2.38) 

And,  of  course,  if  (a,b)  *  ( k,l )  then  (2.21)  is  violated;  if  (e.b)-  (*,/)  then  (2.22)  is  violated 
for  node  l.  Hence,  our  cycle  is  Hamiltonian. 

Finally,  we  wish  to  argue  that  r*  is  in  fact  equal  to  d(tJ).  If  so.  then  we  wiU  have  shown 
that  the  intermediate  product  transfer  into  the  sink  node  is  the  indicator  of  the  length  of  the  Hamil¬ 
tonian  cycle. 

By  (2.26),  we  must  have  that 

l^u,/>(y<*,/)*^u,/))J(f  )“1.  (2.39) 

Since  Wik  0  —  t  +^(/,„)  +  <*<„,„)  +  •  •  •  d{9ik)  (2.33),  by  definition  of  (see  2.24)  it  is  dear 
that  either  t  •md{sj)  or  t  —  r  +  </<*>n  +  Ju%m)+  • +  d{w^.  Since  d(tJ)  >  0,  V(/  j)  €  A  we 
have  that 

The  conclusion  is  that  a  feasible  flow  corresponds  to  a  Hamiltonian  cyde.  It  is  immediate 
that  a  Hamiltonian  cyde  corresponds  to  a  feasible  flow.  One  smply  chooses  one  arc  in  the 
cycle  as  a  starting  point,  say  (*,/),  and  sets  y<*,,)(0)  -  1.  If  (/,m)  is  the  next  arc  in  the  cycle, 
then  one  sets 


Vlk.l)Al.m)(di,m))  —  1 
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and  continues  in  the  obvious  way.  That  the  inventory  balance  constraints  (2.1M2.11)  are  met 
is  immediate.  The  proof  of  Proposition  (2.2)  is  now  complete. 

Let  LN  refer  to  the  network  dynamic  production  correspondence  fot  the  production  sys¬ 
tem  described  in  Proposition  (2.2).  Then  as  an  immediate' corollary  we  have  the  following  fact. 

Corollary  (2.3) 

For  all  r  6  /?+,  if  u  —  1|T)  then  LN(u)  ^  0  if  and  only  if  there  exists  a  Hamiltonian  cycle 
in  H  of  length  less  than  or  equal  to  r. 

As  a  consequence  of  our  above  developments,  the  following  problem  is  equivalent  to  the 
Traveling  Salesman  Problem: 

P:  min  j  J  udp,  I  3x  with  x  €  LN(u)j . 

Hence,  this  last  problem  is  NP-complete.  Since  finding  the  longest  Hamiltonian  cycle  is  also 
NP-complete,  then  the  general  problem  of  maximizing  scalar  output  for  a  production  system  is 
NP-complete. 

2J.3.  Material  Requirements  Planning 

Underlying  Assumptions 

Material  Requirements  Planning  (MRP)  is  a  production  planning  tool  for  the  discrete  parts 
manufacturing  environment.2  To  use  MRP,  the  production  planner  assumes  that: 

(1)  All  Sows,  except  the  intermediate  product  transfer  flows,  are  step  functions 
associated  with  an  equal  length  period  time  gnd.  The  intermediate  product 

transfer  flows  are  assumed  to  be  event-based  flows.  The  number  of  periods  p 
1  For  a  comprchwuiv  pnmnuUon  m  OrUcky  (1970,1975)  or  Ptont,  WipM  (1971). 
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is  assumed  finite. 

(2)  Bill-of-Materiah  Coefficients  fl»,  /«1,2 . m,  k  — 1,2. . . . ,  m  exist  such 

that  to  produce  1  unit  of  product  k  requires  at*  units  of  product  '. 

(3)  Constants  Lk ,  k  -  1,2 . m  exist  such  that  it  takes  no  more  than  Lk  periods 

to  produce  an  arbitrary  quantity  of  product  k  regardless  of  the  magnitude  of 
that  quantity. 

(4)  If  a  production  planner  uses  MRP,  then  the  planner  is  tacitly  assuming,  by  (2) 
and  (3),  that  the  production  network  associated  with  the  discrete  parts 
manufacturing  environment  is  one  for  which  activity  Ak  represents  the  produc¬ 
tion  of  product  k. 

Constraints  «  the  later  mediate  Predict  Transfers 

If  activity  Ak  is  to  produce  J  V*dn  units  of  product  k  in  period  /„,  then  it  is  required 

O 

that  all  intermediate  products  necessary  to  produce  J  Vkdfi  units  of  product  k  for  A,  must  be 

f 'J 

available  Lk  periods  earlier.  That  is,  for  1  <  r  <  p,  1  <y  <  m,  1  <  /  <  m, 

±a„fv*dti<  t  J  |Z^W.  (2  40) 

— »  r?j  — t.  ('-4(i|  >  j 

(2.40)  could  be  suitably  modified  to  incorporate  transfer  lap  if  necessary. 

Cxplesiea  ef  Require  meats 

Let  v  denote  a  fine!  output  vector.  If  ooe  insists  on  equality  in  Constraint  (2.40),  then  it 

is  easy  to  check  that  Assumption  4  implies  that  all  K/s  are  determined.  This  determination  is 

usually  referred  to  as  the  material  requirement's  explosion  or  the  explosion  of  requirements. 

The  Fundamental  PreMcm  with  MRP 
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The  use  of  MRP  only  determines  a  schedule  for  the  P^'s.  It  says  nothing  about  how  the 
individual  activities  which  comprise  the  production  system  will  be  able  to  produce  the  require¬ 
ments.  Stated  differently,  by  Assumption  3,  MRP  ignores  the  activity  production  functions 
since  it  assumes  that  each  activity  may  be  able  to  produce  any  quantity  required.  By  choosing 
to  ignore  the  activity  production  functions,  the  production  lags,  the  Lk' s,  are  introduced  and 
inflated  to  insure  feasibility  of  the  schedule  of  the  Kv’s.  It  is  well-known  that  this  is  the  funda¬ 
mental  problem  with  MRP.3 

2.3.4.  Single-Project  Production  Systems 

A  single-project  production  system  is  one  for  which  activities  are  assumed  to  perform  a  one¬ 
time  job,  uninterrupted  from  start  to  finish,  and  for  a  known  fixed  duration.  A  project  is  said  to 
be  completed  when  all  of  the  activities  are  finished.  Following  Leachman  [1983],  it  is  assumed 
that  each  exogenous  resource  is  non-storable,  that  the  total  amount  of  each  exogenous  resource 
used  to  complete  an  activity  is  known,  and  that  the  resources  are  applied  at  constant  rates. 

Associated  with  a  single-project  production  system  is  an  aanity-on-node  precedence  network 
denoted  by  H.  H  is  an  acyclic,  directed  graph  on  N  nodes.  Node  A,  in  H  corresponds  to 

activity  i,  i  — 1,2 . N.  The  arcs  indicate  stria  precedence  relationships  between  the  activities. 

That  is,  if  there  is  an  arc  from  A,  then  activity  J  cannot  Mart  until  activity  i  has  finished.  Let 
W  denote  the  set  of  whole  numbers.  Let  d,  denote  the  duration  of  activity  i,  d,$.  W  for  each 
i. 

Definition  (2.41) 

A  feasible  schedule  of  sun-times  for  the  activities  is  a  vector  S»(S\ . S»)  €  such 

that  if  arc  (ij)  exists  in  //  then  S,  +  d,  >  S,.  We  note  that  by  taking  S  €  Wl,  we  are  assuming 
that  a  start-time  for  an  activity  always  begins  at  the  beginning  of  a  period. 

1  For  •  deutiod  dwcnpUow  of  ihr  yoOUnvi  with  MRP,  too  Kaapor  [1H31. 
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In  this  section,  a  production  system  is  constructed  such  that,  for  an  appropriate  choice  for 
the.  finaJ  output  vector  u,  there  is  a  I- 1  correspondence  between  the  set  of  feasible  flows  to 
support  output  level  u  and  the  set  of  feasible  schedules  of  start-times.  Furthermore,  certain 
feasible  flows  may  be  identified  with  the  well-known  early-start  and  late-start  schedules  deter¬ 
mined  by  the  ordinary  critical  path  method.4 

We  first  proceed  to  construct  the  production  network  G  associated  with  the  production 
system.  To  the  activity  on-node  network  H  add  sink  node  AN+\  and  add  arcs  (/,\+ 1)  from 
each  node  A,  to  node  AN+ j.  This  network  will  be  taken  to  be  the  production  network  G . 

To  define  the  production  system  associated  with  G ,  we  need  to  specify  the  limitations  on 
the  primitive  elements.  We  list  these  below: 

The  time  grid T  is  the  set  {0,1,2,  •  •  • ). 

(2)  All  flows  are  step-functions  except  the  final  output  vector  which  is  event- 
based. 

(3)  No  disposal  of  any  kind,  is  allowed. 

N 

(4)  Attention  is  restricted  to  the  finite  horizon  [0,A]  where  *  —  £</,. 

.-i 

The  additional  constraints  imposed  *oo  the  flows-other  than  those  constraints  already 
imposed  by  the  general  model-are  listed  below. 


Ob  the  Applications  rf  System  Exagnuas  Inputs 

_  As  notation,  let  A*,,  *—1,2, . . .  ,n  denote  the  total  amount  of  resource  *  that  A, 
requires.  As  notation,  let,  for  /- 1,2. . . .  ,N, 


and  define  for  i-N+1 

4  Sac  Modtr  md  Phillip*  (1970),  for  cumpto. 


,  €  ir  (m)  I  3$  €  such  that 


(2.42) 
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2n+ i  “  {*v+i  ^  L“  (aO  I  BSWi  €  /?+  such  that  z*+|™  1  csw^,.*> 


(2.43) 


Then  it  is  required  that  there  exist  a  z,  €  Z,  such  that 


v*-A,2,,  1  <*<«•  (2.44) 

Note  that  the  description  of  the  application  vector  fits  the  form  as  required  by  DLAAM 
presented  in  Section  2.3.1.  The  variable  z,  is  referred  to  as  an  operating  intensity.  The  interpre- 

t 

tation  of  z,  is  that  f  z.du.  expresses  the  fraction  of  the  total  amount  of  each  resource  required 
o 

to  complete  A,  up  to  time  r  or  expresses  the  fraction  of  the  job  completed  up  to  time  r. 

On  the  Activity  Production  Functions 

To  define  F,(y,,W,),  one  needs  to  know  what  kinds  of  'productsf  A,  is  ’producing.*  It  is 
assumed  that  each  activity  produces  a  distinct  product  for  each  of  the  activity’s  immediate  suc¬ 
cessor!.  Following  DLAAM,  the  output  rue  is  modeled  as  the  rate  of  the  utilization  of  the  sys¬ 
tem  exogenous  inputs. .  Thu  is. 


2,  if  arc  (/,/)  is  in  G  and  y*-  bk,z,  for  each  k 
WM-  0  otherwise. 


(2.45) 


In  effect,  (2.45)  says  that,  up  to  time  r,  of  an  activity  has  been  completed  if,  up  to  time  r, 
x%  of  the  resources  reqttired  to  complete  the  activity  have  been  utilized. 


On  the  Intermediate  Preduct  Transfers 

Each  activity  A„  i  U  N+ 1,  obtains  a  distinct  intermediate  product  input  from  each  of  its 
immediate  successors. 


Doftnitiea  (2.44) 
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An  intermediate  product  transfer  functional  is  a  map  fJt  :Zj  —  Z,  defined  by 

df ' 1  •*/*/*.  ’  '  *  N+ 1 

It  is  then  required  that  if 

>!/*- Vy 
y*  -  bk,2, , 

then  for  l,2,...,n,  i«-l,2, . . .  .AM-1, 

^ -/>(*,)  (250) 

»V-2(.  (2.51) 

This  completes  the  description  of  a  single-project  production  system.  A  routine  check 
verifies  that  the  axioms  are  satisfied.  As  notation,  let  F  denote  the  subset  of  indices 

(1,2 . m)  associated  with  the  intermediate  product  transfen  into  the  sink  node.  The  fol- 

lowinj  proposition  and  its  corollaries  show  why  we  have  chosen  to  define  the  primitive  ele¬ 
ments  in  the  manner  given  above. 

PrapositlM  (2.4) 

If  the  final  output  vector  w  is  defined  by  u'  —  0  if  i  t  F  and  u'  —  1(A)  if  i  €  F,  then  there  is 

a  1-1  correspondence  between  the  set  of  feasible  flows  to  support  u  and  the  set  of  feasible 

schedules  of  start-times. 

PreW  af  Propositi**  (2.4) 

Assume  feasible  flows  to  support  output  u  exist  and  select  one.  Each  y,,  by  constraint 

(2.44),  has  an  operating  intensity  associated  with  it.  By  definition  (2.42),  each  i,  has  an 


1  <  k  <  a 

(2.48) 

1  <*  <n 

(2.49) 

S,  €  R+  associated  with  it.  We  claim  that  the  vector  S  -  (Sj . S>  )  so  generated  is  a  feasible 

schedule  of  start-times. 

Suppose  Aj  is  an  immediate  predecessor  to  A,.  By  the  definition  of  a  feasible  flow  the  fol¬ 
lowing  inventory  balance  constraint  must  hold: 

T 

I  \Vj~WJ\diL  >  0,  V/ , Vr €  10 ,h  1 .  (2.52) 

Jo 

By  the  constraints  imposed  on  VJt  and  W{  (see  (2.50),  (2.51))  it  follows  that  upon  substitution 


t 

Ji {  4"  “  df  >  0  •  Vt  €  [0,hl . 


(2.53) 


Gearly  this  implies  that 


S,>Sj  +  dj.  (2.54) 

Hence,  5  is  a  feasible  schedule. 

Conversely,  if  5  is  a  feasible  schedule  of  start-times  then  define  the  induced  operating 
intensity  for  A,,  denoted  by  r, ,  to  be 


tf-1  ti.W  (2-55) 

A  routine  check  verifies  that  the  primitive  elements  defined  from  the  operating  intensities  given 
in  (2.55)  satisfy  the  inventory  balance  constraints  of  the  general  model  (2.1-2.11).  Thus,  a 
feasible  schedule  5  induces  a  feasible  flow  in  a  natural  way. 

The  association  between  feasible  flows  and  feasible  schedules  given  above  is  easily  seen  to 
be  1-1.  The  proof  is  now  complete.* 


Definition  (2.54) 
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An  Early- start  schedule,  denoted  by  £,  is  a  feasible  schedule  such  that  if  5  is  any  other 
feasible  schedule  then  5,  >  E,  for  /-  1,2 . N. 

Definition  (2.57) 

A  Late-stan  schedule,  denoted  by  L ,  is  a  feasible  schedule  such  that  if  5  is  any  other 
feasible  schedule  then  S,  <  L,  for  /  — 1,2 . N. 

The  following  two  corollaries  are  immediate: 

Corollary  (2.5) 

A  feasible  flow  minimizes  the  sum,  over  all  activities,  of  the  cumulative  intermediate 
product  inventories  at  the  sink  node  if  and  only  if  the  schedule  of  start-times  associated  with 
feasible  flow  is  the  late-start  schedule. 

Corollary  (2.6) 

A  feasible  flow  maximizes  the  sum,  over  all  activities,  of  the  cumulative  intermediate  pro¬ 
duct  inventories  at  the  sink  node  if  and  only  if  the  schedule  of  start-times  associated  with  the 
feasible  flow  is  the  early-start  schedule.  For  ease  of  prwentation,we  excluded  it. 


2.4.  Comments  on  the  General  Model 


Comments  about  the  general  model  and  extensions  to  the  model  are  made  below. 

On  i  Stochastic  Framework 

By  our  assumption  and  definition  of  the  activity  production  functions,  our  framework  is 
clearly  deterministic.  We  are  assuming  that  if  we  knew  the  applications  of  input  that  we  would 
know  the  realized  output  obtained  from  production.  Environments  such  as  agriculture  where 
weather  plays  an  important  role  in  determining  yield  do  not  fit  this  assumption. 

One  can  make,  the  production  functions  random  functionals  to  account  for  such  environ¬ 
ments.  If  we  let  ft  denote  some  probability  space1  then  we  can  define  F,  as  a  rasp  from 

L,t  xLWiX  Cl -~R+ 

where  one  interprets  F/0',,lF,,w)  to  he  the  actual  output  if  was  the  vector  of  inputs 

applied  and  w  was  the  state  of  nature  observed.  We  are  tacitly  assuming  here  that  ft  is 
independent  of  L,  x  L*  ,  i.e.,  ft  does  not  vary  with  the  choice  of  (y,,lF().  Axioms  on  F,  would 
have  to  be  suitably  modified  to  account  for  the  stochastic  formulation  of  F,.2  To  model  produc¬ 
tion  systems  where  it  is  more  appropriate  to  view  ft  as  dependent  on  the  choice  of  0^, WJ  is 
considerably  more  difficult  but  worth  investigating. 

Lags 

As  we  have  already  mentioned,  no  explicit  mention  of  time  lags  was  given  in  the  general 
model.  Such  lags  could  be  easily  incorporated  but  would  require  the  introduction  ,  of  new  vari¬ 
ables  in  a  manner  similar  to  that  treated  in  Section  11.4.  There,  variables  K*  were  introduced; 


1  ft  -  tY,Bf)  where  Y  is  a  ret.  B  is  •  w— algebra  o t  subsets  of  Y  end  P  is  •  measure  such  that  P  is 
w-flnite  with  HY) »  1. 

2  Mak  11981]  develops  a  stochastic  theory  of  Dynamic  Production  Correspondences. 


vs 
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Dependence  Between  Activity  Primitive  Elements 

As  illustrated  in  our  description  of  the  Traveling  Salesman  Problem,  specific  production 
systems  might  have  additional  constraints  linking  primitive  elements  belonging  to  different 
activities.  By  changing  our  definition  of  the  domains,  we  could  allow  for  this  flexibility.  For 
ease  of  presentation,  we  excluded  it. 


3.  TOPICS  IN  PRODUCTION  THEORY 

This  chapter  addresses  some  of  the  issues  found  in  Production  Theory.  In  Section  3.1  we 
analyse  the  technically  efficient  subset.  In  Section  3.2  we  carefully  analyse  Shephard's  duality  in 
both  the  steady-state  and  dynamic  cases.  Section  3.3  addresses  2  versions  of  the  Law  of  Dimin¬ 
ishing  Returns  as  first  formulated  by  Shephard.  Finally,  in  the  Appendix,  Section  3.4,  some  of 
the  technical  propositions  are  proved. 

3.1.  On  The  Technically  Efficient  Subset 

3.1.1.  On  the  Nonemptiness  of  the  Technically  Efficient  Subset 
The  technically  efficient  subset,  denoted  by  EN(«),  is  defined  to  be 

EN(u)  —  |x  €  LN(u)  |  ify  then y  i  LN(u)  .!  (3.1) 

The  technically  efficient  points  are  those  for  which  it  is  not  possible  to  lower  the  input  rates  and 
still  achieve  the  same  output  rates.  Thus  it  is  important  to  know  that  our  general  framework  is 
consistent  with  a  desirable  property  that  the  efficient  subset  is  nonempty  whenever  the  level  set 
itself  is  nonempty.  This  result  is  also  the  first  stepping  stone  towards  proving  Shephard’s  Dual¬ 
ity  Theorem  in  the  finite  horizon  dynamic  case  which  is  presented  in  the  next  section. 

’  We  remark  that  we  prove  this  theorem  under  the  assumption  of  a  finite  horizon.2  That  is, 
for  some  h  >  0  the  underlying  space  of  flows  of  goods  and  services  as  presented  in  Axiom  1, 
Chapter  3  will  be  changed  to  LZ  ([0,Al,.B,X  + v)  where  B  is  the  restriction  of  the  Lebesgue 
cr~ field  to  [0,A],  X  is  Lebesgue  measure  restricted  to  (10,A),B)  arid  *  the  counting  measure 
on  U*  I  tk  <  A} .  As  notation,  let  LZ  0*)  denote  LZ  (10,A),J9,X  +  *). 

1  U>  <  x,  then  V«, 

1  We  do  not  have  a  proof  in  the  infinite  horizon  case.  Addi  tonal  hypotheses  on  the  primitive  elements  we 
believe  would  have  to  be  taken. 


Theorem  3.1.1 


Fcr  the  finite  horizon  case,  if  LN(u)  ^  Z  then  EN(u)  ^  Z. 


Proof  of  Theorem  (3.1.1) 

Pick  an  arbitrary  z  in  the  nonempty  LN(u).  Consider  the  following  optimization  prob¬ 
lem. 

Piz):  infJlX^J^M 

j 

subject  to: 

(1)  x€LM(u) 

(2)  x  <  z 

It  is  clear  that  if  we  can  show  that  there  exists  an  x  €  LN(«),  x  <  z  achieving  the  infimum  in 
P(z)  then  x  €  EN(w)  and  EN(u)  is  nonempty.  What  we  have  to  show  is  that  this  problem  has 
a  solution. 

To  show  this,  endow  each  L,  (the  subset  of  L?  (41)  satisfying  Axioms  2-4  for  the  /rt  flow 
type)  with  the  relative  weak-star  topology  (viewing  it  as  a  subset  of  L"0*))  and 
Lx  x  •  •  •  x  L;  with  the  product  topology.  Letting  Lx  stand  for  Lx  x  •  •  •  x  I£  then  the 

objective  function  4  viewed  as  a  map  from  L„  —  R+  is  contiguous.3  As  an  application  of  the 


theorems  of  Banach-Alaoglu  and  Tychonoff,  the  set  of  x  in  Lx  satisfying  (2)  is  contained  in  a 
compact  subset  of  (JL**0i))".4  In  the  Appendix  we  verify  that  each  Lg.  is  closed  in  LmUi)  (Pro¬ 
position  3.4.2);  hence,  the  set  of  x  in  Lx  satisfying  (2)  is  compact  in  Lx.  Since  a  continuous 
real-valued  function  on  compact  set  achieves  its  infimum  then  in  otJ-.,  u,  complete  the  proof 
that  EN(u)  is  nonempty  it  suffices  to  show  that  LN(ii)  is  a  weak-star  dosed  subset  of 
(I-Qi))*.3 


1  The  operation  of  addition  is  continuous.  The  integrable  function  here  isp> 

4  The  theorem  of  Banach- AlaocJ  j  states  that  the  dosed  twit  ball  in  the  dual  space  of  an  arbitrary  topologi- 
cal  vector  space  is  weak -star  compact  (Rudin  (19731,  p.  <>6).  Since  4”  is  the  dud  to  l1  we  have  that 
U I  llrrll.  <  1}  is  weak -star  compact.  The  theorem  of  Tycho  no  fT  is  that  an  arbitrary  product  of  compact  spaces 
is  compact  in  the  product  topolofy  (Dujrndgi  (1966),  p.  224). 

1  See  Dugundgi.  p.  227. 
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To  show  that  LN(u)  is  closed  let  {xj  be  a  net  in  LN(u)  converging  to  x.6  To  say  that 
xa  €  LN(u)  for  all  a  means  that  for  all  a  there  exists  a  feasible  flow  for  xa  to  support  output 
level  u.  In  order  to  show  that  x  is  in  LN(u)  (which  is  our  goal)  we  need  to  exhibit  a  feasible 
flow  for  x  to  support  output  level  u.  To  arrive  at  this  end  we  will  extract  "limit"  functions  for 
each  flow  type  from  the  feasible  flows  for  each  x„. 

By  the  very  definition  of  a  feasible  flow  the  net  {xj  induces  a  net  of  functions  for  each 
flow  type.  What  we  will  in  effect  show  is  that  there  exists  a  subnet  of  the  original  net  (xj  such 
that  the  associated  sub-nets  of  functions  for  each  flow  type  have  the  property  that  they  are  uni¬ 
formly  bounded  in  the  norm.  Hence,  by  compactness,  we  may  extract  these  sub-nets  to 

obtain  limit  functions.7  By  duality  between  Ll(p.)  and  the  inventory  balance  constraints 

will  hold  for  these  functions.8  The  proof  would  then  be  complete. 

The  inventory  balance  constraints  which  must  be  satisfied  by  the  various  pairs  of  flow 
types  which  when  taken  together  make  a  flow  feasible  have  the  following  structural  form: 

r 

0^  S  +  /{x-r~y}d/*  <  C,  Vr€l0,/t].  (3.2) 

0 

Here  5  represents  the  initial  stock,  if  any,  C  the  capacity  (perhaps  infinite),  x  the  flow  of  input 
into  th6  "system",  T  the  flow  of  disposal  out  of  the  system,  and  y  the  flow  of  input  out  of  the 
system.  (See  Section  3.2.3.) 

To  arrive  at  our  main  result  we  first  prove  a  simple  proposition  which  paves  the  way  for  a 
"verbal"  proof  of  the  rest  of  the  theorem. 


*  We  use  ne*'  instead  of  sequences  since  1“(m)  as  topologixed  does  not  have  a  countable  neighborhood 
base.  For  definition  of  a  net  see  Halmos  [1974),  p.  65.  For  those  unfamiliar  with  nets  think  of  sequences  in¬ 
stead. 

7  We  know  that  the  closed  unit  bell  is  compact.  Operations  of  scaling  are  homeomorphisms  in  any  topo¬ 
logical  vector  space. 

*  We  are  using  the  fact  that  L'Ui)  is  the  pre-dual  of  LmUi)  hence  with  l#*]  the  convergence  is  main¬ 
tained. 
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Proposition  (3.1.2) 

Let  {*„).  IvJ.  and  { 7'„)  be  nets  in  the  (appropriate)  spaces  Lx,  Ly,  and  LT  which  for  each 
a  satisfy  (3-2).  Further  assume  that  {xj  x.  Then  subnets  of  the  original  nets  exist  which 
are  uniformly  bounded  in  the  Lm Ui)  norm. 

Proof  of  Proposition  (3.1.2) 

Let  e  >  0.  For  some  0  we  have  that  eventually  Va  >  0,  f  x.rf/x  <  J xdp  +«  by  conver¬ 
gence.  Thus,  eventually,  J*  yadn  <  J  xdn  +  c  and  J  Ttdn  <  J*  x + « .  By  Axiom  2  it  now  fol¬ 
lows  easily  that  eventually  the  norms  of  ya  and  Tm  (and  of  course  xj  must  be  uniformly 
bounded.* 

Now  on  to  the  conclusion  of  the  proof  of  the  theorem.  First  treat  the  storage  case.  Since 
x{  —  xJ  for  each  j  ;t  follows  easily  by  Proposition  (3.1.2)  that  we  may  find  an  appropriate  sub¬ 
net  so  that  the  x&  and  yf  sub-nets  are  uniformly  bounded.  By  Axiom  8,  each  of  the  induced 
Wf)  sub-nets  are  uniformly  bounded  too.  Since  there  is  an  inventory  balance  constraint 

.  J-N+ 1 

of  type  (3.2)  linking  the  with  the  J  K*  's  it  follows  by  Proposition  (3.1.2) 

>-i 

that  a  further  refinement  of  the  original  net  gives  us  a  uniform  bound  on  the  nets.9  The 
same  argument  may  be  used  to  deduce  the  same  result  for  the  W,  and  T,  sub-nets.  Now  since 
all  nets  of  functions  of  each  flow  type  are  contained  in  compact  subsets  we  may  extract  conver¬ 
gent  sub-nets. 

If  y,  and  W,  denote,  for  each  i,  the  limit  vectors  of  the  sub-nets  (yj,  { Wf),  then  by  Pro¬ 
position  (3.4.6)  in  the  Appendix  to  this  Chapter  we  see  that 

<y, 

*  Actually  the  tame  remit  could  be  obtained  If  we  invoked  Axiom  12.  Thix  would  immediately  imply  that 
the  nets  [Fflyf,  ■*'•)),  I  VJ\,  and  1  •*',•)  are  uniformly  bounded  (since  we  are  restricted  to  a  finite  horizon). 

See  the  comments  at  the  end  of  Axiom  12. 


(3.3) 
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In  view  of  (3.3)  it  now  follows  easily  that  the  limit  functions  together  satisfy  all  of  the  inven¬ 
tory  balance  constraints  needed  to  define  a  feasible  flow  for  x  to  support  output  level  u. 
Hence,  x  €  LN(u)  which  implies  that  LN(u)  is  closed  as  desired  in  the  storage  case. 

If  inventories  (in  some  cases  )  were  not  allowed,  then  the  argument  is  the  same  except 
that  it  is  much  easier  to  obtain  uniform  bounds  for  the  sub-nets.  One  needs  to  insure  however 
that  if  {xj,  {yj  are  nets  such  that  x„  < y„  with  xtt  —  x,  ya  —  y  for  all  a  then  x  is  indeed  less 
than  or  equal  to  y.  This  simple  verification  may  be  found  in  the  Appendix,  Proposition  (3.4.1). 
Our  proof  to  Theorem  3.1.1  is  now  complete.* 


An  immediate  corollary,  which  will  be  needed  later,  is  given  below. 


Corollary  to  Theorem  (3.1.1) 

LN(u)  C  EN(u)  +  (i.*  (#*))"  (3.4) 

We  remark  that  we  used  (L“  0*))*  instead  of  Lx  since  the  difference  of  two  functions  in  Lx 
need  not  be  in  Lx. 


3.1.2.  Extensions  to  Theorem  (3.1.1) 


If  x‘  were 


a  solution  to  problem  Pit),  then  x’  would  be  a  vector  of  system  exogenous 


inputs  which  minimizes  the  total  amount  of  resource  over  all  x  in  LN(u)  p)  (x  I  x  <  z).  It  is 
certainly  possible!  that  x*  might  not  minimize  the  total  amount  of  resource  over  all  x  in  LN(u). 
If  <*:  /,,  —  R+  denotes  the  objective  function,  then  it  would  be  desirable  to  show  that  the  prob¬ 
lem  defined  by 

P:  it\f+(x) 

subject  to:  x€LN(u) 

has  a  solution. 


Actually,  this  is  an  easy  consequence  of  Axiom  2  as  we  now  show.  Suppose  one  couid 
find  an  x  which  was  a  solution  to  P.  Let  {^},  { denote  the  set  of  gk's  and  hk's  (as  defined 
in  Axiom  2)  associated  with  the  flow  type  [xJ).  Since  x  is  assumed  to  be  a  solution  to  P  and  in 
view  of  Axiom  2  it  is  easily  seen  that 

Ibr'II.  <  tnxx  \gU*(2))M(*(2))) 

I*  I  ik  <  ri 

3/(*(z)). 

Hence,  ILxIL-maxlIxML  <  /(d(*))-  If  we  were  to  define  the  following  problem 

J 

P’\  inf$ Or) 
subject  to: 

(1)  x  €  LN(u) 

(2)  llxll.  <  /(4(a)), 

then  by  the  above  it  may  be  seen  that  any  solution  to  P'  is  a  solution  to  P.  But  P‘  has  at  least 
one  solution;  hence,  so  does  P.w 

We  may  say  more  about  problem  P  if  we  note  that  the  objective  function  ^  in  P  has  the 
following  structural  form: 

St  (p,‘XJ)dp  (3.5) 

/-i 

where  the  p J's  were  all  taken  to  be  the  function  For  any  p  —  (pl . p")€  (I'O*))" 

the  functional  expression  given  in  (3.S)  viewed  as  a  map  from  LM  —  R  is  continuous  (by  dual* 
ity).  Hence,  a  solution  to  P  for  general  p  €  exists.  If  p  €  (LxUt))*  was  strictly  posi¬ 

tive,  i.e.,  p'  >  0  for  each  /,  then  any  point  achieving  the  infimum  would  clearly  be  efficient ; 

10  Um  the  proof  of  Tbaoram  (4.1.1).  W«  only  uaad  Um  fact  tht!  moot  x  <  z  thru  «a  could  raatrict 
LN(«)  to  •  compact  wbcct.  W«  git  Um  mom  raatih  by  dtractly  raatridias  the  tat  LN(u)  to  •  dotad  ball. 
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Let  (L'(m))"+  denote  the  subset  of  (I'(,t))*  of  strictly  positive  vectors  and  let 
Lu-Lu |X  •••  *i>ir  Let  us  define  a  correspondence  from 
< p,u )  €  (!’(/*))++  x  Ly  —  <l>(p,u)  €  Ix  where  one  defines  $(p,u)  as  the  set  of  all  minimizers 
to  Problem  P  with  the  objective  function  as  in  (3-5).  The  p*'*  may  be  thought  of  as  weighting 
factors  for  the  resources  which  take  into  consideration  their  relative  value  to  the  production 
planner  and  the  fact  that  future  costs  may  need  to  be  discounted.11  By  our  above  comments  we 
know  that  each  point  in  (J  +(j>,u)  is  efficient  for  LN(u)  and  has  the  property  that  it 

can  be  obtained  through  minimization  of  an  appropriate  optimization  problem.  Alternatively, 
and  on  more  geometric  grounds,  one  can  say  that  each  point  in  (J  +(p,u)  has  the  pro- 

perty  that  there  exists  a  hyperplane  which  supports  LN(u)  at  this  point.12  An  immediate 
theoretical  question  is  whether  all  efficient  points  can  be  obtained  in  this  way,  i.e.,  is 

\J  *(p»-EN(u)? 

If  this  were  true,  then  we  would  have  a  convenient  characterization  of  the  efficient  subset. 

A  simpler  theoretical  question  to  address  is  whether  or  not 

U  +(p,u)  O  EN(u), 

i.e.,  does  each  efficient  point  have  a  hyperplane  which  supports  LN(ti)  at  the  point  (with  p  not 
necessarily  strictly  positive)?  In  the  next  section,  we  construct  an  example  which  shows  that  in 
the infinite  horizon  case  the  answer  to  this  question  is  no-  even  if  we  assume  that  the  level  set 
is  convex  and  compact.  The  counter-example,  being  constructed  for  the  infinite  horizon  case, 
shows  that  even  if  additional  hypotheses  were  imposed  on  the  primitive  elements  so  that 

11  In  fact  future  com*  were  not  discounted  in  the  previous  example.  Further,  all  of  the  resources  were 
given  equal  ‘weights*. 

11  If  X  is  a  Topological  Vector  Space,  then  a  hyperplane  associated  with  X  is  a  subset  of  X  generated  by 
the  following  equation:  U  I  /Ot)-c|  where  /  f  X'  (the  set  of  continuous  tiacer  functionals  on  X)  and 


52 


Theorem  (3.1.1)  could  be  extended  to  the  infinite  horizon  case  efficient  points  need  not  be 
obtainable  through  an  optimization  problem  of  type  (3.5).  We  remark  that  the  level  sets  would 
be  convex  if  the  production  functions  were  assumed  to  be  concave.  In  the  Dynamic  Linear 
Activity  Analysis  Model  of  Production  (DLAAM)  discussed  in  Section  (2.3.1),  the  production 
functions  were  concave. 


3.1.3.  Counter-Example:  Efficient  Points  Need  Not  Be  Supportable13 

The  mathematical  statement  to  the  question  raised  in  the  previous  section  may  be  phrased 
as  follows.  Let  5  C  CL“  (p))”  be  a  (weak-star)  compact,  convex,  monotonic  subset.  By  mono¬ 
tonicity  we  mean  that  if  x  €  5  and  y  >  x  then  >  €  S.  Let  x  be  an  efficient  point  of  S,  denoted 
by  x  €  ECS).  Does  there  exist  ap  €  (L'Cp))"  such  that 

17  P>  x^dp  <  1  f  P*  z>dp ,  Vz  €  S? 


As  mentioned  in  the  previous  section,  the  counter-example  is  constructed  in  the  infinite  hor¬ 
izon  case  (that  is,  we  now  work  with  the  measurable  space  CR+JB)  instead  of  (lO.AlJI).  To 
construct  the  counter-example  we  will  first  work  in  the  space  f*  and  then  'translate*  into 
LT  0*) .  (It  suffices  to  construct  an  example  in  the  case  when  n  —  1). 

Define  «„  in  r  by 


I  ify»s 
-1  if  y  ■■  rr— 1 
0  otherwise 


(3.6) 


Let  S— colOf*  denote  the  convex  hull  of  the  «/s.14  Since  the  closure  of  a  convex  set  in  a 
topological  vector  space  is  also  convex  then  5  is  closed  and  convex.1?  We  will  show  below  that 


,J  I  wish  to  gratefully  acknowledge  the  aaetstance  of  Stephen  Boyd.  Department  of  Electrical  Engineering 
and  Computer  Science,  Uruveraty  of  Califomie  et  Berkeley,  in  the  construction  of  the  counter -exam  pie. 

14  If  x  €  S,  then  there  easts  e  At/",  mil,- 1  with  the  number  of  poenive  indices  finite  *uch  that 
j-Tk4«*.  i.e.,  for  each>,  *;-£**«!• 

■’Kudin  119731,  p  II 
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S  is  bounded  in  norm  by  1;  it  would  then  follow  that  5  is  o  mpact.  We  will  show  that  0  €  E  (S) 
but  0  is  nor  supportable  in  the  sense  described  in  the  previous  section.16 

Proposition  (3.1.3) 

0  €£(S). 

Proof  of  Proposition  (3.1.3) 

Suppose  this  were  not  true.  Then  this  would  mean  'hat  we  could  find  a  z  €  S  such  that 

i 

zk  <  0  for  all  k  and  for  some  i,  z,  <  0.  We  can  find  a  uet  {x J  in  S  such  that  x.  —  z.  By  the 

definition  of  convex  hull  we  have  that,  for  some  X.  6  /“  with  /'  norm  equal  to  1  whose 

number  of  positive  indices  are  finite,  xa~£  X  ***  for  each  a.  By  the  definition  of  c„  and  the 

* 

fact  that  x.  —  z  it  follows  that  |(A^-Ai'fl)-zM  —  0  eventually.  (This  immediately  implies 
that  I?'  |  <  1;  thus  llxll*.  <  1  which  shows  now  that  S  is  bounded  in  norm  by  1.) 

The  nets  {xj  and  (x;+ll  are  each  contained  in  the  unit  interval  which  is  compact.  Hence, 
each  net  has  a  convergent  sub-net.  Without  changing  notation  for  the  sub-net  let  c'  and  c'+l 
denote  the  limiting  values.  Dearly,  (c'-c,+l)  —  z'.  Since  z‘ <  0  this  implies  that 
c'+l  >  U  I  >  0.  Find  integer  N  so  that  A'lr'l  >  1. 

Now  by  sequentially  repeating  this  process  for  the  nets  (X^+r),  r—1,2 . N  and  not 

changing  noution  for  sub-nets  we  have  that  there  exist  constants  c'*',  r  —  1,2, . .  satisfy¬ 
ing 

1—1,2,...,*  (3.7) 

-s'*',  r—1,2 _ *  (3.8) 

Since  r  <  0  the  constants  c1*' ,  r- 1,2, . . .  ,*  are  non-decreasihg.  Since  c‘*f  >  |z'|  this  in 

■*  By  Krein-MiUman’t  Theorem  (tee  Rudin,  11973),  p.  70),  a  doted,  convex  compact  tubtet  of  t  topolof- 
teal  vector  space  has  to  be  the  closure  of  the  convex  hull  of  its  extreme  point*.  Here,  the  «,  j  are  the  ex¬ 
treme  points  of  5. 
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turn  implies  that 

£c'+'  >  Nc'+'  >  Mr' |  >  1  (3.9) 

fm\  r~N 

Since  £  x;+f—  £  c'*'  (by  (3.7))  eventually  (in  view  of  (3.9))  £  Xi+'>  1.  But  then  this 

r—  I  r«l  r«l 

means  that  eventually  x„  1 5,  an  obvious  contradiction.  The  result  follows.* 

Proposition  (3.1.4) 

0  is  not  supportable,  i.e.,  there  does  not  exist  a  p  €  /*,  p  #  0  such  that  O-p-O  <  p  x ,  V 
x  €  S. 

Proof  of  Proposition  (3.1.4) 

Suppose  such  a  p  exists.  Then  p-«  i  >  0  which  implies  that  pl  >  0.  Since  />■«„  >  0  for  each 
n  this  in  turn  implies  that  pm+l  > p*  for  each  it.  Thus  the  sequence  of  p“s  is  nondecreasing 
and  nonnegative.  Therefore,  either  p  -  0  or  p  t  /*.  In  either  case  we  would  have  a  contradiction. 
The  result  follows.* 

We  now  turn  to  constructing  our  counter-example  in  L“  (p).  Let 

l+S-ly€r|  3x  €  S  such  that p,-x,  +  l,  Vi).  (3.10) 

Since  translations  in  topological  vector  spaces  are  homeomorphisms  which  preserve  convexity 
and  efficiency  then  1  is  an  efficient  point  of  1  +  J  and  is  not  supportable  in  the  sense  of  Propo¬ 
sition  (3.1.4).  Now  we  may  identify  1  +  5  with  a  compact  subset  of  LZ  (m)  in  the  obvious  way 
by  associating  each  x  €  1 + S  with  x(t)  €  LZ  (p)  defined  by 

*  -  (3.11) 

i-i 

(Here,  the  set  of  time  points  defining  the  possible  times  of  event-based  transfers  is  the  set  of 
nonnegative  integers.)  Let  5*  denote  the  identification  of  1+S  in  LZ  0*).  Then  S*  has  the 


desired  properties  we  were  seeking  in  our  construction  of  a  counter-example. 

This  example  illustrates  the  problem  of  additional  flexibility  available  by  infinite  time  sub¬ 
stitution.  The  step  functions  c,  that  generate  S'  are  those  which  are  1  everywhere  except  in  two 
consecutive  periods  (/-l,/),  (i,/'+l)  where  the  function  is  0  on  0-1,/)  and  2  on  0V+1). 
Interpret  the  functions  in  Ll(n)  as  price  functions  and  the  optimization  problem  given  in  (3.5) 
as  cost  minimization.  For  any  price  function  p  6  L](p.)t  0,  the  cost  of  the  function 

I  -  X  l(,—i.,)  is  llpllj.  Furthermore  for  some  i,  J p*l((_ii()  >  J* p-l(,.,+i)  (if  not,  then  the  proof 

f—1 

of  Proposition  (3.1.4)  shows  that  either  p  —  Ootpi  L'ip.)).  Hence 
Cost  of  X  1<i-i,i)“  X  J" P‘ltk-i.k)dn 

<-t  *-i 

“  X  /  P’l(k-ljt)dn  +  J*  +  J  p\U'H.\)dp 

k  *  1,1+1 

>  X  f  p-U.k-l,k)dfi  +  2 f  p-l{i^t)dp. 

*««.,/♦  t 

—  Cost  of  . 

Thus,  c,  would  be  cheaper  than  1  so  that  1  could  never  be  obtained  from  cost  minimization, 
i.e.,  1  is  never  supportable.  The  reason  such  a  cost  tradeoff  as  illustrated  above  can  always  be 
done  is  that  tne  generators  c,  of  5*' "go  out  towards  infinity." 

Next,  we  investigate  a  property  of  our  counter-example  which  is  of  independent  interest. 
Consider  the  following  theoretical  question:  is  it  always  true  that  one  can  separate  two  compact, 
convex  subsets  of  Rl  which  meet  at  only  one  point?  Answer:  No;  counter-example:  two  line 
segments  meeting  at  a  single  point.  This  example  may  be  constructed  because  we  did  not  insist 
that  the  two  sets  in  question  have  non-empty  interior.  It  turns  out  that  our  example  also  has 
this  property. 


Proposition  (3.1.5) 


5  is  nowhere  dense  in  /**,  i.e.,  the  interior  of  5-0. 


Proof  of  Proposition  (3.1.5) 

Let  x  €  5.  Let  {jc„)  be  a  net  in  5  converging  to  x.  Write  where  for  each  a, 

* 

X°€/\  IIX“llj  —  1,  and  the  number  of  positive  indices  finite.  Since  x.  — x  this  means  that 
Xf-Xf  —  xx.  Extract  a  sub-net  of  the  net  {xj  so  that  the  sub-net  {xj}  converges,  say  to  pj. 
Since  Xf-X?  —  xi  it  follows  that  kf  —  p\-xh  In  fact,  the  following  simple  claim  may  be 
made. 


Claim 

VX. 

(-1 


Proof  of  Claim 

We  have  shown  this  to  be  true  when  k  - 1  and  k  —  1.  Suppose  it  is  true  for  all  X  <  m. 

m—l 

Since  X*-XJ+i  —  xm  anr1  by  induction  X£  — pi— £x,  it  follows  that 

«- 1 

m 

X“+1  —  XJ— x*,—  pj— Vx,.  The  induction  step  is  shown  and  the  result  follows. ■ 


N 

For  all  a  and  for  each  N  we  have  that  1  >  £  Xf  which  unifies  by  our  claim  that  for  all 

*-i 

N  t-l  -  *-l  *— 1 

Af,  1  >  X  <Pi~  X*')  hence  1  >  X  <Pt“  X  *.)•  So  if  we  define  p*-pi-  X  then 

*“l  i«*l  *“  1  i«l 

may  be  easily  ran  that  p  €  77,  llpllj  <  1  and  the  expression  x-£pti4  is  well-defined  and 

a 

correct. 

We  have  thus  shown  by  our  above  developments  that 
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The  converse  conclusion  is  simple  to  show;  hence,  we  may  conclude  that 


S-lxi 


(3.12) 


With  this  we  will  now  turn  to  proving  the  proposition. 

If  the  interior  of  S  ^  0,  then  we  can  find  an  x  €  S  and  a  neighborhood  Nx  about  x  such 
that  jVx  C  S.  Find  a  net  in  S,  such  that  x*  -  x.  Since  the  weak-star  topology  is  weaker 
than  the  norm  topology  then  for  some  c  >  0  we  have  that  eventually,  for  some  o, 

xaSB,(xa)CNxCS. 


G early  this  means  there  is  some  x 


€  S  such  that  z  <  x*  (we  fix  some  a  large  enough)  and 


If  z  €  S,  then  by  (4.12)  we  have  that  for  some  p  €  /+ ,  llplli  ^  1  (and  number  of  positive 

indices  finite)  z-J>*«*,  i.e.,  z,-P,-P,+ t-  Since  x“€S  eventually  x“-0.  So,  eventually, 
* 

p,  <  p,+1  since  z  <  x-.  Since  z  *  0,  llplli  <  1  we  have  that  eventually  p,  <  0.  So  the  conclusion 

is  that  the  number  of  positive  indices  of  p  is  finite  and  non-empty  (since  z  s*  0).  It  is  not  hard 
to  see  that  if  p  has  the  above  properties  then  llplli  - 1  so  that  z  €  S. 

Write  x*  as  £n**  for  s°me  appropriate  y.  Since  x°  >  z  we  have 


y<  -y,+i  ^  k/  **  k (>i  i''  Vi , 


(3.13) 


or  that 


y,-*k,  >y/+i-kj+i,  Vi. 


(3.14) 


Since  the  number  of  positive  indices  for  both  y  and  p  is  finite  eventually  yl+,-ki+,-0  so  that 
y  >  A,  Vi  in  view  of  (3.14).  This  in  turn  implies  that  since  £y,-£k,-l  (due  to  the 

finiteness  of  the  number  of  positive  indices)  that  y,«A,,Vi.  Hence,  which  is  a 
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finiteness  of  the  number  of  positive  indice?)  that  y,  -X, ,  V/.  Hence,  x„— z  which  is  a  contrad¬ 
iction.  The  result  follows.* 

We  make  one  remark.  What  we  hpve  shown  is  that  5  C  E(S),  that  is,  all  points  in  S  are 
efficient.  Also,  by  (M2)  it  is  easy  to  see  that  not  all  points  in  S  are  efficient.  For  example,  take 
y,»VS\  Xj  «■%(%)'  and  letx-£y*«*,>-£x*ejk.  Thenx  is  not  efficient. 

3.1.4.  Compactness  of  the  Gosure  of  the  Efficient  Subset 

To  prove  Shephard’s  duality,  it  is  necessary  that  the  closure  of  the  efficient  subset  be 
compact.  In  the  axiomatic  framework  of  the  Dynamic  Theory  of  Production  Correspondences  as 
developed  by  Shephard  and  Fare  (1980),  this  property  of  compactness  of  the  closure  of  the 
efficient  subset  was  assumed  as  an  axiom.  We  however  prove  this  property  from  our  axiomatic 
framework  as  presented  in  Section  2.2. 

Proposition  (3.1.6) 

EN(u)  is.  compact. 

Proof  of  Proposition  (3.1.6) 

Assume  that  EN  (u )  is  nonempty  otherwise  there  is  nothing  to  show.  Let  x  €  EN(u)  and 
let  (x")  be  a  net  in  EN(u)  converging  to  x.  Associate  to  each  x"  a  feasible  flow  for  x"  to  sup* 
port  output  level  u.  From  the  proof  of  Theorem  (3.1.1)  we  know  that  we  may  extract  a  sub-net 
of  the  original  net  {x“)  so  that  the  sub-nets  associated  with  each  of  the  flow  type  nets  induced 
from  the  feasible  flows  for  x*  to  support  output  level  u  converge  and  are  uniformly  bounded. 
We  denote  the  sub-net  again  by  {x°}. 

Since  x"  €  EN(u)  this  means  that  there  could  not  be  any  disposal  of  system  exogenous 
input  in  the  feasible  flow  for  x*  to  support  output  level  u.  Thus,  for  each  J  we  must  have  that 
for  each  a 

j“  CxO-rfM  - J  {£  (x&  )“)</**  -  /  II  W)dv . 

0  0  /  0  ( 


(3.15) 


Furthermore  in  view  of  Axiom  10  it  follows  easily  that 


3,2.  On  Shephard's  Duality 


3.2.1.  Statement  of  Shephard's  Duality  Theorem  (Steady-State) 

In  the  Shephard  model  (Shephard,  [1970]),  the  production  technology  is  characterized  by 
the  Production  Input  sets  or  level  sets  L(u),  u  €  R+.  Liu)  is  to  be  interpreted  as  the  set  of  all 
input  vectors  x  €RJ.  which  yield  at  least  output  rate  u.  The  level  sets  are  assumed  to  satisfy 
the  Allowing  axioms:1 

P.l  L(0)-/?;,0*L(«)  forw>0. 

P.2  x€L(u)  and  x’>x,  imply  x'€L(w). 

P.3  If  (a)  x>  0,  or  (A)  x  >  0  and  (X-x)  €  L(«)  for  some  X  >  0  and  u  >  0,  the 
ray  {X-x  I  X  >  01  intersects  Liu)  for  all  u  €  (0,<»). 

P.4  >  «i  >  0  implies  L(u2)  C  L(ut). 

P.5  D  L(u)-L(uo)  for  u0>  0. 

0<v<w0 

1  >  n  L(u)  is  empty. 

n€ 

P.7  L(u)  is  closed  for  all  u  €  [Q,“>). 

P.8  Liu)  is  convex  for  all  u  €  [0,«). 

P.9  E(u)  —  {x  €  L(u)  |  if  y  <  x  then  y  f£  L(u)}  is  bounded  for  all  u  €  10,<»)}. 

From  the  description  of  the  production  technology  and  the  axm;**  taken  on  the  level  sets, 
the  production,  distance,  and  factor  minimal  cost  Junctions  can  hr.  derived. 

The  Production  Function 

The  production  function  <J>:  RZ  — •  R+  is  defined  by 

4>(x)  -  max{u  |  x  €  L(u)}.  (3.18) 

The  production  function  4>  measures  the  largest  output  rate  obtainable  with  x.2 
1  Shephard  (1970a],  p.  14. 

J  See  Chapter  2,  Shephard  11970]  for  the  proof  that  ♦  is  well-defined. 
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The  Distance  Function 

The  distance  function  4»'R+  x  Rl  —  R+  is  defined  by 

u  >  0  and  {A  |X\x  €  L(u)}>-0 
u  >  0  and  {X  |X-x  €  L(u)}  ^  0  (3.19) 

u-0 

Essentially  4>(uyx)  measures  how  much  one  has  to  scale  x,  up  or  down,  in  order  to  just  obtain 
output  level  u;  that  is,  to  put  ■■■■/-;■  on  the  boundary  of  L(u)  if  possible.  Observe  that  by  the 

4 '(UrX) 

definition  of  4> (u^c)  one  has  that 

L(tt)  —  {x  |  4i(ujc)  >  1} .  (3.20) 

Hence,  <b  is  useful  in  that  it  may  be  used  to  characterize  the  level  sets.3 

The  Factor  Minimal  Cost  Function 

Interpret  p  €  Rl  as  prices  and  p  x  as  the  cost  of  input  vector  x.  If  one  defines  the  factor 
minimal  cost  function  Q.R+  x  Rl  —  R+  by 

Qiujt)  -  min(p-x  |  x  €  L(u)),  (3.21) 

thv'i  Q(uj)  would  be  the  minimal  cost  of  obtaining  output  level  u.4  The  statement  of 
Shephard’s  Duality  begins  to  unfold  when  it  is  realized  that  Q(uj>)  has  all  of  the  properties 
that  +(uyx)  has.3  In  effect,  Q(uj>)  acts  like  a  distance  function.  If  it  would  be  thought  of  as 
such,  then  the  candidate  "level  sets"  for  which  it  would  be  a  distance  function  would  have  to  be 

L0(u)-{p|  C(«J>)>1).  (3.22) 

Thinking  of  L$  (u)  as  price  sets  of  the  cost  structure  much  in  the  same  way  one  thinks  of  L(u)  as 

*  ft  open**  of  4>  may  be  found  in  Chapter  3,  Shephard  11970).  ' 

*  Aa  an  eaey  consequence  of  the  axioms  ft  l-ft.9.  L(e)-E(u)  +  ft;  «£(«)*  Jt$.  Since  E(w)  is  com* 
pad,  it  foUowt  that  Qiuj)  is  well-deAned,  i.e.,  the  minimum  is  obtained. 

*  See  Proposition  16,  Shephard  11970], 


0 

*(w  x)  -  lmin{X  |x-x  €  L(u)}]~‘ 


the  production  possibility  sets  of  the  production  structure,  Shephard's  Duality  Theorem,  loosely 
worded,  is  that  we  may  derive  one  structure  (cost  or  production)  from  the  other.  To  quote 
Shephard  (p.8,  (19701): 

"...  production  possibility  sets  of  the  production  structure  and  the  price  sets  of  the 
cost  structure  are  shown  to  be  duals,  derivable  from  each  other  by  dual  cost  minimi¬ 
zations  which  determine  the  factor  minimal  and  price  minimal  cost  functions  as  dual 
distance  functions." 

To  mathematically  formulate  the  above  statement,  define  the  price  minimal  cost  function 
as6 


+  ‘(upc) 


0  if  L(u)*-0 

iof{px  |  p  €  L0(u)}  ifn>0,L (u)*0 
oo  if  u-0 


Shephard’s  Duality  Theorem  is  then  the  statement  that 

♦(x)  “  max  {u  |  if>’(vpc)  >  1) . 


(3.23) 


(3.24) 


3.2.2.  Shephard’s  Approach  to  the  Proof  of  Duality 

Let  L*(i/)  —  {x  I  +’(urx)  >  1}.  Shephard  first  shows  that  ift'  had  the  properties  of  a  dis¬ 
tance  function.7  Next,  he  shows  that  the  sets  L*(«)  defined  a  production  technology,  i.e., 
L‘(u)  satisfied  axioms  P.  l-P.  9.®  Then,  he  shows  that  L*(ii)  —  L(a)  from  which  the  duality 
statement  may  be  easily  derived  and  from  which  d*  may  be  seen  to  be  equal  to  <fi. 

The  next  two  sections  provide  two  alternate  ways  to  proving  Shephard’s  theorem.  Version 
1  we  believe  is  far  more  direct  and  simpler  to  understand  than  the  original  approach.  Version  2 
motivates  the  definitions  of  d  and  d*  and  in  some  sense  derives  the  identity.  It  shows  that  to  a 

*  Shephard  defined  **  differently.  There  wu  no  mention  of  the  (pedal  cue  when  L(u)  *  •  end  u  >  0. 

7  Proportion  39,  Shc^rd  (1970). 

1  Proportion  40,  Shephard  11970]. 
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large  degree  Shephard’s  theorem  is  a  statement  about  supporting  hyperplanes  of  boundary 
points  of  a  particular  class  of  convex  subsets  of  Rl . 

3.2.3.  Proof  of  Shephard’s  Theorem:  Version  1 

The  proof  of  Shephard’s  theorem  is  immediate  if  we  can  prove  first  that «!» -«/»*.  Note  that 
if  then  <k*  is  a  distance  function,  L*  — I,  and  L*  is  a  production  technology. 

Before  we  prove  that  1 we  must  first  show  that  if>'  is  well-defined  which  in  this  case 
means  that  we  must  show  that  if  u  >  0,  L then  Lp(u)  is  nonempty.  Since  L(u)  is 
closed  and  convex  with  0  £  L(u)  it  follows  by  the  well-known  strong  separation  theorem  for 
convex  sets  in  Euclidean  space  that  we  can  find  a  p  €  R ",  a  €  R  such  that  0< a  < p-z  for  all 
z  €  L(u).9  Since  L(u)  is  monotonic  it  follows  that  p  €  /?"  .  Dividing  p  by  a  produces  an  ele¬ 
ment  in  Lq  (u).  We  are  now  ready  to  prove  the  theorem. 

Shephard’s  Duality  Theorem 

The  Distance  Function  equals  the  Price  Minimal  Cost  Function,  i.e.,  '  • 

Proof  of  Shephard’s  Duality  Theorem  (Version  1)  i 

If  either  u-0  or  u  >0  but  L(u)-0,  then  the  result  follow?  by  definition.  If  we  let 

Rix)  -  (Xx  |  X  >  0)  be  the  ray  generated  by  x  €  RJ.  .  then  only  two  possibilities  are  left  to  be 

considered: 

(1)  u  >  0,  L(u)  ^  0,  and  R  Cx)  p|  L(w)-0, 

(2)  u  >  0,  L(u)  *  0,  and  R  Or)  pj  L(«)  0. 

Let  us  analyse  case  (1)  first.  By  definition  *(u,x)-0.  Hence  we  must  show  that 
^ ’ (u pc)  - 0;  that  is,  we  must  find  a  p  €  Lp («)  such  that  p  x-fl. 

*  See,  for  example,  Bazarra  and  Shetty  119731,  p.51. 


64 


Let  A  Or)  —  {/ 1  x,  —  0).  Since  L(u)  is  monotonic,  it  is  easy  to  see  that  A  (x)  is  non-empty. 
Set  up  the  following  optimization  problem: 


P:  min  T  y, 

*  €  L(“)<  (A  lx) 


Since  Liu )  -  Eiu )  +  RJ.  it  follows  that  problem  P  is  equivalent  to  problem  P*: 


p':  BjfL,  I  y<- 

1  «  EW;  €  A  U) 


r\ 


E(u)  was  assumed  compact;  the  objective  function  is  continuous  so  we  may  conclude  that  a 
minimrm  is  obtained.  Let  >*  denote  such  a  minimum  and  6  denote  the  objective  function 
value  evaluated  at  y’. 

If  5-0,  then  for  X  large  enough  X-x  > y‘  which  would  imply  that  X-x  €  L(«).  This  in 
turn  would  mean  that  R  be)  meets  L (#)  which  is  not  the  case  we  are  considering.  Hence, 
5  >  0. 

Define  p  €  Rl  by 


j  1  €  A 

0  it  A  tx) 


(3.26) 


Observe  that  p-x-0  and  p-z  >  1  if  z  €  Liu).  Hence  p  €  Lq (u).  with  the  desired  proper¬ 
ties.  We  have  shown  that  ifr*(w,x)  — 0  as  desired. 

Let  us  now  turn  to  case  (2).  In  this  case  we  have  that  0<i>(urx)  <  <».  For  any 
p  €  L q(u)  it  is  immediate  that 


’’"iiujt)  *  1 


0.27) 


since  €  Hu).  Since  (3.27)  holds  for  any  p  €  Lq(u),  by  definition  of  #*(u,x)  (and 

multiplying  through  by  fiujc)  in  (3.27))  we  have  that  *‘(«,x)  >  #(i/,x). 
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If  tp’(u^c)  >  ifi(uyx),  then  by  definition  of  *l>(upc)  (as  a  distance  function)  it  must  be  the 

case  that  — 77 - -  i  L(u).  By  the  Separation  theorem  we  may  find  a  p  €  R"  and  an  a  €  R 

1(1  (u,x) 

such  that 

P'~~rr — r  <  a  <p-z ,  Vz€L(u).  (3;28) 

1 1>  (u,x) 

As  before  we  note  that  p€Rl;  thus  a  >  0.  By  dividing  p  by  a  we  have  that  is  in  Lq  ( u ) 

OL 

with  -2-x  <  1 1>" (.Uyx)  by  (3.28).  By  taking  the  infimum  over  all  p  €  L0 (u)  we  would  obtain  the 
a 

obvious  contradiction  that  ^'(u^c)  <>lt‘(upc).  Hence  (upc)  which  shows  that 

<l/‘(upc)~<Ji(urx)  in  this  case, 
ur  proof  is  complete.* 


3.2.4.  Proof  of  Shephard’s  Theorem:  Version  2 

Let  S  denote  the  set  of  all  closed,  convex,  monotonic  subsets  of  R%  not  containing  the 
origin.  Let  L  €  S.  Let  R  (x)  denote  the  ray  generated  by  an  x  €  Rl  ,  i.e.,  /?  (x)  -  {X-x  I  X  >  0). 
If  x  €  Rl ,  then  either  R  (x)  Q  L  -0  or  R  (x)  p)  L  *  0.  If  R  ix)  P)  L  *  0,  then  there 

exists  a  well-defined  constant  <i»(x)  such  that  -f  -r  is  on  the  boundary  of  L  (relative  loRl). 

*(x) 

_ It R  (x)  L  —0,  then  define  ifrOc)  «— 

Now  pick  an  x  €  Rl  and  suppose  R(x)f)Li*f8.  Let  «  >  0  and  let 7  denote  the  vector 

(*,« . a)  €  Rl .  Since  *  is  continuous  (Proposition  (3.4.8),  Appendix)  then  for  some  y  >  0 

(to  be  determined  later)  we  can  find  an  <  >  0  so  that 


*±I  f  a  ) 
*<x+?)  "V*)' 


(3.29) 


(BA-—--)  denotes  the  open  ball  of  radius  y  about  -;  r ■?)■  Since  — f  **-  ■  is  strictly  positive 
MU)  *(*)  4»(x+7) 

it  follows  by  an  application  of  the  separation  theorem  for  convex  subsets  of  Rl  that  a 
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Y 

supporting  hyperplane  for  L  at  — - - -  exists.  That  is,  there  exists  a  p  €  R% ,  p  s*  0  such  that 

♦  Or +7) 


0<p- <pz,  Vz€I, 
♦  0r  +  «) 


(3.30) 


Re-normalizing  we  may  assume  that  p- 


x+7 
♦  Or +7) 


—  1.  It  follows  therefore  that 


,  _  X,+t  X,+«  I  _ 

i_z a-  >  z  p, • -- > I  z  h 

,  ♦  (*  +  «)  ♦  Or  +  e)  |(i u  > o) 


min  x.  +  « 

1<  f  jr,  >  0| 


♦  Or +7) 


Let 


Then  we  have  that 


/Or,.) - t<£±IL- 

I,  +  * 


Z  A  <  fix,*). 

(/U,  >0) 


(3.31) 


In  view  of  (3.29)  and  (3.31)  we  have  that 


x  Jr+7  ^  r  *.  *<+«  ,,  , 

'  +  W*0r+T)  I  *♦(,+«  <  /(X'‘)  T' 


(3.32) 


Since  p-  -  -  1  rearranging  (3.32)  we  obtain  that 

♦  Or +7) 


p  x  <  ♦(*)  +  ♦  0c)/ (x  ,e)y. 


(3.33) 


Now  let 


//^-{(p.a)  €  RlxR++  |  3x  €  L  such  that  a~p-x  <p-r ,  Vr  €  1} . 


Since  y  is  arbitrary  and  /0r,c)  is  bounded  for  <  small  enough,  it  is  easy  to  see  that  we  have 
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just  shown  that 

inf  -^■•x<4r(x).  (3.34) 

« 

For  (p,a)  €  H+,  -2-  •  —?■  -.-  >  1  so  that,  as  an  easy  consequence, 
a  4>v x) 

inf  £-x>^(x).  (3.35) 

{(l,a)  €  H*)  « 

Putting  (3.34)  together  with  (3.35)  gives  us 

inf  •^••x  — tf»(x).  (3.36) 

The  interpretation  of  If*  is  precisely  that  it  is  the  collection  of  supporting  hyperplanes  for 
L  which  separate  the  origin  from  L.  By  what  we  have  just  shown  those  points  in  the  boundary 
which  have  supporting  hyperplanes  in  H*  are  dense  in  the  boundary.  Further  we  have  also 
shown  that  given  any  point  2  on  the  boundary  there  is  a  point  y  on  the  boundary  arbitrarily 
close  which  has  a  supporting  hyperplane  (p,a)  €  H*  whose  "value"  at  z  is  arbitrarily  close,  i.e., 

1 .  £..y  >  2..z  $  small.10 

at  a 

We  comment  that  it  is  possible  that  there  may  exist  a  point  on  the  boundary  which  has  no 
supporting'  hyperplane  in  H*  to  support  L  at  this  point.  To  see  this,  we  construct  the  following 

example.  Let  x„  -  (1  -  — ,  l-~,  V4")  for  a  >  1  and  let  L  be  the  smallest  closed,  convex, 

/I  ft  , 

monotonic  subset  of  Rl  containing  the  %,'s.  Since  x„  —  (1,1,0),  it  follows  that  (1,1,0)  €  L. 
We  claim  that  no  (/>,<*)  €  H*  could  exist  which  could  support  L  at  (1,1,0).  If  there  were  such  a 
(p,a)  €  H+,  then  it  is  easy  to  se?  that  (1)  at  least  one  of  p\  or  pi  must  be  positive,  and  (2) 
Pi  >  0.11  Now  we  have  that 

y  mf 

19  To  *ee  this,  let  *  ^  0,  and  y  >  0.  Think  of  2  u  ,  y  u  — ■ — — ,  end  (3  u  t.  Review  equations 

♦  u)  *Gr+T) 

0.29)  and  0.32). 

11  If  (1)  were  not  true,  then  o-O;  if  (2)  were  not  true  then  px,  <#  (1,1,0)  for  any  n. 


68 


PX„-(\--)(j>l  +  P2)  +  '/i”P3-(j>l+P2)  +  [iA'P3--^(Pl+P2)}-  (3.37) 

n  n 

Since  —  —  0  as  n  —  »  we  may  find  an  n  large  enough  so  that  te"p3-  —  (P1+P2)  <  0  imply* 
2  n 

ing  that  p-x„  <  p- (1,1,0),  a  contradiction.12  We  remark  that  this  example  dearly  illustrates 
directly  why  the  infimum  is  used  in  (3.36)  as  opposed  to  the  minimum. 

It  is  easy  to  see  that  since  ip, a)  €  H*  if  and  only  if  ,  1)  €  H*  we  may  identify 

a 

with  the  set13 

#*“((p,l)  €  /?*  x/I-m.  I  3x  €  L  such  that  1  <p*x  <p-z ,  Vz  €  L). 

This  last  set  of  course  may  be  identified  with  the  set 

[p  €  R+  |  3x  €  L  such  that  1  <  p  x  < pz ,  Vz  €  L) . 

Thinking  of  L  as  a  level  set  L(u),  this  last  set  is  precisely  Lg(u)  and  ^r(x)  is  really  +(.ujc).  So 
looking  at  (3.36),  we  see  that  we  have  shown  Shephard's  Duality  in  the  case  identified  in  the  previous 
section  as  case  (2). 

So,  by  identifying  L0(w)  with  H*  (which  has  already  been  interpreted)  Shephard’s  Dual¬ 
ity  theorem  in  this  special  case  reduces  to  a  statement  about  supporting  hyperplanes  of  closed, 
convex,  monotonic  subsets  of  /?" .  In  other  words,  loosely  worded,  the  level  set  is  character- 
ired  by  the  boundary  which  is  characterized  by  the  distance  function  on  one  hand  (the  boun¬ 
dary  of  £  is  the  set  of  x  which  attains  the  "distance”  of  1)  and  by  the  price  minimal  cost  func¬ 
tion  on  the  other  hand  (which  approximates  the  boundary  with  supporting  hyperpianes). 

The  reason  why  +  and  **  are  defined  differently  for  special  cases  is  simply  that  one 

desires  to  extend  the  equality  of  +  and  to  the  entire  domain  of  R” .  Let  us  see  how  we  may 
11  This  can  be  done  became  the  terms  pj  and  (pj+pji  ire  both  positive. 

11  This  follows  because  p  can  be  arbitrarily  re-normalized.  Ia  one  case. we  associate  for  the 

point  (-£,1)  €  H‘.  la  the  other  case,  we  associate  for  (p,l)  €  H’  the  point  (p,p-x)  €  H*  where  x  is  a  point 

of  support  in  L  for  p.  Rigorously,  one  could  define  equivalence  desses  of  the  form  l(p,l)l  where 
(pxrK(q^)  if  and  only  if  3r>0  such  that 


motivate  their  respective  definitions  from  this  perspective. 

First,  suppose  u-0,  thus  L (u)  —  R+.  The  idea  of  +(Uyx)  is  that  — .  *  should  be  a 

point  on  the  boundary  of  L(u)  on  the  ray  generated  by  x  "closest*  to  the  origin.  Since  0€  L(0) 
a  natural  choice  for  ^(O^x)  is  to  take  it  to  be  equal  to  ».  Hence,  we  take  »fr*(0,jr)-oo. 

Now  suppose  x  6  Rl ,  L(u)  ^0  but  R  Or)  f  \  L(u)-0  (case  (1)).  From  pr  perties  on 
convex  subsets  of  Rl  we  showed  in  the  previous  section  that  there  is  a  €  Rl  such  that 
min{p-z  | z  €  L(i/)|- 1  and  p\x-0;  that  is,  i|» ’ (« ) — 0  in  this  case.  Hence,  we  take  tittyx)  to 
be  0  in  this  case. 

Finally,  suppose  L(u)  — 0.  In  short,  none  of  the  functions  are  well-defined  for  this  case. 
Since  L(u)-0  implies  that  RCx)  L(u)-0  to  be  consistent;  we  therefore  take  both  $  and 
«k*  to  be  0  in  this  case. 

The  upshot  of  this  development  and  proof  of  Version  2  is  that  Shephard’s  Duality 
Theorem  may  be  motivated  by  studying  properties  of  supporting  hyperplanes  of  a  particular 
class  of  convex  subsets.14  One  does  utilize  heavily  the  simple  and  useful  characterization  of  the 
boundary  by  the  distance  function.  One  obtains  the  Duality  Statement  as  stated  by  Shephard 
when  one  interprets  the  functions  in  an  economic  light. 


3.2.5.  Extensions  to  the  Dynamic  Case 

To  extend  the  proof,  as  given  in  Version  1,  of  Shephard’s  Duality  Statement  for  the  finite 
horizon  Dynamic  Case  we  make  some  modifications.  Instead  of  defining  A  Or)  -  {/ 1  x,  -0)  one 
defines  A,  (x) — {r  I  x,  (t )  —  0} .  Instead  of  solving  the  optimization  problem  as  given  in  the  proof 
of  case  (1)  we  re-write  it  as 


min  TJ{yi,4}dM- 


14  We  remirk  that  we  'proved*  tint  0-0*  m  this  version  by  extending  the  definitions  of  0  tod  0 
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After  making  the  appropriate  modifications  for  the  statement  of  Shephard's  Duality  in  the 
Dynamic  (Function  Space)  case  the  proof  of  Version  1  still  applies  because 

(a)  LN(u)  is  (weak-star)  closed  (see  Theorem  (3.1.1)). 

(b)  EN(«)is  (weak-star)  compact  (see  Proposition  (3.1.6)). 

(c)  LN(u)  C  ENOJT  +  (m»"  (see  Theorem  (3.1.1)). 

(d)  Separation  Theorem  for  Locally  Convex  Topological  Vector  Spaces  applies. 

(We  only  proved  (a) -(c)  under  the  assumption  of  a  finite  horizon.) 

We  make  two  final  remarks.  First,  to  extend  the  proof  of  Shephard's  Duality  Statement, 
as  given  in  Version  1,  for  the  infinite  horizon  case  additional  hypotheses  would  have  to  be 
assumed  to  insure  that  LN(u)  is  (weak-star)  compact.  The  axiomatic  system  presented  in  Sec¬ 
tion  2.2  is  not  strong  enough  to  guarantee  this  property.  Second,  it  is  not  possible  to  extend 
the  proof  of  Shephard’s  Duality  Statement,  as  given  in  Version  2,  because  our  counter-example 
(Section  3.1.3)  shows  that  strictly  positive  points  need  not  be  supportable. 


3.3.  On  Laws  of  Diminishing  Returns 

In  this  section  we  show  how  our  general  framework  of  a  production  process  may  be  used 
to  deduce  two  versions  of  the  Law  of  Diminishing  Returns  as  formulated  by  Shephard.  By 
proving  such  laws  from  the  axiomatic  framework,  we  reduce  the  question  of  whether  such  laws 
are  'true*  to  the  question  of  whether  the  axiomatic  framework  is  reasonable.  In  this  sence  we 
hope  to  get  *at  the  root"  of  such  laws. 

The  two  general  types  of  laws  of  retum'posed  by  Shephard  and  considered  here  are:  (1) 
Laws  of  Return  for  bounded  input  rates  of  Essential  factors,  and  (2)  Laws  of  Return  for 
bounded  intervals.  The  first  type  is  defined  and  proved  in  Section  3.3.1.  The  second  type  is 
defined  and  proved  in  Section  3.3.2. 

3.3:1.  Laws  of  Return  for  Bounded  Input  Rates  of  Essential  Factors 

Taken  from  Shephard  and  Fare  [1980],  p.98  we  describe  this  law  in  their  own  words: 

"A  law  of  return  so  expressed  for  the  static  model  of  production  is  one  of  a  law  of 
bounded  output  rate.  It  is  suggested  for  input  and  output  rates  which  are  not  con- 
,  stant,  i.e.,  for  the  dynamic  structure  of  production,  that  a  law  of  bounded  output 
rate  may  hold,  i.e.,  if  time  histories  for  essential  factors  are  subject  to  an  upper 
bound  on  input  rate,  the  related  output  rate  histories  will  be  bounded  in  some  way 
under  unlimited  increase  in  the  maximal  time  rate  of  the  input  rate  histories  of  the 
other  factors.’ 

We  prove  this  statement,  solely  in  words,  from  the  axiomatic  framework  for  the  finite  horizon 
case. 

Proof  of  the  First  Type  of  the  Law  of  Diminishing  Returns 

Suppose  the  vector  of  system  exogenous  inputs  is  norm  bounded.  As  an  easy  consequence 
it  follows  that  the  cumulative  amounts  of  the  applications  of  exogenous  resources  to  any 
activity  must  be  bounded  in  magnitude  (in  the  finite  horizon).  By  Axiom  2  on  the  flows  of 


goods  and  services  it  follows  that  the  flow  types  associated  with  the  applications  of  exogenous 
resources  is  north  bounded.  By  Axiom  8  bn  the  activity  Production  Functions  this  in  turn 
means  that  the  outputs  realized  through  the  production  process  for  any  activity  must  be  norm 
bounded.  As  an  easy  consequence  of  this  fact  the  cumulative  amounts  of  the  outputs 
transferred  to  the  other  activities  from  a  given  activity  must  also  be  bounded  in  magnitude. 
This  means  that  the  possible  cumulative  amounts  of  final  outputs  obtainable  through  produc¬ 
tion  is  bounded  in  magnitude  hence  in  norm  by  Axiom  2. 

If  we  restrict  attention  to  just  bounding  in  norm  an  essential  set  of  system  exogenous 
inputs,  then  in  view  of  Axiom  8  the  comments  above  still  apply.  Thus  regardless  of  the  magni¬ 
tude  of  the  other  types  of  ‘nonessentiaT  inputs  the  output  rate  would  be  bounded.  And  this  is 
the  statement  of  the  first  version  of  the  Law  (in  the  finite  horizon  case). 

3.3.2.  Laws  of  Return  for  Bounded  Intervals  of  Essential  Factor  Application 

Taken  from  Shephard  and  Fare  [1980],  p.100  we  describe  this  law  in  their  own  words: 

'The  time  spans  over  which  essential  input  rates  may  be  or  are  applied  positively 
need  not  be  infinite,  that  is  the  support  of  an  input  may  be  bounded,  and 
unbounded  time  substitutions  for  resources  may  not  be  permitted.  Then  the  ques¬ 
tion  arises  how  outputs  may  be  limited  by  limitations  on  the  intervals  of  time  over 
which  essential  factors  are  applied.  Propositions  of  this  type  are  laws  of  return  for 
bounded  intervals  of  application  of  essential  factors.' 

To  prove  this  law,  we  will  first  prove  a  stronger  statement. 

Proposition  (3.3.1) 

If  x  €  LN(u)  such  that  for  some  T  >  0,  [xJ  >  0)  C  10,7*1  for  all  J,  then  there  exists  some 
f*  such  that  {u*  >  0}  C  lO.f’J  for  all  k.1 
1  h  is  sufficient  to  prove  the  Law  for  the  Essential  Set  of  inputs  constating  of  <df  inputs. 
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To  pave  the  way  for  a  'verbal*  proof  of  Proposition  (3.3.1),  we  prove  a  simple  proposition 
below.  For  this  proposition,  let  L  denote  a  space  associated  with  a  particular  flow  type.  Let  e 
denote  the  'lower  bound'  on  the  flows  in  L  as  dictated  by  Axiom  3.  Let  5  denote  the  'minimal 
spacing"  between  the  set-up  times  of  a  flow  in  L  as  dictated  by  Axiom  4. 

Lemma  (3.3.2) 

T  T 

If  x,z  €  L  such  that  J* xdp.  <  J*  zdn  for  all  r  €,R+  and  [z  >  0}  C  [0,7"1  for  some  T,  then 

o  o 

there  exists  a  T‘  such  that  {*  >  0}  C  [0,7"*J. 

Proof  of  Lemma  (3.3.2) 

If  the  cardinality  of  the  number  of  set-up  times  in  Six)  (see  2.19)  is  finite,  then  it  is  easy 
to  see  that  T'  —  sup{t  \  t  €  S' Or))  has  the  desired  properties. 

Suppose  then  that  the  cardinality  of  Six)  is  infinite.  By  Axiom  4  we  may  index  S (x)  by  a 
set  of  time  points  i  — 1,2 . such  that 

Os/0<i|<r2<  •  •  •  such  that  |r,  -  r,_, )  >  & ,  Vi .  (3.38) 

By  the  definition  of  the  t/s  it  »'4ows  that  either 

f  xdfi-0  or  f  xdp.  >  min{«(tj+i -/,),«), 


the  latter  inequality  due  to  Axiom  3.2  But  then  it  follows  that 

/  xdfi  >  min{«-8,e)  if  J  xd(i  *  0.  (3.39) 

I',  .',♦!>  W 


1  We  must  account  for  the  possibility  that  L  is  an  event-based  Dow  type.  If  this  were  not  the  case,  then 
the  r/s  would  have  to  be  titne-trid  points;  so  J  which  is  either  0  or  at  least «. 


1  *-*  mm*  Ve  *  a  *  a  *  » ’*  *  *  *  *  ,  •  ***  %  %  .  *_*  *  *  •  *  w  *  *  *  *  VV  * 


Now  let  N-  {/ 1  J*  xdfi^O).  Since  the  cardinality  of  S(x)  is  infinite  it  follows  that 
the  cardinality  of  N  is  infinite.  But  this  could  not  be  the  case  because  we  would  have  that 

+°°  >  J*  zdp.  ^  f  xdp.  >  X  /  “X  /  >  X  nhn{«‘6»«)  “  +“ 

0  )t+  (-0  (/,./, +1)  /  €  N  I»,./j+|)  UN 

So  the  proof  of  Lemma  (3.3.2)  is  complete. ■ 

Proof  of  Proposition  (3.3.1) 

Suppose  the  concf'tions  of  the  theorem  held  Find  a  feasible  flow  for  x  to  support  output 
level  u.  By  Lemma  (3.3.2),  we  may  easily  deduce  that  the  supports  of  the  flow  type  associated 
with  the  applications  of  exogenous  resources  to  the  activities  are  finite  By  Axiom  11,  it  follows 
that  the  realized  output  of  production  for  each  activity  is  restricted  to  this  finite  horizon.  By 
Axiom  12,  only  a  finite  amount  of  output  may  be  realized  through  production  in  this  finite  hor¬ 
izon;  hence  by  Lemma  (3.3.2)  the  supports  of  the  flow  types  associated  with  the  intermediate 
product  flows  are  finite.  By  a  final  application  of  Lemma  (3.3.2)  the  supports  of  the  flows  types 
associated  with  the  final  output  variables  must  be  finite.  That  is,  eventually  output  must  stop  and 
this  is  what  we  wanted  to  show.* 

1  i  t 

Proof  of  the  Second  Type  of  the  Law  of  Diminishing  Returns 

Since  the  supports  of  the  flow  types  associated  with  the  find  output  variables  is  finite,  it 
follows  immediately  by  Axiom  2  that  such  flows  must  be  norm  bounded.  This  is  Shephard’s 
version  of  the  Second  Type  of  the  Law  of  Diminishing  Returns. 


3.4.  Appendix 

In  this  appendix,  we  prove  the  technical  propositions  used  in  the  main  body  of  the  text. 
As  notation,  let  Lf  denote  the  subset  of  L“  (p)  which  satisfies  Axiom  2,3,  and  4  for  an 
event-based  flow  type,  let  denote  the  subset  of  L+  (p)  which  satisfies  Axioms  2,3,  and  4 
for  a  continuous  flow  type,  and  let  LE»Li  (^\  I/,  Ic  —  L£  p|  L*-  It  is  under¬ 

stood  that  these  symbols  are  generic  and  stand  for  a  particular  flow  type.  Finally,  throughout 
this  section,  L?  (p)  is  topologized  with  the  relative  weak-star  topology. 


Proposition  (3.4.1) 

If  {xj  and  lyj  are  nets  in  L~  (p)  such  that  xa  < ya  for  all  a  and  if  x„  — »  x,  ya  —  y , 
then  x  <>'. 


Proof  of  Proposition  (3.4.1) 


Suppose  this  is  not  true.  This  implies  that  we  find  a  set  A  with  positive,  finite  measure 
and  an  <  >  0  such  that  the  function  x  exceeds  the  function  y  by  c  on  the  set  A .  Let  p  —  1^ . 
Clearly,  p  €  Lx.  Let  denote  the  continuous  linear  functional  induced  from  p.  By  the  con- 
ti  iuity  of  p,  +,txa)  —  *,Cx),  +>(ya)  —  ifr,^),  and  +pbcn)  <  +p(ya)  for  all  a.  Hence, 
*,(x)  <  ♦,(y).  But  by  definition  of  p ,  4»,(x)-4>,(v)>«mC4)>0.  The  contradiction  is 
ruadily  apparent.* 


Proposition  (3.4.2) 

Le  and  Lc  »re  each  closed  in  L?  0*)- 

Before  we  prove  Proposition  (3.4.2),  we  first  show  that  neither  L£  nor  Lc  is  by  itself 
dosed  as  the  following  two  counter-examples  show. 


n*er-Exaaple  to  the  Closure  of 
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Let  c  denote  the  lower  bound  described  in  Axiom  3.  Let  the  time  grid  T  (as  described  in 
Axiom  1)  be  the  set  of  natural  numbers.  Define  for  n  ^  1, 


<-l  *  2i I’m 


It  may  be  verified  that 


x2n  1(0,1)  - 

Of  course,  the  function  yl^i)  violates  Axiom  3. 

Counter-Example  to  the  Closure  of  Lc 

Let  8  denote  the  minimal  distance  between  set-up  times  as  described  in  Axiom  4.  Fix 

h  >  0  and  choose  K  >  •£.  Again,  let  T  be  the  set  of  natural  numbers.  Define  for  n  >  1, 
o 


K  K  1 

i-l  X  IK'  k’  i- 1"  X 


J _ ]_ 

K  2 K'  K 


The  set  of  set-up  times  for  each  x„,  S(jc„),  is  clearly  empty;  thus,  each  x„  satisfies  Axiom  4. 
But,  clearly,  . 


<-l  X  2X’X 


which  has  K  set-up  times  in  the  horizon  [0,/r  ].  Thus,  J  1  , _ j_  .  violates  Axiom  4. 

I-I  X  U'X 


Proof  of  Proposition  (3.4.2) 

In  light  of  the  two  counter-examples,  we  will  prove  Proposition  (3.4.2)  by  first  proving 
that  (and  L}  )  is  closed  (Proposition  3  4.3)  and  then  prove  that  Li  f\Lc  (and  Li  f\  Li) 
is  closed  (Proposition  3.4.4). 
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Proposition  (3.4.3) 

and  L}  are  each  closed  in  LX  ip.)-  , 

Proof  of  Proposition  (3.4.3) 

By  the  continuity  of  parameters  {/»j)  used  to  define  Axiom  2  for  an  event-based  flow  type, 
it  is  immediate  that  Ll  is  closed  in  LX  (m)-  To  show  the  closure  of  Lt2,  we  will  show  that  the 
complement  of  L}  in  LX  (p)  is  open. 

Pii  a  2  in  the  complement  of  I<?  in  LX  (a*)-  It  follows  by  definition  of  Axiom  2  that  for 
some  A  €  R+  and  Ik  that 

but  llr-l/JL  >  gk(A).  (3.40) 

Let  8-llz-ljJL  — gk(A).  Since  gk  is  continuous,  we  may  find  an  e  >  0  such  that 

|g*(2J)-g*G4)|  <  4  whenever  \B-A\  <  «.  (3.41) 

4 

Now  define 

a  -  (/  pdp)j 

v  -  (z + tfo>,<*))  n<*+ iV('i/^,c))  n  /  r  o*) 

where  for  y  €  *.'(m),  B  >  0 

Z  +  N(qji)  -  |y  €  17  Qi)l  lJ*fl(z-,y)dMl  <  0 
( N  is  a  neighborhood  about  z  in  this  topology.) 
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Pick  any  y  fc  N.  Let  3- f  yllkdn.  Since  y  €  z  +  N(p, a)  it  is  easy  to  verify  that 

Itrl^lL  >  |. 

By  (3.40)  and  (3.41),  we  then  have  that 

Hyl/JL>llrl-.IL-4"&^>+T*  (3.42) 

*  *4  i 

There  are  two  cases  two  consider. 

Case  1:  \B—A  |  <  e. 

Since  —  A  |  <  «,  by  (3.41)  and  (3.42),  it  is  immediate  that 

Hrl4IL  >  **(20  +  -| 

which  implies  that  y  is  in  the  complement  of  in  LT  0*)- 

Case  2:  \B-A  I  >«. 

Since  y  is  also  in  z-f  Af(l/4,e),  we  have  that 

/ z  l,kdn-«  <  J y  l,kd(i  <  J* <><+«• 

Hence,  B  <  A  +«.  Since  \B  -  A  I  >  «,  it  follows  that  B  <  A  By  monotonicity  of  gk ,  we 
then  have  that 

gk(JB)  <  &C4  -«)  <  itlA).  (3.43) 

It  is  immediate  by  (3.42)  and  (3.43)  that  y  is  in  the  complement  of  L £  in  L?  0*)- 

In  either  case  1  or  case  2  we  have  shown  that  Jy  €  S,  then  y  is  in  the  complement  of 
in  L?  Ox).  Thus,  the  complement  of  Q  in  L~  Ox)  is  open  proving  the  desired  result. ■ 
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Proposition  (3.4.4) 

L\  O  Le  and  Pl  are  each  closed  in  L?  ip.). 

Proof  of  Proposition  (3.4.4) 

By  the  definition  of  an  event-based,  flow  type  and  the  properties  imposed  on  the  lime  grid 
T,  it  is  easily  verified  that  Lj!  and  1/  are  each  closed  in  17  ip)-  Hence,  so  is  Li  P|  I/. 

To  facilitate  the  proof  of  closure  of  P|  L*,  we  first  prove  a  simple  lemma. 

Lemma  (3.4.5) 

Let  x  €  Lc  and  let  S(x)  denote  the  set  of  set-up  times  for  x  (see  2.19).  If  7—  (a ,b]  such 
thatx{7  P)  (x-O)J  >  0  and  X{/  P\  (x  >  0)}  >  0,  then  S(x)  P)  7^0. 

Proof  of  Lemma  (3.4.5) 

If  a  €  SC*),  uen  the  result  follows  immediately.  So,  suppose  a  i  SGr).  By  definition  of 
SC*),  there  must  be  an  interval  of  the  form  (a,?),  c>a ,  such  that  either 
X(ta,cl  P|  Cx- 0)}-0  or  X{Ia,c]  p|  (x  >  0))-0.  Without  loss  of  generality  assume  that 
Xlla.cl  P|  Or  >  0))-0.  Let  t'-sup[d  €  /  I  X{la,dl  P)  Or  >  0)1-0). 

We  now  claim  that  r*  €  SC*).  If  t'  t  Six),  then  by  definition  of  SCx)  we  could  either 
find  a  /“  >  r*.  t"  €  la,bl,  which  has  the  property  that  X([d,/**1  p|  Gr  >0))-0  or  t’  must 
equal  b.  By  definition  of  t\  it  must  be  the  second  case,  i.e.,  f*  — 6.  But  then  this  means  that 
X{7  Pj  C*  >  0)}—0  which  violates  our  original  assumption.  The  restilt  follows.* 

We  resume  the  proof  of  Proposition  (3.4.4).  Let  (xj  be  a  net  in  Lj  f)  Lj  converging  to 
some  x.  Suppose  xiL$.  Then  (0  <  x  <  a)  has  positive  Lebesgue  measure.  Since  Lebesgue 
measure  is  regular,  it  follows  that  there  is  a  finite  interval  la,blc  (0<x<«).  Choose  a 
t  €  (a,b).  Let  8  denote  the  minimal  distance  between  set-up  times  u  described  in  Axiom  4. 

Choose  «€  (0,-j)  so  that  if  7|- lr-~,r  +  -jl  and  72-lf+«,t  +  2«)  then  7,  C  («,b),y-l,2. 
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We  make  the  following  claim:  for  some  0,  S(x„ )  P)  Ij  *  0,  y  — 1,2,  for  all  a  >  If  this 
were  not  true,  then  we  could  extract  a  countable  sub-net  {xj  such  that  S(x„)  P)  /,  «*0, 
y-  1,2.  By  Lemma  (3.4.5),  it  follows  that  for  each  n  either  y  —  1,2.  This  in  turn  would  imply 
that  either  or  J*x,,-1  idu'Z  «X(/y)  >  f  x-l/dft  >  0,  ./— 1,2.  For  ^-1;, 

j  - 1,2,  we  would  then  have  that  x„  does  not  converge  to  or -an  obvious  contradiction. 

So  we  have  shown  that  for  some  0,  S Car.)  p|  /,  *  0,  y- 1,2,  for  all «  >  0.  By  construc¬ 
tion  of  each  /,  it  is  immediate  then  that,  for  a  >  0,  x.  violates  Axiom  4.  And  this  is  a  contrad¬ 
iction.  Hence,  x  must  be  in  !<?. 

Suppose  xiL*.  Then  we  may  find  a  rt  and  i2  each  in  S(x)  such  that  |r2-fil<6. 

Choose  «€(0,- — ^ — ~)  so  that  if  /,-  y«l,2,  then  /,  P) /2-0.  By 

definition  of  S(x),  we  have  that  X{/y  (x>0)>0,  /- 1,2.  Since  x„  —  x,  we  then  have 
that,  for  all  a, 

f*>\ n u-onrfM  —  / x\ n u-wi •  ,JmW 
/*.1u/n<*>on^-/x-,(/)nt»>o)i  PI  <*>o)l  >  o,  y-1,2.  0.45) 

(We  are  using  the  fact  that  x  €  in  (3.45).)  By  (3.44)  and  the  fact  that  each  x„  €  I<?  it  fol¬ 
lows  that  for  some  0y,  >- 1,2,  X|/y  P|  (x.-0)}>0  for  all  «>0y,  y-1,2.  Similarly,  by 
(3.45)  and  the  fact  that  each  x.  €  L^,  it  follows  that  for  some  yj%  J~  1,2, 
x(/y  p|  (x.  >  0)}  >  0  for  all  a  >  yJt  J  — 1,2.  Thus,  by  Lemma  (3.4.5)  and  the  above,  we  have 
that,  for  all  a  >  max  (0;,yyJ,  S(x„)  f)  lt  a*  0,  >—1,2.  This  implies  that  xm  (  Lc  which  is  an 

obvious  contradiction.  Hence,  x  must  be  in  Lc- 

The  above  argument  shows  that  x  €  L<*  p)  The  result  follows.* 

Proposition  (3.4.4) 
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Let  {*„),  W  be  nets  in  Lc  such  that,  for  all  a,  MCxa  >  0)  P|  (ya  >  )}-0.  If  xa  —  x 
and  y„  -*  y ,  then  X{(x  >  0)  p|  (y  >  0)}  - 0. 


In  the  discussion  on  Axiom  5  (closure  of  the  domain),  Section  2.2.2. 1  we  remarked  that 
there  may  be  additional  constraints  linking  the  domains  of  the  inputs  applied  in  production. 
One  example  we  gave  is  when  an  activity  utilizes  one  machine  to  produce  two  similar  types  of 
products.  For  this  example,  the  functions  which  define  the  rate  of  machine  hours  applied  to 
each  type  of  product  are  linked  in  that  both  functions  cannot  be  positive  at  the  same  time.  It  is 
desirable  to  show  that  this  property  preserves  closure  so  as  to  give  a  plausible  basis  for  accept¬ 
ing  Axiom  5. 

Before  proving  Proposition  (3.4.6),  we  remark  that  if  the  nets  (xj  and  (yj  were  only 
assumed  to  belong  to  L?  (p)  then  Proposition  (3.4.6)  would  be  false.  To  see  this,  define  for 
n>l 


*2-  *  1  l.±_ ±  1) 
<-l  U’h 


yim  •  l<o,n-*i»  • 


It  may  be  verified  that 


*2»  — '  1«M>  •O**  yu  "*  1«U> 

and  thus  clearly  the  nets  (x}„),  {yjJ  satisfy  the  conditions  of  Proposition  (3.4.6)  but  violate  the 
conclusion. 


Prwf  af  Propaaiti—  (3.4.0 

Suppose  X'tr  >  0)  f)  (y  >  0)}  >  0.  Find  a  finite  interval  (a,6l  C  {Or  >  0)  ft  (y  >  0)). 
Choose  a  t  €  (e  ,6).  Let  8  denote  the  minimal  distance  between  set-up  times  as  described  in 
Axiom  4.  Choose  <  €  (0,-j)  so  that  if  +  /2»  {ir+t,r  +  ?*J  then  Ij  C  (o,6l. 
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;-u. 

Since  xa  — *  x ,  x„  €  L$  for  each  a,  and  X{(x  >  0)  P|  /,)  >  0,  /— 1,2,  it  follows  easily  that 
eventually,  for  some  fy,  X{(x„  >  0)  P|  7, }  >  0,  Va  >  /3 j,  j— 1,2.  Since  ya  —  y,  y„  €  Q  for 
each  a,  and  X{(y„  >  0)  P)  Jj)  >0  J-  1,2,  it  follows  again  that  eventually,  for  some  6j,  j-  1,2, 
X((y0  >  6)  f)  /;}  >  0,  Va  >  6j,  j  — 1,2.  Since,  for  all  a,  X{(xa>  0)  f)  (ya>  0)}-0  we  have 
that  for  a  > 8 ;,  j— 1,2,  X{(xa  > 0)  p|  Ij)  >  0.  By  Lemma  (3.4.5),  S(xa)  j- 1,2. 

But  this  implies  that  for  a  >  max  {/9,,8;),  xa  i  L£  by  our  construction  of  Ij,  j- 1,2.  This  con¬ 
tradiction  proves  the  result.* 

Proposition  (3.4.7) 

Let  F:L?  (p)  “*  IT  (p)  such  that  for  each  x  €  Z,"  <p),  f  is  bounded  (in  nomi)  in  a 
neighborhood  of  x.  Assume  further  that  for  each  A  €  /?+,  7  satisfies  the  following  property:  V 

«  >0, 

38  >  0  such  that  |  f  (F(x)-F(y))dftl  <  •  whenever  if  (jc-y)dn  |  <  8. 

(DM  (DM 

Then  F  is  weak-star  continuous. 

Proof  of  Proposition  (3.4.7) 

Pick  an  x  €  ZJ*  (p).  Let  U Or)  denote  the  neighborhood  about  x  in  which  F  is  bounded. 
Let  D(x)  denote  the  bound.  Let  F(x)  +  N(p* . p*;«)  denote  a  basic  open  neighborhood 

about  f(x).  Let  8 —^Bix).  Find  A  large  enough  so  that  if  /-(A,oo)  then  llp'  l/ll|  <8  for 
/  — 1,2, Find  integrate,  simple  functions  Z-1,2.  . . .  such  that  llp'-f'llj  <8 
for  i-l,2 . a.  Finally,  let  p- min  • 

We  may  find  a  8>0  that  the  property  on  F  bolds  when  *-p.  Let 
P(x)-p  €  (x  +  JVd/,,#))  n  n  /.r  0*>-  For  any  p,,/- 1,2...,*,  we  have  that 
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llp'(FOc)  — FO')>!l,  <  llp'(F(x)-F(y))  l,ll1  +  2B(x)llp'-l/ll, 


<  lip'  - q‘W,\\F(x)  —  F(y)lloo  +  ||«'ll.-ll(F(x)-.F(y)M,llI  +  IBix) llp'1,11, 


+  “  t*><«S77T>  -  *  • 


3II«'II. 


6B{x) 


Hence,  we  have  found  a  neighborhood  about  x,  K(x),  such  that  if  y  €  K(x), 
Fiy)  €  F(x)  + . p";«).  Thus,  F  is  weak-star  continuous.* 


Proposition  (3.4.8) 

,  Let  I  denote  a  closed,  monotonic  but  not  necessarily  convex  proper  subset  of  it"  which  is 
non-empty.  Let  (jc)  denote  the  distance  function  (dropping  the  u  in  +(ujc)).  +  is  continu¬ 
ous  on  it”.1 

Before  we  prove  this  proposition,  we  remark  that  the  proof  of  Shephard’s  Duality  given  in 
Section  (3.2.4)  utilized  the  fact  that  the  distance  function  is  continuous.  Shephard  11970] 
proves  that  4>(u,x)  is  continuous  on  x  R1  by  first  citing  a  theorem  that  states  that  a  con¬ 
vex  function  defined  on  a  convex  open  subset  of  R*  is  continous  on  this  open  subset  and  then 
proving  that  +  is  both  upper  and  lower  semi -continous  on  the  boundary.2  Our  proof  below  is 
simpler  and  more  direct;  furthermore,  we  do  not  assume  that  the  level  set  is  convex. 

Proof  of  Proposition  (3.4.8) 

To  prove  this  proposition,  we  first  prove  a  very  simple  lemma. 

Lemma  (3.4.9) 

1  Since  0  is  constant  and  hence  continuous  when  L-R  J  (corresponding  to  the  case  *«*0)  or  when  L  •• 
we  omit  these  cases. 

]  Proposition  16,  Shephard  (19701.  Berge  (19631,  p.  193. 


T 
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if  or  >>  then  >  ^(y). 

Proof  of  Lemma  (3.4.9) 

If  this  were  not  true,  then  for  some  x  >y  +(x)<  ♦(y).  As  ♦O')  >  0  this  means  that 
R  O')  P)  L  *  Z.  Gearly  this  implies  that  R  (x)  P|  L  *  0  since  x&y.  By  monotonicity  of  L 

x  11 

and  by  definition  of  ♦  we  have  that  r  €  I  with  —7- r  <  ~-r—r.  This  contradicts  the 

♦  O')  ♦  O')  ♦  (*) 

definition  of  ♦U).  The  result  follows.* 

We  now  proceed  with  the  proof  of  Proposition  (3.4.8).  First  we  show  that  ♦  is  continu¬ 
ous  on  {or  I  ♦0r)-*0).  Let  x  be  such  that  ♦  0r)  —  0.  Define  y"  byy,**-*,  +  —  for  each  i. 

ft 


Claim  (3.4.10 

♦  O'”)  —  0  as  n  —  00. 


Pntf  ef  Claim  (3.4.10) 

If  this  were  not  true,  then  for  some  <  >  0  we  could  extract  a  subsequence  y*k  such  that 
♦  O'**)  >  c  for  all  4.  By  Lemma  (3.4.9)  it  follows  that,  for  each  *,  <  <  ♦  (y"*)  <  ♦(y*1).  Thus 
(♦O'**))  >>  contained  in  a  compact  subset  and  therefore  has  a  convergent  subsequence.  (We  will 
not  change  notation  for  the  subsequence.) 

Hence  for  some  p  >  0  we  have  that 


(3.46) 


"ft 

But  €  L  for  each  k  and  since  L  is  closed  this  means  that  —  €  L  too.  But  then  R  Or) 

♦O'*)  p 

would  meet  L  implying  that  ♦Or)  >  0,  a  contradiction.  The  result  follows.* 


A  '.V- -V-  A  A  A  A  *.  A/..-*.  A  a;>va 


[»1»J 


if  x,  —  0  then  /?  Oc)  p|  1*0  implies  that  R  (y)  P|  L  *  0.  From  this  we  can  now  say  that 

0  <  ♦  (2"1)  <  <  ♦  Oc)(l-a) ,  V*. 


As  before,  extract  a  convergent  subsequence  so  that  for  some  p  >  0  we  have  that 


♦  Or) ' 


(3.47) 


2  *  w 

Since - — ■  is  in  the  boundary  of  L  for  each  k  (which  is  closed)  then  —  must  be  in  the 

*(2"*)  P 

boundary.  But  by  (3.48)  this  is  clearly  not  the  case.  A  contradiction  is  reached  thus  proving  the 

claim.a 

Let  a  >  0.  By  Claims  (3.4.11)  and  (3.4.12),  we  can  find  n(a),m(a)  so  that 
♦  (y"<a))  <^(x)(l+a)  and  +(zm(a))  >  *(x)(l— a).  Let  tf(a)«*max(n(a),m(a)}.  Let 
h  €  B  |  Or)  p|  R J .  By  Lemma  (3.4.9)  we  have  that 

*<a) 


♦  0r)(l  — a)  <  ♦  (2"<a))  <  ♦  (A)  <  +(yN(a))  <  ♦  0r)(l+a) 

Therefore  |«p  (A )  —  Cx)  |  <<*♦(*).  As  a  is  arbitrary,  it  follows  that  ♦  is  continuous  on 
(x  |  ♦Or)  >  0}.  Hence,  ♦  is  continuous.*. 
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4.  APPLICATION  OF  THE  GENERAL  MODEL  TO  MULTI-PROJECT  RESOURCE- 
USE  PLANNING 

In  this  chapter,  we  use  the  general  model  as  a  tool  to  provide  a  systematic  analysis  of  a 
heuristic  solution  proposed  by  Leachman  and  Boysen  [1983]  for  the  problem  of  multi-project 
resource-use  planning  for  a  multi-project  production  system.  The  problem  is  to  determine  explicit 
resource  allocations  through  time  to  projects  to  insure  that  schedules  are  met.  Our  systematic 
analysis  not  only  provides  a  logical  foundation  for  their  approach  but  also  shows  how  their 
approach  can  be  extended  and  improved.  More  importantly,  the  analyses  carried  out  in  this 
chapter  illustrate  the  value  of  using  a  general  conceptual  framework  of  a  production  system  to 
evaluate  proposed  heuristic  solutions  to  production  planning  problems. 

4.0.  Introduction 

A  multi-project  production  system  0  is  a  production  system  comprised  of  a  number  of 
.  single-project  production  systems  GUG2 . Gk  each  utilizing  the  same  set  of  system  exo¬ 

genous  inputs.  In  these  organizations,  project  managers  are  responsible  for  keeping  projects  on 
schedule  and  within  budget.  Without  an  effective  method  for  allocating  scarce  resources  to  the 
projects  frequent  project  delays  ensue.  Such  delays  could  be  avoided  if  an  effective  method  for 
multi-project  resource-use  planning  were  available;  that  is,  a  method  for  determining  explicit 
resource  allocations  to  projects  through  time  to  insure  that  schedules  were  met. 

Traditionally,  decision  support  systems  for  multi-project  planning  develop  resource- 
constrained  schedules  of  the  activities  within  each  single-project  production  system  by  treating 
the  multi-project  system  as  if  it  were  one,  large  single-project  system.1  In  our  model  of  a 
single-project  production  system  (see  Section  2.3.4),  once  the  schedules  for  the  activities  are 
determined  the  resource  allocations  for  the  activities  and  hence  the  resource  allocations  to  the 
projects  are  determined.  The  problem  with  this  approach  is  twofold:  first,  developing  resource- 
constrained  scheduling  of  the  activities  for  large  single-project  production  systems  is 


1  See.  for  example,  Kurtuius  and  Davis  [19821. 


computationally  unattractive;  second,  higher-level  manage'  ent  is  not  responsible  for  detailed 
scheduling  of  individual  activities.  For  similar  reasons,  .t  is  also  unsatisfactory  to  approach  the 
problem  by  determining  resource  allocations  to  the  activities  A  new  approach  is  needed. 

Lcachman  and  Boysen  in  their  1983  paper,  ’An  i  ggregate  Model  for  Multi-Project 
Resource  Allocation,"  present  a  novel  approach  to  solving  the  problem  of  multi-project 
resource-use  planning.2  The  general  idea  is  to  combine  activities  in  each  single-project  produc¬ 
tion  system  into  "aggregate  activities"  and  then  allocate  resources  to  the  aggregate  activities*  It  is 
envisioned  that  the  allocations  of  resources  to  the  aggregate  activities  would  serve  as  resource 
constraints  for  the  resource-constrained  scheduling  to  follow.  For  the  method  to  work  the 
choice  of  allocations  to  aggregate  activities  must  insure  that  the  set  of  schedules  consistent  with 
the  allocations  includes  schedules  which  are  on  time  and  within  budget. 

To  develop  the  general  idea  of  aggregation,  let  G  denote  a  production  network  associated 
with  a  single-project  production  system.  When  activities  are  combined  into  aggregate  activities, 
an  "aggregate  production  network"  G'  derived  through  structural  aggregation  from  G  is  created. 

Definition  (4.1) 

A  network  G'  with  M  nodes  is  said  to  te  derived  through  structural  aggregation  from  a  net¬ 
work  G  with  N  nodes  if  G’  is  acyclic,  directed,  and  with  the  following  property:  there  exists  a 

surjective  map  {1,2 . JV}— >{1,2 . M)  such  that  arc  Uj)  is  in  G'  if  and  only  if  there 

exists  an  arc  (*,/)  in  G  with  fc€tj»-,(/),  /€gf  *(/). 

For  example.  Figures  (4-la)  and  (4-lb)  show  two  detailed  subnetworks  and  G2  associated 
with  a  single-project  production  system.  The  circles  indicate  which  activities  have  been  com¬ 
bined  into  aggregate  activities.  The  resulting  aggregate  subnetworks  G\  and  G'2  of  the  aggre¬ 
gate  production  networks  derived  through  a  structural  aggregation  .  re  shown  in  Figures  (4-lc) 
and  (4*ld).  Note  that  Figures  (4-lc)  and  (4-ld)  are  identical.  This  example  shows  that  two 


2  See  also  Leachman  and  Boysen  (19821. 


(a)  Detailed  Subnetwork  G1  (b)  Detailed  Subnetwork 


(c)  Aggregate  Subnetwork  G^  (d)  Aggregate  Subnetwork  G 

FIGURE  (4-1) 

EXAMPLES  OF  STRUCTURE  UNDERLYING  AGGREGATE  SUBNETWORKS 


Aggregate  production  networks  may  be  identical  but  the  underlying  production  networks  from 
which  they  were  derived  may  be  fundamentally  different. 

For  this  example,  it  is  intuitively  clear  that  the  choice  for  resource  allocations  to  aggregate 
activities  B  and  C  (in  G\  or  <7'2)  should  be  dependent  in  some  way  on  the  choice  for  the 
resource  allocations  to  aggregate  A.  Since  G\  is  fundamentally  different  from  £2>  it  is  also 
intuitively  clear  that  the  dependence  between  choices  for  resource  allocations  of  aggregate 
activities  A,  B ,  and  C  in  G\  is  different  from  the  dependence  between  choices  fct  resource 
allocations  of  aggregate  activities  >4, 2?,  and  C  in  C'2.  Essentially,  Leachman  and  Boysen’s 
method  for  determining  resource  allocations  to  aggregate  activities  to  facilitate  scheduling  of  the 
activities  is  a  method  for  moS'ling  the  dependence  between  the  choices  for  resource  allocations 
to  aggregate  activities.  The  example  hopefully  motivates  why  the  production  network  and  the 
aggregate  production  network  must  both  be  considered  when  developing  a  method  for  modeling 
the  dependence  between  choices  for  resource  allocations  to  aggregate  activities. 

In  their  paper,  Leachman  and  Boysen  give  an  example  of  a  subnetwork  Gx  of  G  and  a 
subnetwork  G\  derived  through  structural  aggregation  from  G\  for  which  it  would  be  difficult 
to  model  the  "dependence  relationships"  between  the  aggregate  activities  in  G\  so  as  to  facili¬ 
tate  the  detailed  scheduling.  Hence,  they  present  only  those  subnetworks  G\  and  G\  for  which 
they  could  model  the  dependence  relationships  using  linear  constraints.  They  then  formulated  a 
linear  program  to  accomplish  the  multi-project  resource-use  planning. 

the  ideas  presented  in  Leachman  and  Boysen’s  paper  are  novel,  innovative,  and  intui¬ 
tively  appealing.  They  recognized  that  dependence  relationships  between  aggregate  activities 
exist  and  that  modeling  these  dependence  relationships  is  intimately  related  to  the  structure  of 
the  subnetworks  at  both  levels.  The  fundamental  problem  with  their  approach  is  that  is  does 
not  present  a  methodology  for  attacking  the  general  problem  of  how  one  should  model  the 
dependence  relationships  for  an  arbitrary  subnetwork.  Furthermore,  certain  modeling  tech¬ 
niques  were  employed  but  were  hot  adequately  justified. 


In  this  chapter  we  use  the  general  model  as  a  tool  to  provide  a  systematic  approach  to 
analyzing  how  one  should  model  the  dependence  relationships.  To  a  lesser  extent,  our  sys¬ 
tematic  approach  provides  a  mathematical  basis  for  accepting  some  of  their  modeling  tech¬ 
niques.  The  approach,  being  structured  and  systematic,  enhances  clarity  and  thus  provides  valu¬ 
able  insight  into  their  approach.  To  a  larger  extent,  our  systematic  approach  provides  the 
means  for  modeling  dependence  relationships  for  a  larger  class  of  subnetworks.  Furthermore, 
for  one  important  class  of  subnetworks,  our  approach  differs  substantially  from  theirs.  For  this 
class  of  subnetworks,  we  feel  our  approach  is  more  sensible. 

The  key  idea  to  developing  a  systematic  approach  to  modeling  the  dependence  relation¬ 
ships  between  aggregate  activities  is  to  realize  that  Leachman  and  Boysen’s  approach  for  multi- 
project  resource-use  planning  is  a  production  planning  technique  for  an  'aggregate’  production ' 
system  (albeit,  a  conceptual  one).  As  mentioned  in  Chapter  1,  to  do  production  planning  for 
an  (aggregate)  production  system  the  input-to-output  transformation  must  be  modeled  for  the 
(aggregate)  production  system.  Our  systematic  approach  begins  by  assuming  that  the  correspon-, 
dence  which  defines  the  input-to-output  transformation  for  the  aggregate  production  system 
satisfies  the  axioms  in  Chapter  2.  This  assumption  implies  that  to  model  the  input-to-output 
transformation  the  flow  types  associated  with  the  aggregate  production  system  must  be  modeled 
(i.e.,  one  must  model  the  aggregate  production  functions,  the  intermediate  product  transfers 
between  aggregate  activities,  and  the  applications  of  system  exogenous  and  intermediate  pro¬ 
duct  inputs  to  the  aggregate  activities).  Once  the  flow  types  have  been  modeled,  the  set  Z  of 
feasible  choices  for  the  allocations  of  resources  to  the  aggregate  activities  has  been  determined. 
Since  the  set  Z  determines  any  dependence  between  the  choices  for  the  applications  of 
resources  between  aggregate  activities,  if  we  provide  a  systematic  approach  to  modeling  the  flow 
types  then  we  will  have  presented  a  systematic  approach  to  modeling  the  dependence  relation¬ 
ships. 

To  develop  models  of  the  flow  types  associated  with  the  subnetworks  to  be  analyzed  in 
this  Chapter,  we  first  delimit  the  class  of  aggregate  production  networks  that  we  will  analyze 
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(Section  4.1).  The  definition  of  this  class  will  allow  us  to  define  aa  induced  operating  intensity  for 
each  aggregate  activity.  The  induced  operating  intensity  will  be  a  weighted  average  of  the  indivi¬ 
dual  activity  operating  intensities  and  will  have  the  interpretation  of  measuring  the  'progress*  of 
an  aggregate  activity. 

In  Section  4.2,  we  present  two  classes  of  subnetworks  for  which  we  develop  models  of  the 
flow  types  associated  with  the  aggregate  activities  within  the  subnetwork.  In  a  manner  similar  to 
how  the  flow  types  were  modeled  at  the  detailed  level,  the  flow  types  for  each  aggregate  activity 
will  be  determined  by  the  induced  operating  intensity.  By  an  analysis  of  the  inventory  balance 
constraints  associated  with  the  flow  types  at  the  detailed  level,  the  flow  types  so  modeled  will  be 
shown  to  be  consistent,  i.e.,  they  satisfy  the  appropriate  inventory  balance  constraints,  and  rea¬ 
sonable,  i.e.,  the  set  Z  induced  from  the  models  of  the  flow  types  contains  collections  of  alloca¬ 
tions  to  the  aggregate  activities  which  insure  that  there  are  schedules  for  the  detailed  activities 
consistent  with  the  allocations  which  are  or  time  and  within  budget. 

The  models  developed  are  determined  from  the  induced  operating  intensity.  Hence,  they 
are,  in  eflect,  dependent  on  the  knowledge  of  the  schedules  for  the  underlying  activities.  To 
allow  the  models  to  be  I’t^endent  on  such  knowledge  would  dear'y  defeat  die  purpose  of  aggre¬ 
gation.  Thus,  as  an  absolute  necessity,  the  models  constructed  must  be  independent  of  any 
knowledge  of  the  schedules  for  the  underlying  activities  and  therefore  cannot  be  determined 
froin  the  induced  operating  intensity. 

In  Section  4.3,  we  abstract  from  models  of  the  flow  types  for  the  subnetwork*  analyzed  in 
Section  4.2  to  obtain  models  for  flow  types  which  are  independent  This  is  done  by  modeling 
the  domains  of  the  induced  operating  intensities.  Functions  belonging  to  these  domains  are 
referred  to  as  aggregate  operating  intensities  and  are  not  necessarily  induced  from  a  schedule  for 
the  underlying  activities.  Hence,  the  models  of  the  flow  types  given  in  this  section  are 
‘independent”  from  the  underlying  network.  It  will  be  shown  that  these  models  are  'reason¬ 
able.*  The  definitions  of  the  models  given  in  Section  4J  serve  to  motivate  the  general 
definition  for  the  aggregate  network. 
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In  Section  4.4,  we  introduce  the  technique  of  replication  of  detailed  activities.  This  tech¬ 
nique  will  enable  us  to  identify  classes  of  subnetworks  which  are,  in  fact,  equivalent  to  subnet¬ 
works  already  analyzed.  Two  additional  classes  of  subnetworks  are  analyzed  through  this  tech¬ 
nique. 

In  Section  4.5,  we  present  some  concluding  remarks.  We  show  how  the  models  of  the 
flow  types  can  be  approximated  so  that  the  constraints  which  define  the  set  Z  (of  feasible  allo¬ 
cations  to  the  aggregate  activities)  are  linear.  We  then  point  out  how  a  Linear  Program  could  be 
formulated  to  accomplish  multi-project  resource-use  planning.  Next,  we  give  an  example  which 
shows  that  the  models  developed  in  Section  4.3  need  not  be  "consistent."  The  example  points 
out  the  need  to  further  restrict  the  domain  of  the  aggregate  operating  intensities.  We  also  point 
out  how  the  analyses  carried  out  for  the  specific  classes  of  subnetworks  treated  in  this  chapter 
can  be  extended  for  wider  classes  of  networks.  Finally,  suggestions  for  future  research  are  pro¬ 


vided. 
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4.1.  The  Operating  Intensity  of  an  Aggregate  Activity 

In  this  section,  we  delimit  the  class  of  aggregate  production  networks  G'  that  we  will 
analyze.  Within  this  class,  we  can  define  an  operating  intensiy  for  each  aggregate  activity.  The 
operating  intensity  will  be  used  to  define  the  domains  of  the  applications  of  exogenous  resource 
flow  types  (the  yl”s).  In  a  manner  similar  to  Section  2J.4,  we  will  also  use  the  operating  inten¬ 
sity  to  define  the  domains  of  the  other  flow  types. 

Notation 

To  differentiate  between  activities  at  the  detailed  level  from  those  at  the  aggregate  level, 
we  use  lower  case  letters  (possibly  with  subscripts)  to  denote  detailed  activities  and  upper  case 
letters  (possibly  with  subscripts)  to  denote  aggregate  activities.  The  symbol  which  denotes  an 
aggregate  activity  will  also  be  used  to  denote  the  set  of  detailed  activities  within  the  aggregate 
activity. 

The  class  of  aggregate  production  networks  that  we  will  analyze  must  first  satisfy  the  fol¬ 
lowing  property. 


Property  I 

For  each  aggregate  activity  A ,  there  ,'xist  numbers  aAh  I  €  A  such  that  if  £  b*,  >  0  for 

H.A 

some  k  then 


<*AI 


hi 

I  bu  ’ 

HA 


1  <*<». 


(4.2) 


Property  I  insists  that  each  detailed  activity  within  an  aggregate  activity  must  utilize  the  same 
percentage  of  the  total  amount  of  each  resource  required  to  complete  all  of  the  activities  within 
the  aggregate  activity. 

The  Application  of  System  Exogenous  Resources 
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Let  S  denote  a  feasible  schedule  of  start-time  for  the  activities  in  G  and  let  A  denote  an 
aggregate  activity  in  G'.  The  application  vector  of  system  exogenous  resources  tc  A  is  easily 
seen  to  satisfy 

yHMW.  (4-3) 

ItA  HA 

(The  last  equality  holds  by  (2.45).  The  definition  of  z? is  given  by  (2.56).) 

Definition  (4.4) 

The  induced  operating  intensity  for  A  derived  from  S,  denoted  by  2a,  is  given  by 

z|--  X«M/Z/S- 

It  A 


If  we  let 


"  21  bu 

it  A 


then  by  4.3-4. 5  it  is  simrle  to  verify  the  identity 


(4.5) 


yj-awxi-  *4-6) 

Any  application  vector  of  system  exogenc  s  inputs  at  the  aggregate  level  satisfies  the  form 
given  in  (4.6).  This  is  precisely  the  usual  restriction  imposed  on  the  application  vector  of  sys¬ 
tem  exogenous  inputs  imposed  by  DLAAM  (see  Section  2.3.1).  The  interpretation  of  is  that 


the  f  z*dn  expresses  the  fraction  of  the  total  resources  required  to  complete  all  of  the  detailed 

activities  within  4  up  to  time  r.  In  effect,  x?  is  a  way  of  measuring  the  progress  of  A  towards 
completing  the  detailed  activities  within  it. 

Finally,  in  order  to  insure  that  the  flow  types  at  the  aggregate  level  are  not  already  deter¬ 
mined  it  is  necessary  to  insist  that  the  class  of  aggregate  production  networks  that  we  wiH 

analyze  satisfies  the  following  property. 
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Property  II 

For  each  aggregate  activity  A,  zfe  zi 

We  insist  on  Property  II  because  if  if,  the  induced  operating  intensity  derived  from  the  early- 
start  schedule,  were  equal  to  z%,  the  induced  operating  intensity  derived  from  the  late-start 
schedule,  then  clearly  all  flow  types  at  the  aggregate  level  associated  with  A  would  be  deter¬ 
mined. 
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4.2.  Developing  Models  of  the  Flow  Types  From  the  Induced  Operating  Intensity 

In  this  section,  models  of  the  flow  types  associated  with  two  classes  of  subnetworks,  the 
Parallel  A  :B  Subnetwork  and  the  Parallel  A  :BC  Subnetwork,  are  provided.  The  models  of  the 
flow  types  ere  shown  to  be  consistent  and  reasonable  as  defined  in  Section  4.0.  The  models  of 
the  flow  types  will  be  determined  by  the  induced  operating  intensity  which  also  determines  the 
application  vector  of  system  exogenous  inputs  to  an  aggregate  activity. 

4.2.1.  The  Parallel  A  B  Subnetwork 

The  Parallel  A:B  Subnetwork  is  shown  in  Figure  (4-2).  Here,  detailed  activities 

auai . aM  are  aggregated  into  aggregate  A  and  detailed  activities  bubj,  .  .  . ,  bM  are 

aggregated  into  aggregate  B.  In  the  aggregate  subnetwork,  there  would  be  one  arc  from  A  to 
B. 

In  order  to  present  our  models,  we  first  make  two  definitions. 

Definition  (4.7) 

Let  A  denote  an  aggregate  activity  in  G'.  The  set  of  induced  operating  intensities  for  A 
derived  from  the  set  of  feasible  schedule  of  start-times  for  the  detailed  activities  within  A ,  denoted  by 
ZA  v  is  defined  by 

ZA  m  1  *a  *  L?  (m)  I  3  a  feasible  schedule  S  such  that  zA  *»  ]£  aAtz?  . 

1  1(4 

Definition  (4.8) 

The  Parallel  A  B  intermediate  product  transfer  functional  is  a  map 
fAt '•  Z«,  *  •  "  *  Z.w  —  ZB  defined  oy 

M 

fe-.t  k«,.  . .  •  2  • 


99 


See  (2.45)  in  Section  (2.3.4)  for  the  definition  c.~  the  function  4 

The  models  of  the  flow  types  associated  with  the  Parallel  A  :B  Subnetwork  are  determined 
by  the  induced  operating  intensity  in  the  following  manner.  If,  for  i  ^  k  <  n , 

M 

ykA  “  e„zA  where  zA  - 


ya  -  a,j:B  where  zB  -  £  “  ***», 

/—I 

then  the  flow  types  FA(yA,  WA),  VAB,  and  WB  associated  with  the  Parallel  A  :B  Subnetwork  are 
modeled  by 

1-1 

-low,:.,  (f  y  2.46) 

•/'  1 

i . *»v) 

(by  4.8) 

<-i 

<by2.si) 

1 

■  (4.iD 

i-i  1 

JW 

•  ]£<*«»  z*.  (by  4.4) 

i-i  ' 


(4.9) 


(4.10) 
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Thus,  we  have  shown  that  the  models  of  the  flow  types  FA  {yA ,  WA  >  and  VAB  are  consistent. 

For  the  Parallel  A  :B  Subnetwork,  there  is  only  one  inventory  balance  constraint  relevant 
to  both  A  and  B.  This  is  the  constraint 


[VM-WB)diL>  0,  Vt  €  R+ . 


(4.13) 


Substituting  in  the  definitions  for  VAB  and  WB  into  (4.13)  we  obtain  that,  Vr  €  R+, 


J  l  VaB~  WB}dfl  —  J  | X  aBbl  K'k,  ~  X  aSi,2A, 
0  *0  l'-1  »-* 


•.»&  ! 

4 

m ! 


1-1  n 


(4.14) 


Constraint  (4.14)  imposes  a  restriction  on  the  choices  for  the  induced  operating  intensities  for 
A  and  B.  Since  the  choices  for  the  induced  operating  intensities  for  A  and  B  determine  the 
application?  of  system  exogenous  inputs  to  A  and  3,  constraint  (4.14)  is  the  model  of  the  depen¬ 
dence  relationship  between  A  and  B  for  the  Parallel  A  :B  Subnetwork.  Thus,  the  question  of  rea¬ 
sonableness  of  the  models  of  the  flow  types  for  the  Parallel  A  :B  Subnetwork  reduces  to  the 
question  of  how  reasonable  is  (4.14)  as  a  model  of  the  dependence  relationship  between  A  and 


At  the  detailed  level,  M  inventory  balance  constraints 


\V,,>,-*i>)dP  >  Vt€*+./-U,  .  ..M 


ffl 


2 


insure  that  the  start-times  for  b,,  — 1,2 . A/,  are  consistent  with  the  finish  times  for  a,, 

— 1,2, . . .  ,Af.  At  the  aggregate  end,  only  one  inventory  balance  constraint  (4.14)  exists  to 
constrain  the  possible  start-times  for  6,,  /- 1,2, , . .  ,M,  given  the  finish  times  for  a,, 
<  —  1,2, ...  ,M  (as  reflected  by  zA).  Hence,  with  only  one  constraint  it  will  not  be  possible  to 
model  the  dependence  relationships  exactly.  But,  for  our  models  of  the  flow  types,  it  is 
immediate  by  (4.14)  that,  for  a  fixed  zA ,  all  zB's  whose  start-times  for  the  b,'s  are  consistent 
with  zA  do  satisfy  (4.14)  Thus,  our  models  of  the  flow  types  are  "reasonable*  in  the  sense 
described  in  Section  4.0. 

We  make  two  comments  about  the  subnetwork  just  analyzed.  First,  this  one  example 
illustrates  that  mf'deling  the  so-called  dependence  relationships  between  the  applications  of  sys¬ 
tem  exogenous  inputs  to  aggregate  activities  is  encompassed  by  modeling  the  flow  types  associ¬ 
ated  with  the  intermediate  product  transfers  between  aggregate  activities.  Second,  for  the  sub¬ 
network  just  analyzed,  the  production  function  FA(yA,WA)  was  not  equal  to  the  induced 
operating  intensity  zA.  This  differs  from  the  detailed  case  where  Fa_(ya,Wa)  -  z„  , 
<-1,2 . ,M. 

The  reason  why  we  cannot  model  FA  (yA ,  WA )  as  zA  is  because  zA  is  a  measurement  of  the 
rate  A  utilizes  its  resources,  not  a  measurement  of  when  the  activities  within  A  have  been 
completed.  To  illustrate  this  point,  consider  the  following  example  of  a  Parallel  A  B  Subnet¬ 
work  with  2  activities: 

(a)  the  duration  of  each  of  the  4  activities  is  one  period, 

(b)  aa,  -  atj  -  .05  ,  aAj  95, 

(c)  r,,*  1(0,1)  .  *»,  "1(2,3)  • 

Suppose  Fa  (yA,WA)  were  modeled  as  zA.  In  this  example  at  r-2,  zA  tells  us  that  A  has  com¬ 
pleted  5%  of  its  work.  As  measured  by  resource  use,  this  is  true.  However,  at  r-2,  a,  has 
finished  and  we  would  desire  to  allow  b\  to  start.  If  hj  did  start  at  r— 2,  then  the  following 
problem  would  emerge: 

2  2 

J  FA(yA,WA)dv-  f  zA<t»- 


2 

/ 1«* 
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2  2  /• 

J  >  J* 2gdn~  J  [aBt.f^  +  aab/bJdp 


-  .95. 


So,  for  no  choice  of  VAB  would  it  be  possible  to  begin  bj  at  time  r-2.  This  eliminates  a  possi¬ 
bility  for  b\  which  is  not  desirable. 


4.2.2.  The  Parallel  A  :BC  Subnetwork 

The  Parallel  A:BC  Subnetwork  is  shown  in  Figure  (4-3).  Here,  detailed  activities 

a1,a2, . . .  ,aK,  a*+i . were  aggregated  into  A,  detailed  activities  bub2, . . .  ,bK  were 

aggregated  into  aggregate  B,  and  cltC2, . . .  ,cL  were  aggregated  into  aggregate  C. 

The  flow  types  associated  with  the  Parallel  A -BC  Subnetwork  a n  FA{yA,WA),  VAB-  VAC, 
WB,  and  Wc.  The  inventory  balance  constraints  associated  with  these  flow  types  are: 

T  * 

+  Vr€J?+,  (4.15) 


[VAB-WB)dp>0,  Vr  €*+, 


(4.16) 


(K«r~  WcUn  >  0,  Vt  €*+. 


(4.17) 


Our  analysis  of  the  Parallel  A  B  Subnetwork  motivates  the  following  structure  for  the  models  of 
the  flow  types: 

For  some  i\  $  dj€  A+*  +  1,  if,  for  1  </c  < 


»-i 


<r  ^  j-  0 


•  .  •  r-  o  ’.V 
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yB  -  a^zg  for  zB  -  Yt<*Bb,h, 

i- 1 


yc  “  akC2C  f°r  2C  “  XaCc,zr, 

i- 1 


(4.18) 

•  (z*,)f 

(4.19) 

t-1  J 

(4.20) 

-  dxzB 

(4.21) 

Wc-diZc. 

(4.22) 

Jt  is  easy  to  venfy  that  the  models  of  the  flow  types  defined  by  (4.18-4.22)  satisfy  (4.15-4.17) 
and  thus  are  consistent. 

Let  A i  denote  the  aggregate  activity  comprised  by  at, <12, .  .a*  “d  'et  ^2  denote  the 

aggregate  activity  comprised  of  oK+\, . . .  ,au.  If  FAj{yA),WA^,  J+ 1,2,  is  defined  by  (4.9)  and 
fAiBt  fijC  are  defined  by  (4.8)  then  (4.18-4.20)  become 


Fa  -  ^FAl(yAt.WAx)  + 

VAB  “  d\fAxB (r«j, . . .  »$•*)  ““  d\VA^.g 
Vac  m  •  •  •  ■  m  ^iVajC 


(4.23) 

(4.24) 

(4.25) 


What  4.23-4.25  show  is  that  the  Parallel  A  :BC  Subnetwork  (K  +  L-M)  is,  in  effect,  a  combi- 
nation  of  a  Parallel  A\.B  Subnetwork  and  an  Parallel  Ay.C  Subnetwork.  What  4.23-4.25  also 


v  o 


,  ai*  .  V_  »  -w  •»  >  m  #  .■»  «'•  .*•  •%  .%  JN  /*  »  *•’*>  •**  4“*  **>  .**  *  *  ** 

••  -*  >*  *  *•  .  “  ••  **  *•  -  **  *•  •*  «  .  *'  v*  ~ <  af*  «*  #•.  -  «y  .  v*  .*,  if,  r.  »  .  »  .  «  , «  ,  «  ,  < 
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shows  is  that  it  is  not  possible  to  arrive  at  independent  models  for  the  flow  types  VAB  and  VAC. 
This  is  because  VAB  (VAC)  is  determined  only  from  the  zA]  (zA:)  part  of  zA.  Hence,  we  will 
model  this  sum  VA  -  VAB  +  VAC.  It  will  be  convenient  for  later  purposes  to  introduce  the  fol¬ 
lowing  definition. 

Definition  (4.26) 

The  Parallel  A  :BC  intermediate  product  Junctional  is  a  map 
/ A  BC '■  Za,  *  ' '  ‘  x  —  d\ZB  +  djZc  "• 


d\zB  +  djZc  I  zB  €  ZB ,  zc  €  Zc 


defined  by 

. *aA J  “  dyfA^B(.ZCy  .  .  .  >r8jt)  "f  d^f a^q(za^ j.  .  .  . 

Of  -  jurse,  VA  -fA:BC(z. . *.*)• 

The  inventory  balance  constraint  associated  with  VA,  WB,  and  Wc  is 

r 

C 

0<  J  { VA-(WB+Wc))dti,  Vt € /?,. 

which  reduces  to 


l\ 


°<  J  Mi - +  . 


dp..  (4.27) 


I  .a  zA  s  ZA.  It  is  immediate  by  (4.27)  that  all  start-times  for  bx . bK,  C\, ...  ,cL  which 

are  consistent  with  the  finish  times  for  a\, . . .  ,aM  satisfy  (4.27).  Hence,  our  models  of  the 
flow  types  are  ‘reasonable’  in  the  sense  described  in  Section  4.0.  The  appropriate  choice  for  dx 
and  </2  is  deferred  until  the  next  section. 


4.3.  Constructing  Independent  Models  of  the  Flow  Types 

The  models  of  the  flow  types  in  the  previous  section  were  determined  by  the  induced 
operating  intensity  and  hence  were  not  independent.  To  construct  independent  models,  we  first 
model  the  domain  of  the  induced  operating  intensity.  Functions  belonging  to  this  domain  will 
be  referred  to  as  aggregate  operating  intensities  The  models  of  the  flow  types  provided  in  this 
section  are  determined  from  the  aggregate  operating  intensities. 


4.3.1.  Constructing  the  Domain  of  the  Aggregate  Operating  Intensity 


Lw.  A  represent  an  aggregate  activity.  As  notation,  let 


Ea 


min £, ,  I,. -maxi,,  W. 

UA  A  HA  A 


\Ea,La). 


(4.28) 


(£,  is  the  early-start  time  for  activity  /  and  L,  is  the  late-start  time  for  activity  /.)  It  is  immedi¬ 
ate  by  the  definition  of  the  induced  operating  intensity  (4.7)  that  if  zA  €  ZA  then  zA  would 
satisfy  the  following  boundary  conditions: 

zA  is  a  step— function  (4.29) 


zAdp  < 


h 

ZAdfl  -  1 


Vr€*V 


(4.30) 


(4.31) 


Let  Da  denote  the  set  of  all  functions  hi  I*  (m)  satisfying  4.29-4.31.  It  is  immediate  by  the 
definition  of  ZA  that  DA  D  ZA.  DA  will  be  taken  as  the  model  of  ZA.  Functions  belonging  to 
Da  are  referred  to  as  aggregate  operating  intensities  for  aggregate  A.  Clearly,  the  functions 
belonging  to  DA  are  independent  (in  the  sense  described  in  Section  4.0). 
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4.3.2.  Constructing  an  Independent  Model  of  the  Application  Vector  of  System  Exogenous 
Inputs 

In  Section  4.2,  yA ,  the  application  vector  of  system  exogenous  inputs  to  aggregate  A ,  was 
determined  by  an  induced  operating  intensity  for  aggregate  A,  zA  €ZA.  We  will  model  the 
application  vector  yA  in  a  similar  manner.  That  is,  if  yA  denotes  an  application  vector  of  sys¬ 
tem  exogenous  inputs  to  aggregate  A  then  it  is  assumed  that,  for  1  <  k  <  n , 

for  some  zA  iDA.  (4.32) 


4.3.3.  Constructing  an  Independent  Model  of  VAB  tot  the  Parallel  A  B  Subnetwork 

For  the  Parallel  A  :B  Subnetwork,  the  intermediate  product  transfer  variable  VAB  was 
determined  from  the  induced  operating  intensity  via  the  intermediate  product  transfer  func¬ 
tional  • 

Ta  b'-  x  ’  "  *  —  Z, . 


(see  (4.8)  for  the  definition  of  fAB).  In  this  section,  we  will  construct  an  independent  model 
for  VAt  which  we  will  denote  by  VAB.  The  essential  idea  is  to  construct  an  /ab-Da  ~~  &b 
which  reasonably  approximates  fAa  on  Z#1  x  •  •  •  x  Zty  and  then  define  VAB  —  fAj(zA)  for 
*a*Da. 


An  inspection  of  the  definition  for  the  intermediate  product  transfer  functional  fAB 
shows  that  it  satisfies  the  following  boundary  conditions,  V(rJ(, . . .  ,z,y)  €  Z,,  x  •  x  Z,y, 
and  Vr  €  WB: 


Lb 

. 


-  1 


(4.33) 


f  -  f  '  f 

J /a  b^.  • .  •  .ziy)dit  <  l  f A'-B^Zfyt  ...  »^)^^  ^  f  a  3  | . zfy  )dp  (4.34) 
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It  is  also  easy  tc  see  that 


/a  Mx . z°u )  “  zk 


(4.35) 


fA  :B  ^ax . ZQy  >  "  ZS • 


(4.36) 


Our  goal  is  to  construct  an  JA:B  on  DA  which  reasonably  approximates  fAB  on 
Zax  x  ...  •  x  Z„  .  In  view  of  the  definition  of  DA  (4.29-4.31)  and  the  boundary  conditions 
4.33-4.36,  it  seems  reasonable  to  insist  that  fAB  satisfies  the  following  boundary  conditions,  V 
zA  €  Da  ,  and  Vr  €  WB : 


•A)dlt  -  1 


(4.37) 


T 


:B(zA)d».  <  I  fAB(z* . z‘)dn  (4.38) 


which  by  4.35-4.36  is  equivalent  to 


t  T 

jf  zkdn  <  Jf/; 


B(zB)dfi  <  I  z§dn 


fA.tbf)-fA,  f(zf)  -  :f 


-  zk  ■ 


(4.39) 


Fix  tA  €  Da  .  The  only  real  information  that  we  know  about  zA  is  given  in  the  boundary 
conditions  4.29-4.31  and  this  is  not  much.  Since  fA  B(xA)  satisfies  4.37-4.39,  an  intuitively 
appealing  idea  for  the  definition  of  fA  B(zA)  is  that  fA  B(zA)  should  satisfy  the  following  pro- 


For  some  p:  WB  — •  WA,  fAB(zA )  must  satisfy  the  following  equation: 


B(zA)-zb\dn 


pit) 


r  pit) 

r  r  {.t/- 

‘l  ra 


,  Vr  6  W/f  - 


(4.40) 


A  map  p  is  required  since  WB*WA.  Of  course,  to  transform  the  idea  of  (4.40)  into  a  well- 
defined  mathematically  correct  definition  for  /J:B(zA)  requires  us  to  impose  certain  restrictions 
on  p  and  the  right-hand  side  in  (4.40) 

To  arrive  at  the  restrictions,  it  will  be  convenient  to  introduce  the  following  definition.1 


Definition  (4.41) 

The  relative  progress  Junctional  for  an  aggregate  activity  A  is  a  map  pA:  DA  x  /?+  “*  10*1) 
defined  by 


Pa  Of a  ,t)  - 


I.'- 


0  if  tiWA 

zhUn 


if / €  WAt*A*iS 


1  if  r€  WA  ,zA-z$ 


1  First  introduced  in  Leachmtn,  Boyeen  [I9S21. 


pA  (zA  ,t)  is  referred  to  as  ’the  relative  progress’  of  zA  at  time  t. 

Let  us  re-write  (4.40)  as 

I  zbdu  +  pA  (zaj>(t))  Jf  [z§-zjf)dft ,  Vr  €  I 'Vg.  U* 


2) 


Ill 


Equation  (4.42)  is  an  implicit  definition  for  fA:B(zA).  Since  indefinite  integrals  of  integrable 
functions  are  absolutely  continuous  and  hence  differentiable2,  if  p  were  differentiable  then  one 
could  differentiate  both  sides  of  (4.42)  with  respect  to  r  to  obtain  a  definition  for  fAB  (2^). 
The  criterion  of  differentiability  on  p  is  not  too  restrictive.  Equation  (4.42)  and  the  definition 
for  pA  (zA  „o(r))  suggest  that  it  is  reasonable  to  insist  that  p  should  be 

continuous  and  increasing  (4.43) 


and  satisfy  the  boundary  conditions 

I  (/T,)  -  Ea  ,  p(LB)  -La  (4.44) 

If  so,  then  p  would  be  automatically  differentiable.3  Hence,  to  complete  the  definition  for 
/A:b(x, 4)  *s  given  in  (4.42)  and  hence  V'AB  it  is  sufficient  to  select  a  p  which  satisfies  4.43-4.44. 
A  natural  choice  for  p  which  does  satisfy  4.43-4.44  is  to  insist  that 


p(r)~E, 1  __  t-Eb 
Ea  —Ea  Lb  —  Eb 


Vr€  WB 


which  re-written  becomes 

p(r)  •  Ea  +  Vt£Wb.  (4.45) 

However,  this  will  not  be  a  good  choice  for  p.  To  motivate  why,  consider  the  special  case  when 
LA  -Ea  —Lb-Eb.  This  would  occur,  for  example,  if  the  detailed  activities  within  A  and  B 
had  equal  duratidns.  Under  this  special  case,  (4.45)  reduces  to  the  requirement  that 

p(r)  —  Ea  +  ( t-Eb )  -  r-  (Eb-Ea)  ,  Vt  €  WB. 


1  tlmott  everywhere. 
5  See  Royden  p.  96. 
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Hence,  p  would  be  a  simple  time  leg  of  length  EB-  EA.  Since  the  definition  for  p  as  suggested 
in  (4.45)  does  not  in  any  way  consider  the  boundaries  zA,  zA.  zj,  and  z|,  it  follows  that  under 
the  special  case  cited  above,  p  would  always  be  a  time  lag  length  EB  —  EA  regardless  of  the 
boundaries.  In  general,  the  boundaries  zjj  and  2g  are  not  determined  from  zA  and  zA  by  a  sim¬ 
ple  time  lag.  Thus,  (4.45)  should  not  be  a  good  choice  for  p. 

To  illustrate  further  why  the  choice  for  p  as  given  by  (4.45)  is  not  a  good  one  and  to 
motivate  our  choice  for  p,  consider  the  example  of  a  2-Parallel  A  :B  Subnetwork  shown  in  Fig¬ 
ure  (4-4).  The  dotted  line  in  Figure  (4-4a)  represents  the  curve  Jf  zAdp.  generated  from  a  par¬ 
ticular  zaSDa.  The  dotted  line  in  Figure  (4-4b)  represents  the  curve  Jj  fABizA )dp  as  defined 

,  ** 

by  (4.42)  with  p  defined  by  (4.45). 

When  we  modeled  VAB  in  Section  4.2.1,  we  motivated  that  it  is  through  the  inventory  bal¬ 
ance  constraint 

>  0,  Vt  €  WB 

that  the  dependence  relationship  between  the  applications  of  system  exogenous  inputs  to  A  and 
B  is  modeled.  Since  we  have  modeled  VAB  by  VAB  it  is  now  through  the  inventory  balance 
constraint 

W*(*4>-*ty»>0,  Vr€^,  (4.46) 

that  the  dependence  relationship  between  A  and  B  is  modeled.  Hence,  for  each  r  €  WBs  the 
restriction  imposed  on  the  choice  for  zB  given  a  choice  for  zA  may  be  measured  by  the  shaded 
area  shown  in  Figure  (4-4b).  In  our  example,  for  any  r  €  WB,  the  shaded  area  evaluated  at  r 
in  proportion  to  the  total  area  evaluated  at  r  is  quite  small.  This  reflects  the  fact  that  by  our 
choice  of  fA  g(zA),  i.e.,  VAB  constraint  (4.46)  is  forcing  z«  to  essentially  *run  late.*  On  the 


0 


5 

0 

1 


(■) 


Duration 
Ear ly-a tart 
Latt-atart 

Eaaourea  uaa  o  •  .2 
1 


FIGURE  (4-4) 

EXAMPLE  OF  INTERMEDIATE  PRODUCT  TRANSFER  CURVE 
DERIVED  FROM  A  SIMPLE  TIME  LAG 
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other  hand,  an  inspection  of  Figure  (4-4a)  reveals  that  the  shaded  area  evaluated  at  p(r)  in 
proportion  to  the  total  area  evaluated  at  p(r)  is  fairly  large.  This  reflects  the  fact  that  zA  is 
essentially  "running  early."  So  what  we  have  is  that  zA  is  "running  early"  but  we  are  constraining 
zb  to  run  "late". 

Based  on  this  one  example,  let  us  insist  that  regardless  of  the  choice  for  p,  fA.n(zA)  must 
satisfy  the  additional  property  that  the  shaded  areas  normalized  by  the  total  areas,  evaluated  at 
the  corresponding  points  in  time,  are  equal.  Mathematically,  we  are  insisting  that  fAB (zA) 
satisfies  the  following  criterion: 


r  p(r> 

J  j [fA  B-zftdp  dx  J  J{zA-zj]dpdx 

T  p(r) 

£  J =&)dp  dx  J  J {zj  -  zfldp  dx 


Vr  €  WB . 


(4.47) 


It  will  be  convenient  to  introduce  the  following  notation:  for  an  aggregate  activity  * , 


X 

^  Of)  -  J  {rf-rjJt/M ,  x  €  /?+ 


(4.48) 


(Note  that  hA  (x)  -  hA  (LA ) ,  Vr  >  LA  .)  , 
Re-write  (4.47),  Vr  €  WB,  as 


r  t 

| a) dp  dx  -  ^  f  zfjdp 


dx  + 


r  p(r) 

J  ha(x)dx  £pA(z, 


■,Apc)hA(x)dx 


p(r) 

f 

i. 


a  »  / 

f  hA(jc)dx 


(4.49) 


Proposition  (4.3.1) 


-4 


is 


r« 


-4 

ess! 


ft'/ 
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If  Vz,,  6  Da  ,  j\  b  (zA )  is  required  to  satisfy  (4.49)  for  some  p  satisfying  4.43-4.44  then 
satisfies  (4.42)  Vz^  €  A<  for  the  same  p  if  and  oniy  if  p  is  the  (unique)  solution  to  the 

equation 


t  pit) 

hB(x)dx  ^  hA(x)dx 

fig  ( Lg ) 


MA«> 


Vr  6  Iffl . 


(4.50) 


Proof  of  Proposition  (4.3.1) 

Suppose  Vz^  €  A<  and  for  some  p  satisfying  4.43-4.44,  fAj(zA)  satisfies  bcth  (4.42)  and 
(4.49).  Differentiate  each  side  of  the  equation  (4.49)  with  respect  to  r  to  obtain  the  equation,  V 

•tCH^Vz^A,, 


/ 

J -  J zidp  +  Ai  (z^  ,{p(r))hA  (p(r))p'(r)  * 


hB(x)dx 


I- 


<*)dx 


p(r) 

I" 


(z4^r)A^(x)<£r 


p(r) 

•l 


hB(,r)  I  hA(x)dx—  fiA(p(r))p'(r)  I  hg(x)dx 


p(r) 

7 - 

1  hA(x)dx 

(4.51) 


Since  fAB(zA)  also  satisfies  (4.42),  if  we  subtract  (4.42)  from  (4.51)  and  rearrange  terms  we 
obtain  the  equation,  Vr  €  WB,  Vz*  €  , 


r 

p(r) 

1  hg(x)dx 

.  /_  /. /-\\  *4 

natTj  —  nA  tpvrjjp  tx;  '(r)  . . 

.  /• 

^(r) 

J  Aj(x)<ftc 

J  hA  (x)dx 

,  T. 

% 

(4.52) 
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4 


Since  hA  Cx)  >  0  (except  at  x  -  EA )  then  it  is  clear  that 


pA(zAp(r))  - 


pit) 

J  pA(zApc)hA 
J  hA(jc)dx 


*0,  Vr  €  Wa 


(4.53) 


pA  (iAj)  is  strictly  increasing.  Qearly,  DA  contains  zA  's  such  that  pA  (zA  jc)  is  strictly  increasing 
(by  Property  I,  Section  4.1).  Since  (4.52)  holds  Vr^  €  DA  and  since  the  expression 


t 

{ 

f>B  (*>  ~  K  <p(t))p'(t)  ~ 

t 


hB(x)dx 


hA  (jc)dx 


is  independent  of  zA ,  it  follows  by  (4.53)  that 


or,  equivalently, 


0  —  ha  (r)  —  hA  (p(r))p'(r) 


l"’ 

b 


hB(x)dx 


( x)dx 


MT>  hA(p(r))p'(r)  „  _ 

- — - »  Vr  €  Wg. 


hB(x)dx 


/- 


(x  )dx 


hg(r) 


hg(jc)dx 


f 

•  ~  In  |  hgbc)dx  ,  Vr  €  WB 


(4.54) 


*.*  '  V\  \  1  *  •\>  ‘.r  \* 


V*'  v  «**  •.*  V  •.*  ,*  *  %*.  V  •••  •/  o .« 
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Differentiate  each  side  of  (4.50)  with  respect  to  r,  we  obtain  the  identity 


hB  (*)  **  T~ 77^T  hA  (p(r))p'(r) ,  Vr  €  . 


(4.58) 


Substituting  (4.58)  into  (4.57)  we  have  that 


T  T  T 

fA  B(2A)dfi  -  zfaii  +  pA(zAxp(r)) ^  [z§ -  zj;}dii  t  Vr  €  JFj 


which  is  equation  (4.42).  This  concludes  the  proof.  ■ 

Our  model  of  V\g  is  now  complete.  We  define  fAj(zA)  implicitly  by  equation  (4.42) 
with  p  chosen  to  satisfy  implicitly  the  equation  (4.50)  and  then  set  Vga-ft-jUg).  As  a  result 
of  this  choice  for  fA:B(zA),  we  know  by  Proposition  (4.3.1)  that  it  also  satisfies  (4.49). 

We  make  one  final  comment.  Consider  the  restriction  imposed  on  the  choice  for  zB  given 
a  particular  choice  for  zA.  As  we  have  stated  before,  this  restriction  is  reflected  in  the  inven¬ 
tory  balance  constraint 


J.***lK 


y;adp,  Vr  €  Wg . 


(4.59) 


Substituting  in  (4.59)  the  definitioh  for  VAB  we  have  that,  Vr  €  WB, 


T  T  .  t  t  T 

J tad*  <  J  VAgdfi">  ^  fA  j(zA)dp.—  ^ zidp  +  pA(zAj>(r)) ^  (z# -  zfrdp 


or,  after  rearranging  terms, 

Pabas)  <  pAizAj>(r)). 

(4.60) 

Essentially,  Leachman  and  Boysen’s  method  for  modeling  the  dependence  relationship  between 
A  and  B  for  a  Parallel  A  :B  Subnetwork  was  to  require  that  the  chokes  for  zA  and  zB  satisfy 
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(4.60). 


4.3.4.  Constructing  an  Independent  Model  of  VA  far  the  Parallel  A  :BC  Subnetwork 

For  the  Parallel  A  BC  Subnetwork,  the  intermediate  product  transfer  variable  VA  was 
determined  from  the  induced  operating  intensity  via  the  intermediate  product  transfer  func¬ 
tional 

/a-.bc:  Z.xx  •  •  x  Z.u  -*  diZg  +  diZc 


(see  4.26  for  the  definition  of  fA:BC )•  In  this  section,  we  will  construct  an  independent  model 
of  VA  which  we  denote  by  VA.  The  essential  idea  is  to  construct  an  fA.ac'-  Da  —  dyDB  +  d2Dc 
which  reasonably  approximates  fA  K  on  ZB,  x  •  •  •  x  ZBy  and  then  define  VA»fABC(zA)  for 


2 a  €  Da  . 

An  inspection  of  the  definition  for  the  intermediate  product  transfer  functional  fA  :BC 
shows  that  it  satisfies  the  following  boundary  conditions,  V(zBl, . . .  ,z„v)  €  Zfl)  x  •  •  •  x  Z„v, 
and  min (£«,£<:)  <  t  4?  ma x{LB,Lc), 


1 


nuiZi^) 


fA*c(*ai . "  1 


(4.61) 


t  T 

I  •  •  •  •2iiu^dtt’  ^  I  /aJK^a j»  •  •  • 

nifiU|^*)  owfc|^) 


<  I  /a  :BC  (^/j  ■  •  •  •  >2fy)dp 

ma{It,£c) 


(4.62) 


For  notational  convenience,  let  £jc*min(£0,£c)  snd  LK  -  max (LttLc)  and 
WK  m  Lfac.Iac).  U  is  easily  verified  by  the  definition  of  fA  ,Bc  the* 


/rfjcfc., . m  d\2i+dizt  (4.63) 
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f a  sc  (*/, . if)  -  dtzf+  d2z£ 


Of  course,  boundary  conditions  4.61*4.62.  are  exactly  analogous  to  the  boundary  condi¬ 
tions  for  f A  g  given  in  4.33*4.34.  Based  on  the  motivation  and  explanation  provided  for  the 
development  of  the  definition  for  fA.g,  we  define  fA.sc  implicitly  by  the  equation,  Vr  €  WBC> 


r  r  r 

f Sa  bc^a  )dn  -  r /a  bc (2A)dfi  +  pA  (zaj(t))  J  \fA  Kiz^)  ~  f A'.SC  (zA)}dfi 

*§C 

t  t 

-  £  [diz^+djz^dfi  +  pA  (zA  [di(z§-zi)  +  d2(z£-zk))dn 


for  some  p:  “•  WA  which  satisfies 


p(.Egc)  “  £* ,  pU-ac)  "  £4 


and  which  is 


continuous  and  increasing. 


Analogous  to  the  additional  requirement  imposed  on  fA.g  (see  4.49),  let  us  insist  that,  V 


zA  €  Da  ,  fA  BC  also  satisfies  the  equation,  Vr  €  W/ac, 


*(r)  r 

izA-zi)dndx  {/AMrbA)-(dxzi+d22$)]dpdx 

p(r)  r  * 

By  a  proof  similar  to  the  proof  of  Proposition  (4.3.1)  it  may  be  readily  verified  that  p  must 


satisfy  the  equation 


I  [d\hB  Or)  +  d2hc  (x))dx  I  hA  (jc)dx 
* _  ..  „  b  . 

dihA(Lg)  +  dJhc(JLc)  MA(LA) 


l- 


(4.64) 


121 


To  complete  the  moJel  of  V\  requires  us  to  select  the  weights  dx  and  d2.  To  motivate 
our  choice  for  </,  and  d2,  let  us  first  consider  the  restriction  imposed  on  the  choices  for  z,  and 
zc  given  a  particular  choice  for  zA .  As  we  have  stated  before,  this  restriction  is  reflected  in  the 
inventory  balance  constraint 


l 


{d\2g  +  d2zc)dfi  4 


f  v'Adnr\  /; 

£  £ 


:BC^A)dn,  Vt  €  • 


(4.65) 


Substituting  into  (4.65)  the  definition  for 


)4p  we  obtain. 


after  re-ananging  terms. 


the  inequality 


T 

J  {</,tZa-z|)  +  </2(zt-2£M<fM 

** . — - *Pa(ZaJ>(  r)>,  Vt€»W. 


(4.66) 


/ 


Irfj (z| -  zj?)  +  d2(.z£  zfc)}d/i 


(4.66)  is  equivalent  to  the  inequality 


„  ^  *M)  *  l" <I-r> 


(4.67) 


Let  <4,  denote  the  set  of  activities  in  A  which  precede  the  activities  in  B,  let  A2  denote 
the  set  of  activities  in  A  which  precede  activities  in  C.  and  define 


I  <*AI 

If  A,  .  m 

<*A,A  -  j-  1.2- 


ItA 


(4.68) 
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It  may  be  easily  derived  algebraically  that  if  zA  -  z|  for  some  feasible  schedule  S,  then 


Pa(za4>W  "  Oa,a 


hAMr)) 


Pa^^M) 


<*AiA 


*aMt)) 


«  A ,  A  f>A  ,  (p  (r  )  )  +  a  A  r*  *A  1  (p  (t  )  > 


/Uj(z|,.pfr))-  (4.69) 


To  motivate  our  choice  for  </t  and  d2,  consider  the  special  case  when,  for  some  /  >  0, 


z/( t)  -  zj^  (r  —  /)  -  rjj  (r  -  /)  -  z*(t) ,  Vr  €  IF*; 


(4.70) 


ri(r)  -  rij(r-/)  -  7^(7-/)  -  z^(t),  VrCH^. 

In  this  case,  an  inspection  of  (4.64)  shows  that  p(t)-t-/  (which  is  intuitively  clear).  Furth¬ 
ermore,  p(t)  - 1  - 1  and  (4.70)  imply  that 


hB(r)  -  A^tpfr))  •  h4j(p(r))  -  hc(r),  Vr  €  IF*-. 


(4.71) 


Substituting  the  identities  found  in  (4.69)  and  (4.71)  into  the  inequality  (4.67)  gives  us  the  ine- 
«f*lity,  Vr  €  WK, 

0  >  l</iPaOra,T)-a/l)/<p^1(zi)  ^>(r))|  +  [dtfC he*)"* a ^Pa^a ,*(*)))  ■  (4.72) 

Si|nce  the  subnetwork  associated  with  A\  and  B  (resp.,  A2  and  C)  comprises  a  Parallel  Ay.B 
Subnetwork  (resp.,  a  Parallel  A2:C  Subnetwork),  it  follows  from  our  work  for  that  subnetwork 
thlat  we  would  like  to  constrain  z*  and  Zc  by  the  inequalities 


P#(z*,r)<p,4|(zi|,p(T))  ,  rCW'jc 
Pc<*c s)  <  4>(r» ,  t  €  WK . 


(4.73) 
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Hence,  if  we  set 

dj  -  aAjA  ,  j  —  1,2  (4.74) 

then  any  zB  €  Db,  zc  £  Dc  which  do  satisfy  (4.73)  will  satisfy  (4.72).  Hence,  for  this  choice  of 
d j  and  d2  (and  for  this  special  case),  we  will  not  eliminate  any  reasonable  choices  for  zB  and  zc- 
Therefore,  we  specify  d\  and  d2  by  (4.74). 

We  make  one  final  remark  about  the  model  for  V\  and  hence  our  model  for  the  depen¬ 
dence  relationship  between  A,  B,  and  C.  Leachman  and  Boysen’s  method  for  modeling  the 
dependence  relationship  between  A,  B,  and  C  is  to  insist  that  zA ,  zB,  and  zc  satisfy  the  ine¬ 
quality 

aAiAPBbBJi'b)) +  aAjApc(zcJ>clM'><PAbA,T),  Vt€IVbc  (4.75) 
where,  for  i  —  B,C,  p,_1(r)  satisfies 

t  *»r'<T> 

hA(x)dx  ^  h,(jc)dx 
*aUa) 

An  inspection  of  the  relationship  between  (4.75)  and  (4.72).  reveals  that,  except  for  the 
differences  in  the  definition  of  p(r),  Leachman  and  Boy  sen  's  method  for  modeling  the  dependence 
relationship  is  our  method  under  the  special  case  cited  above.  For  more  general  cases,  we  feel  our 
approach  is  more  sensible  since  it  will  weight  a,**,  and  a AAj  by  time-varying  factors.  Ultimately, 
it  will  have  to  be  tested  to  see  if  it  performs  better.  (See  Section  (4.5.4)  for  a  discussion  con¬ 
cerning  testing  of  the  approach  for  modeling  the  flow  types.) 
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4  \5.  Constructing  an  Independent  Model  of  the  Production  Fraction 

For  both  subnetworks  analyzed,  FA  (yA ,  WA )  does  satisfy  the  usual  boundary  conditions.  It 
would  therefore  be  appropriate  to  model  it  in  exactly  the  same  manner  as  we  modeled  Kb  or 
VA.  Since  VAB  and  VA  were  modeled  in  detail,  we  omit  the  analysis. 
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4.4.  Analyzing  Subnetworks  Through  the  Technique  of  Replication 

In  Section  4.3,  we  developed  models  of  the  flow  types  for  the  Parallel  A  B  Subnetwork 
and  the  Parallel  A  :BC  Subnetwork,  This  section  introduces  the  technique  of  replication  of 
detailed  activities  and  uses  it  io  develop  models  of  the  flow  types  for  two  further  classes  of  sub¬ 
networks.  We  proceed  to  develop  models  of  the  flow  types  for  two  classes  of  Non-Parallel 
A  :BC  Subnetworks. 

4.4.1.  The  Complete  Precedence  Non-Parallel  ABC  Snbnetwork 

The  Complete  Precedence  Non-Parallel  ABC  Subnetwork  is  shown  in  Figure  (4-5).  Here, 
detailed  activities  a j, a 2. ...  ,au  were  aggregated  into  aggregate  A ,  detailed  activities 

b\,bi, . . .  ,bM  were  aggregated  into  aggregate  B,  and  detailed  activities  cuc2 . cM  were 

aggregated  into  aggregate  C.  Let  (/  denote  the  underlying  production  system  and  let  G* 
denote  the  production  network  associated  with  G.  We  will  construct  a  production  system  H 
which  is,,  in  effect,  ‘equivalent*  to  G  and  analyze  H  to  arrive  at  models  of  the  flow  types  for 
the  subnetwork  shown  in  Figure  (4-5). 

To  construct  H ,  let  us  first  construct  the  production  network  /f*  associated  with  H: 

Step  1:  Add  nodes  aJ2 . au2  to  G.  Re-label  node  a,  to  a,i,  1  - 1,2, . . . ,  M. 

Step  2:  Add  arc  (d,al2)  to  G*  if  there  is  an  arc  id, a,)  in  G,  />■  1,2 . M. 

Add  arc  (a,2,d)  to  G*  if  there  is  an  arc  ( a,,d )  in  G,  1  —  1,2, . ,M. 

H‘  is  the  network  obtained  from  G*  from  Steps  1  and  2.  The  subnetwork  in  H  which 
corresponds  to  the  Complete  Precedence  Non-Parallel  A  :BC  Subnetwork  in  G  is  shown  in  Fig¬ 
ure  (4-6).  As  notation,  we  have  let  the  symbol  Aj,J — 1,2,  denote  the  aggregate  activity  in  H 

comprised  of  activities  a, J% /— 1,2 . M.  Essentially,  //*  is  G*  with  nodes  aj, 

’replicated.* 

The  production  system  H  whose  production  network  is  H’  is  defined  as  follows: 
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(i)  The  duration  and  resource  requirements  of  activity  a,j  in  H  are  identical  to 

the  duration  and  resouice  requirements  of  activity  a,  in 
G,  i-U . A/,; -1,2. 

(ii)  For  an  activity  in  H  other  than  some  a,;,  the  duration  and  resource  require¬ 
ments  are  identical  to  the  duration  and  resource  requirements  of  its 
corresponding  counterpart  in  G. 

(iii)  It  is  required  that  the  operating  intensities  2^  and  z„2  are  equal. 

/'  — 1,2, . . .  ,M. 


Since  the  subnetwork  pictured  in  Figure  (4-6)  is  comprised  of  two  Parallel  A  B  Subnet¬ 
works,  it  follows  that  for  the  production  system  H,  the  variables  FAl(zAl),  F4j(zaJ, 

VAjc,  WB,  and  Wc  are  determined  from  the  induced  operating  intensities  zAj  and  zAj  exactly  as 
described  in  Section  4.2.1.  Given  zAj  €  ZAj,  J  — 1,2,  the  restrictions  imposed  on  the  choices  for 
Zg  and  Zf  are  modeled  by  the  constraints 

(4.76) 

(4.77) 

It  is  immediate  by  our  construction  of  H  and  (ii)  above  in  particular  that  ZA,  the  set  of 
induced  operating  intensities  for  aggregate  A  in  G,  is  equal  to  ZAjiJ- 1,2,  the  set  of  induced 
operating  intensities  for  aggregate  Aj  in  H,  j  -  1,2,  when  these  sets  are  viewed  as  sets  of  func¬ 
tions.  If  one  identifies  ZA  with  Z4,,  j  — 1,2,  then  it  is  immediate  by  our  method  for  modeling 
the  intermediate  product  transfers  presented  in  Section  4.2  that 

Kic-V-  (478) 


ji*v- 

*1 


zg)dn  >  0,  t£Wb 


jW- 

fc 


zc)dfi>  0,  r  €  Wc 


Hence,  it  is  appropriate  to  view  aggregate  A  in  G  as  producing  two  identical  products,  i.e.,  to 
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define  for  zA  €  ZA 

FAiyAM-  [FBAlyA  ,WA),FHyA,WA)}-\FAJyA,WA)s  FA2(yA ,  WA ))  (4.79) 

(where  FA  is  defined  by  (4.9))  and,  in  view  of  (4.79),  to  define  for  zA  €  ZA 

VAa  ”  /J4,:atz«,.  -  • .  >Zau)  ™  K<,a  *  VAc  “  “  ^2c  •  (4.80) 

The  restrictions  on  the  choices  for  za  and  zc  given  zA  would  be  modeled  by  constraints  4.77- 
4.78  (with  the  substitution  of  the  identities  given  in  4.81). 

Since  the  models  of  the  flow  types  determined  from  the  induced  operating  intensities  for 
these  two  subnetworks  are  identical,  it  seems  reasonable  to  insist  that  the  independent  models 
of  the  flow  types  determined  from  the  aggregate  intensities  also  be  identical.  Thus,  we  define 
v;B,  VAC,  and  FA (yA,WA)  to  be 

y'AB  -  v:x„ 

Vlc-V^c 

F'a (yA,WA)-  (f;( (yA,WA),  FAl (^ ,  )] . 

The  variables  FA  iy,,,WA),  /  - 1,2,  VAjg  and  VAjC  were  defined  in  Section  4.3.2.  This  con¬ 
cludes  our  analysis  of  this  subnetwork. 

4.4.2.  The  'Partial  Precedence  Non-Parallel  A  :BC  Subnetwork 

The  Partial  Precedence  Non-Parallel  A  :BC  Subnetwork  is  shown  in  Figure  (4-7).  Here, 

detailed  activities  ax, ...  ,au  were  aggregated  into  aggregate  A ,  detailed  activities  bx . bu 

were  aggregated  into  aggregate  B,  and  detailed  activities  c  t . cL  were  aggregated  into  aggre¬ 

gate  C.  Let  G  denote  the  underlying  production  system  and  let  G *  denote  the  production  net¬ 
work  associated  with  G.  We  anil  construct  a  production  system  H  which  is,  in  effect,  equivalent 


to  G  and  analyze  H  to  arrive  at  models  of  the  flow  types  for  the  subnetwork  shown  in  Figure 
(4-7). 

To  construct  //,  let  us  first  construct  the  production  network  H'  associated  with  H: 

Step  1:  Add  nodes  2  G'  Re-label  node  a,  to  alU 

i~K+l,  . .  ,M. 

Step  2:  Add  arc  (d,al2)  to  G*  if  there  is  an  arc  (d,a,)  in  G,  i  —  K+ 1, . . .  ,M. 

Add  arc  (a,j,d)  to  G’  if  there  is  an  arc  (a^d)  in  G,  i»K+ 1 . M. 

H’  is  the  network  obtained  from  G*  from  Steps  1  and  2.'  The  subnetwork  in  H  which 
corresponds  to  the  Partial  Precedence  Non-Parallel  A  :BC  Subnetwork  in  G  is  shown  in  Figure 
(4-8).  As  notation,  we  have  let  the  symbol  Af,  >—1,2,  denote  the  aggregate  activity  in  H 
comprised  of  activities  a,;,  1— JC+1, .  .  .  ,M.  Essentially,  H'  is  G*  with  nodes  aK+\, .  . .  ,aM 
“replicated." 

The  production  system  H  whose  production  network  is  J/*  is  defined  as  follows: 

(i)  The  duration  and  resource  requirements  of  activity  atj  in  //  are  identical  to 

the  duration  and  resource  requirements  of  activity  a,  in  G , 
/  — A+l. . . .  >— 1,2. 

(ii)  For  an  activity  in  H  other  than  some  a(>,  the  duration  and  resource  require¬ 
ments  are  identical  to  the  duration  and  resource  requirements  of  its 
corresponding  counterpart  in  G. 

(iii)  It  is  required  that  the  operating  intensities  z,a,  are  equal, 

/-*+ 1 . M. 

It  is  immediate  by  our  construction  of  H  that  ZA  ,  the  set  of  induced  operating  intensities 
for  aggregate  A  ‘  in  H,  equals  ZA ,  the  set  of  induced  operating  intensities  for  aggregate  A  in  G, 
and  that  ZAj,  the  set  of  induced  operating  intensities  for  aggregate  /<  3  in  //,  equals  ZAj,  the  set 
of  induced  operating  intensities  for  A 1  in  G  (when  all  sets  are  viewed  as  sets  of  functions).  If 
one  identifies  Z4-  with  ZA  and  ZAj  with  ZAf  then  it  is  immediate  by  our  modeling  of  the  inter¬ 
mediate  product  transfers  presented  in  Section  4.2  that,  for  zA  €  ZA, 


V 
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“  f A  '  J  (*•  | •  •  •  •  #  2* 


«+»r 


•2»m) 


VAC  “  K*,C  “  //<j:C^a(x+,)2*  •  *  • 

O'*  ,  #4  )  -  i  a***,  +  L  aCc2c,  ■ 
i-l  <-*+1 

Since  the  models  of  the  intermediate  product  transfer  flow  types  determined  from  the 
induced  operating  intensities  for  these  two  subnetworks  are  identical,  it  seems  reasonable  to 
insist  that  the  independent  models  of  these  flow  types  determined  from  the  aggregate  operating 
intensity  also  be  identical.  When  At  and  A '  have  been  aggregated,  the  subnetwork  pictured  in 
Figure  (4*8)  is  a  Parallel  A  :BC  Subnetwork.  From  having  studied  this  case  in  Section  4.3,  we 
know  that  while  it  is  not  possible  to  arrive  at  an  independent  model  for  either  VaB  or  VAjC  we 
do  know  how  to  constrain  the  choices  for  zB  and  2c  given  zA-  by  modeling  the  sum 
VA-—  VA  B  +  yA]c-  That  is,  zB  and  zc  must  satisfy  (see  4.67) 

*,.<**«>  > 


+  dj 


hc(r) 


d\hB(r)+ djhcir) 


pc(zc,r),  Vt  6  WK 


(4.81) 


where 


d\ 


HA" 


dim*AVr 


I 


HA 


*S"  -  *AA<1  +‘*AjA"Za, 
•  aAA-zf+aA}A-zfl 


134 


*a-  -  <*AA-zj[  +aAyt*}2 


™  aA  'A  "ZA  +  ®/4  "ri, 


and  p(r)  satisfies  the  equation,  Vr  € 


T  ph) 

^  [d\hB(jc)  +  d2hc(jc)]dx  £  [d\hA(x)  +  d2hAj(x))dx 


d\hB  (LB)  +  djhc  (Lc)  d\hA{LA  ) +d2hAi(LAJ 


1 


\dxhA(x)+d2hAi(x)\dx 


dxhA<JLA)  +  d2hAi(LA} 


Note  that  if  zA--zA-  for  some  feasible  schedule  S  in  H  then  by  (4.69)  we  would  have  that,  V 
t  €  WKt 


PA-(*A’J>b))  -  </i| 


hA(p(r)) 


d\hA>U>(r))  +  d3hA.(p(r)) 


\PaWj>^ 


+  df 


*Aj<p(t)) 


dihA(p(r))  +  djhAi(fi(r)) 


,  f  hA(p(r))  \  .. 

*'  d}hA(p(r))  +  d2hAl(p(r))  **'**+'*” 


+  dj\ 


A.,(p(r)) 


d\hA  (p(r))  +  d2hAJp(r)) 


Intuitively,  pB  and  Pc  are  betas  constrained  by  pA  and  jr4|  which  makes  sense  based  on  our 
analysis  of  Parallel  A  iB  Subnetworks. 
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If  one  used  zA-,  then  (4.81)  would  serve  to  model  the  dependence  relationships.  Unfor¬ 
tunately,  it  is  yA  not  yA-  that  models  resource  use  for  A .  Without  knowing  zA]  we  could  not 
easily  transform  from  yA-  to  yA .  However,  if  we  made  the  simplifying  assumption  that 

zA- “  V  -  iAl  (  -  Za  -  Za,  )  (4.82) 

then 


9a— 9a  voApa--Pa. 

Essentially,  Leachman  and  Boysen’s  method  for  modeling  the  dependence  relationships  for  this 
subnetwork  is  to  use  (4.82)  with  the  simplifying  assumption  (4.83).  While  (4.83)  is  intuitively 
appealing-  without  additional  information,  it  assumes  all  activities  within  aggregate  A  progress 
at  the  same  rate-  it  is  restrictive.  More  than  anything  else,  it  points  out  the  fact  that  modeling 
dependence  relationships  for  complicated  (i.e.,  non-parallel)  subnetworks  is  difficult. 
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4.5.  Concluding  Remarks 


4.5.1.  Extensions 

The  development  of  the  models  of  the  flow  types  for  the  Parallel  A  J?,  Parallel  A  :BC, 
Complete  Precedence  Non-Parallel  A:BC ,  and  the  Partial  Precedence  A  SC  Subnetworks 
presented  in  Sections  4. 1-4.4  may  be  repeated  for  similar  classes  of  subnetworks.  For  example, 
it  is  easy  to  see  by  'symmetry*  how  one  would  model  the  flow  types  for  the  Parallel  BCiA  Sub¬ 
network  (see  Figure  (4-9»,  the  Complete  Precedence  Non-Parallel  $C .A  Subnetwork  (see  Fig¬ 
ure  (4-10)),  and  the  Pariial  Precedence  Non-Parallel  BCA  Subnetwork.  The  subnetworks 
shown  in  Figures  (4-9),  (4-10)  and  those  analyzed  in  detail  in  this  chapter  include  all  of  the 
subnetworks  presented  in  the  leachman  and  Boysen  [1983]  paper. 


4.5.2.  Constructing  Tractable  Models 

In  order  to  develop  a  convenient  method  for  accomplishing  multi-project  resource-use 
planning,  it  is  desirable  to  construct  models  of  the  flow  types  which  are  tractable.  That  is,  the 
constraints  which  define  the  set  Z  of  feasible  choices  for  the  allocations  of  resources  to  the 
aggregate  activities  induced  from  the  models  of  the  flow  types  are  linear.  An  inspection  of  the 
models  of  the  flow  types  developed  in  Section  4.2-4.4  reveals  that  only  the  indefinite  integrals 
of  the  flow  types  and  not  the  flow  types  themselves,  are  required  to  determine  the  constraints 
which  define  the  set  Z.  To  define  models  of  the  flow  types  which  are  tractablei  we  simply  con¬ 
struct  a  piecewise-linear  approximation  to  the  indefinite  integrals  of  the  flow  types. 


Deflnitien  (4.83) 

Let  x  €  L?  (m).  The  piecewise-linear  approximation  to  J  zrfp,  denoted 

o 


3*  ‘jl  ufc*  ‘jl 

J  z'Jm*  J  ♦  Or-Ul)  J  xdn  -  J  xdp. 


by  J x'dn,  is 


defined  by 
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(where  U1  equals  the  greatest  integer  of  x). 

The  models  of  the  flow  types  are  now  understood  to  be  the  piece  wise-linear  approximations  to 
the  previous  models  of  the  flow  types.  (Since  the  aggregate  operating  intensities  w<  e  con¬ 
strained  to  be  step-functions,  the  indefinite  integrals  of  such  flow  types  are  already  piecewise- 
Iinear.) 

An  inspection  of  the  constraints  which  model  the  dependence  relationships  between 
aggregate  activities  given  in  Section  4.3  clearly  shows  that  they  now  become  linear  constraints. 
Since  all  flow  types  and  hence  all  constraints  which  model  dependence  relationships  are  deter¬ 
mined  from  the  aggregate  operating  intensities,  a  Linear  Program  which  accomplishes  multi¬ 
project  resource-use  planning  may  be  formulated  by  having  the  indefinite  integrals  of  the  aggre¬ 
gate  operating  intensities  as  the  decision  variables.  A  suitable  objective  function  (one  for  which 
minimizes  cost,  for  example)  is  all  that  is  required.  See  Leachman  and  Boysen  [1983]  for  their 
choice  for  an  objective  function. 

4.5.3.  Lack  of  Consistency  for  the  Models  of  the  Flow  Types 

The  constraints  which  model  the  dependence  relationships  are  required  to  be  satisfied  and 
appear  in  the  Linear  Program  discussed  in  Section  (4.5.1).  However,  it  is  possible  that  the 
inventory  balance  constraints  which  link  the  production  variable  to  the  intermediate  product 
transfer  variables  are  not  satisfied. 

Consider  the  example  of  a  Parallel  A  :BC  Subnetwork  (where  A  has  two  activities)  shown 
in  Figure  (4*11).  The  inventory  balance  constraint 

|F;^,^)~K;}^  >  0  (4.84) 

is  not  always  satisfied.  Furthermore,  the  indefinite  integral  of  V‘4, 


.«■-  <r.  \ 

.v.v 


./.  «•.  .\v 

■V.v . 


(4.85) 


is  not  always  increasing,  i.e.,  VA  is  negative  at  some  points  in  time. 

The  main  reasons  why  such  an  example  can  be  constructed  are  twofold.  First,  the  domain 
Da  is  too  large,  i.e.,  additional  constraints  must  be  imposed  to  insure  that  the  variable  zA  is 
close  to  an  induced  operating  intensity  (the  particular  zA  in  Figure  (4*11)  is  not).1  Second,  the 
boundaries  z|  and  r|  are  completely  different  in  "shape*  than  the  boundaries  zA  and  zA.  Perhaps 
such  inconsistencies  point  to  the  fact  that  the  underlying  activities  should  not  be  aggregated. 

In  any  event,  models  of  the  flow  types  associated  with  the  production  functions  and  the 
intermediate  product  transfers  should  satisfy  the  inventory  balance  constraints  similar  to  (4.8S). 
In  addition,  an  inconsistency  such  as  (4.86)  should  not  occur  either.  As  a  test  for  the  construc¬ 
tion  of  the  models,  these  criteria  could  be  employed. 

4.5.4.  Suggestions  for  Future  Research 

It  is  desirable  to  extend  the  analyses  provided  in  Sections  4.2-4.4  to  wider  classes  of  more 
complicated  networks.  In  addition,  the  approach  to  constructing  independent  models  (Section 
4.3)  which  is  the  key  to  modeling  the  dependence  relationships  should  be  tested.  Results  of  a 
preliminary  test  of  the  approach  may  be  found  in  Dalebout’s  {1983]  Master’s  thesis.  Based  on 
her  results,  the  approach  due  to  Lcachman  and  Boysen,  extended  in  this  chapter,  seems  to 
work  well.  Lastly,  further  study  into  methods  for  restricting  the  domains  of  the  aggregate 
operating  intensities  so  as  to  insure  consistency  is  desirable  and  might  be  the  key  to  more  real¬ 
istic  models  of  the  dependence  relationships. 


1  Set  LMchman  md  Boywn  1 191}  I  for  iMt  proponit  for  funar  rawricttoru  on  Um  domain  D4. 
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